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A theoretical  analysis  of  the  combustion  of  composite  propellants  (with 
an  emphasis  on  the  heterogeneous  effects)  is  presented.  The  analysis  based 
on  the  combination  of  a unique  statistical  treatment  of  the  burning  surface 


(Sf-Tfr)  with  a very  comprehensive  multiple  flame  type  of  physiochemi cal  model: 
for  the  combustion  process.  | 

The  analysis  is  divided  into  four  parts:  the  steady  state  burning  rateJ 
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erosive  burning  rdte,  the  pressure  coupled  response,  and  the  velocity 
coupled  response.-; 

The  results  of  calculations  with  the  steady  state  burning  rate  model  are 
compared  with  the  experimental  results  of  Miller  (6,7)  for  twenty  one  HTPB/AP 
propellants  containing  a wide  range  of  oxidizer  particle  size  distributions. 

The  calculated  values  of  both  the  burning  rate  and  the  pressure  exponent  are  i i 
good  agreement  with  the  experimental  results.  The  poorest  agreement  with  the 
experimental  results  occurs  for  propellants  that  contain  a very  wide  overall 
particle  size  distribution.  It  is  demonstrated  that  the  poorer  agreement  is 
possible  due  to  the  lack  of  consideration  of  flame  interactions  between  the 
flames  of  individual  oxidizer  particles. 

The  steady  state  burning  rate  model  was  extended  to  include  erosive  burn- 
ing effects  by  coupling  a iiiodi'l  for  the  turbulent  boundary  layer  over  a flat 
plate  with  blowing  to  the  steady  state  burning  rate  model.  The  erosive  burn- 
ing model  is  based  on  the  hypothesis  that  the  erosive  effect  is  due  to  enhance 
mont  of  the  transport  properties  in  the  region  of  the  gas  phase  reaction  zone. 
The  enhancement  of  the  transport  properties  Is  assumed  to  be  due  to  the 
presence  of  a turbulent  boundary  layer. 

The  results  of  calculations  with  the  erosive  burning  model  are  not  com- 
pared with  experimental  data.  However,  calculations  are  presented  which 
illustrate  that  the  model  correctly  predicts  the  experimentally  observed 
trends  in  erosive  burning.  A critical  particle  diameter  for  erosive  burning 
IS  defined  since  the  combustion  of  oxidizer  particles  of  diameter  below  the 
critical  diameter  is  not  influenced  by  the  presence  of  a cross  flow  velocity. 

Two  pressure  c;'upled  response  models  are  presented.  One  is  based  on  the 
Cohen  hypothesis  whicli  related  the  resonant  frequency  and  peak  magnitude  of 
the  pressure  coupled  response  function  of  the  oxidizer  particle  diameter.  The 
other  model  is  based  on  the  Zel ' dovi ch-Novozhi lov  formalism  which  relates  the 
nonsteady  state  propellant  properties  to  steady  state  propellant  properties. 

The  theoretical  pressure  coupled  response  results  compare  favorably  with 
the  experimental  results  of  several  propellants.  The  addition  of  small  par- 
ticles to  the  oxidizer  particle  size  distribution  of  a propellant  is  shown  to 
increase  the  propellant's  resonant  frequency  and  the  magnitude  of  the  response 
Conversely,  the  addition  of  coarse  particles  is  shown  to  decrease  the  resonant 
frequency  and  the  peak  magnitude  of  the  response.  A sample  calculation  is  pre- 
sented which  illustrates  the  dramatic  effect  of  the  oxidizer  particle  size 
distribution  on  the  shape  of  the  curve  which  represents  the  pressure  coupled 
response  as  a function  of  frv.quency. 

The  model  for  the  velocity  coupled  response  is  based  on  the  previously 
described  erosive  burning  model  and  the  relation  between  the  pres.sure  coupled 
response  and  the  velocity  coupled  response  as  derived  by  Lengelle  (43). 

The  results  of  calculations  with  the  velocity  coupled  response  model  in- 
dicate th-»t  the  addition  of  coarse  particles  to  the  oxidizer  particle  size 
distribution  of  a propellant  increases  the  potential  for  velocity  coupling. 
Mcreuver,  it  is  shewn  that  the  velocity  coupled  response  vanishes  as  frequency 
is  increased. 
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NOMENCLATURE 


English  Symbols 


Empirical  constant  in  the  Denison  and  Baum  nonsteady 
analysis 

Arrhenius  frequency  factor  for  the  AP  flame 

Arrhenius  pre-exponential  for  fuel  pyrolysis 

Arrhenius  pre-exponential  for  oxidizer  pyrolysis 

Arrhenius  frequency  factor  for  the  primary  flame 

Empirical  constant  in  the  Denison  and  Baum  nonsteady 
analysis 

Blowing  parameter  defined  by  equation  257 

"Bunsen  burner"  diameter  in  diffusion  flame  analysis 

Empirical  constant  defined  in  equation  135 

Skin  friction  coefficient 

Skin  friction  coefficient  with  no  blowing 


Constant  in  the  equation  for  the  oxidizer  ignition  delay 

Specific  heat  of  the  propellant 

Constant  described  by  equation  281 

Prandtl  mixing  length  constant 

Oxidizer  particle  diameter 

Mean  intersection  diameter  of  an  oxidizer  particle  with 
the  fuel  plane 

Distance  of  the  fuel  regression 
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p,d,k 


^AP 

KpF 
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ox 


NOMENCLATURE  (cont'd.) 

Distance  of  the  oxidizer  regression 

Molecular  gas  diffusion  coefficient 

Turbulent  diffusion  coefficient  of  gas  below  the  final 
flame 

Turbulent  diffusion  coefficient  of  gas  below  the  primary 
flame 

Activation  energy  for  the  AP  flame 

Activation  energy  for  fuel  pyrolysis 

Activation  energy  for  oxidizer  pyrolysis 

Activation  energy  for  the  primary  flame 

Oxidizer  particle  size  distribution  function  given  by 
equation  153 

Oxidizer  particle  size  distribution  function  for  particles 
in  mode  j 

Distribution  function  defined  by  equation  123 

Resonant  frequency 

See  Figure  21 

AP  flame  rate  constant 

Primary  flame  rate  constant 

Thermal  conductivity  of  combustion  gases 

Flame  height  for  determining  the  average  eddy  diffusivity 

The  number  of  distribution  modes  of  oxidizer  specie  k 

Natural  log  of  the  o>;idizer  particle  weight  mean  diameter 

Mass  flux  from  an  oxidizer  particle  fuel  surface  pair 

Mass  flux  from  the  fuel  surface 

Mass  flux  from  the  oxidizer  surface 

Total  mass  flux  from  the  propellant  burning  surface 
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III 


d.k 


III 
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III 


o.k.j 


ni 

dm 


o.d.k 
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0 .d  ,k  • .1 


N 
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AN 


(i,d.k 


dN 


o,d 


dN 


p.d,k 


P 


P' 


NOMLNCLAIURL  (Coiifd) 

AvoiMip'  mass  f lux  per  unit  plaiUAr  biiniliuj  surface  area 

Mass  flux  at  any  point  on  tfie  luirninij  surface 

Mass  flux  per  unit  liurnint.i  surface  area  from  the  subarea 
contaiiiimi  oxidi/er  particles  of  si;e  between 

0 and  U I dll  and  oxidizer  species  k 

Average  of  the  above 

Average  nviss  flux  per  unit,  planar  surface  area  from  a 

siibsuriace  containiinj  oxidizer  crystals  of  size  between 

11  and  11  ♦ dll  and  oxidizer  species  k 

0 0 0 

lota)  mass  of  oxidizer  in  mode  j and  oxidizer  specie  k 

1 luctuatiiu)  mass  flux 

1 he  elemental  mass  of  oxidizt'r  wiDi  diameter  between 

P and  1)  ♦ dl)  and  oxidizer  specie  k 

0 0 0 

Ihe  eleiiuMital  mass  of  oxidizer  with  diameter  between  H 

and  P ♦ dP  in  distribution  mode  J and  oxidizer 
s fieri 8 k 

lota)  number  of  f'articles  at  the  biirniiui  surface 
Pressure  I'xponent 

Number  of  (lartieles  at  tlie  burniiii)  surface  of  size 
between  t\j  and  Po  i dll^  and  of  species  k,  per  unit  of 
planar  burniiu)  surface  area 

Ihe  number  of  oxidizer  crystals  per  unit  volume  with 
diameters  between  P^^  and  P^  * dO^ 

Number  of  particles  at  the  burninu  surface  of  size  between 
l\,  and  1\,  ♦ dP^j  and  species  k per  unit  of  planar  buniinq 
surface  area 

Pressure 

I luctuat  iiui  pressure 

Oxidizer  net  surface  heat  release 

Heat  to  vaporize  tlie  fuel  binder 

Heat  flux  at  the  surface  from  tlie  AP  and  final  flame 
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NOMENCLATURE  (Cont’d) 


Cri 

R 

Re 


P'.r 

V.1 

V.d.k 

^v,d,k 


* 

^o,d,k 


^P.d.k 


^S.d.k 


Heat  flux  at  the  surface  from  the  primary  flame 
See  Figure  21 

Reynolds  number  based  on  down  stream  position 
Pressure  coupled  response 
Velocity  coupled  response 

Real  part  of  the  pressure  coupled  response  function 

Imaginary  part  of  the  pressure  coupled  response  function 

Pressure  coupled  response  of  a pseudopropel lant 

Velocity  coupled  resonse  of  a pseudopropellant 

Average  propellant  burning  rate 

Burning  rate  of  a pseudopropellant 

Burning  surface  area 

Planar  burning  surface  area 

Fuel  surface  area 

Oxidizer  surface  area 

Planar  fuel  surface  area 

Planar  oxidizer  surface  area 

Planar  total  burning  surface  area 

The  exposed  surface  area  of  oxidizer  in  a psuedopropellant 
with  oxidizer  particles  of  size  between  Dq  and  + dll^  and 
oxidizer  specie  k 

Planar  surface  area  of  a psuedopropellant  containing 
('xidizer  particles  of  size  between  and  + dU^  and 
oxidizer  specie  k 

Number  of  oxidizer  types 

Burning  surface  area  of  a psuedopropellant  containing 
oxidizer  particles  of  size  between  0^  and  + dD^  and 
oxidizer  specie  k 
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NOMENCLATURE  (Cont'd) 


PlaniJr  surface  area  of  a psuedoi'ropel  lant  containing 

oyidizer  particles  of  size  betwiuMi  D and  D + dD„  and 
■ . . , . 0 0 0 
oxulwor  species  k 

Average  of  AS  . . 

P tU  » K 

* 

The  average  value 

T emperaturo 

Surface  tempercture 

Adiabatic  AP  flame  temperature 

Adiabatic  final  flame  temperature 

femperaturo  a^  the  location  of  the  AP  flame 

Burn  time 

Oxidizer  particle  ignition  delay 

Uiametcr  exponent  in  the  equation  for  the  oxidizer 
ignition  delay 

Velocity  parallel  to  the  burning  surface 

rime  average  velocity  parallel  to  the  burning  surface 

Rocket  chamber  core  gas  velocity 

fluctuating  velocity  parallel  to  the  burning  surface 

Velocity  normal  to  the  burning  surface 

lime  average  velocity  normal  to  the  burning  surface 

fluctuating  velocity  normal  to  the  burning  surface 

Pressure  exponent  in  the  equation  for  the  oxidizer 
iqnition  delay 

location  of  the  AP  flame  heat  release 
location  of  tlie  final  flame  heat  release 

AP  flame  standoff  distance 


Final  flame  standoff  distance 


NOMENCLATURE  (Cont'd.) 

Primary  flame  diffusion  distance 

Primary  flame  kinetics  distance 

Natural  log  of  the  oxidizer  particle  diameter 

Oxidizer  particle  number  density 


Greek  Symbols 
Oxidizer  mass  fraction 

Mass  fraction  of  oxidizer  specie  kin  the  psueciopropcllant 


Species  concentration  term  in  the  diffusion  flame  analysis 
Erosive  burning  strength 

Oistance  between  the  inner  and  outer  annulus  of  the 
"Bunsen  burner"  configuration 

Eddy  diffusivity  of  momentum 

Volume  fraction  of  oxidizer 


The  volume  fraction  of  a psuedopropel 1 ant  containing 
oxidizer  particles  of  size  between  Dq  atid  Do  + dOo  and 
specie  k 

The  volume  fraction  of  oxidizer  particles  with  diameter 
between  Dq  and  Oq  ♦ dDo 

The  fraction  of  particles  in  the  propellant  with  diameter 
between  Dq  and  Do  + dDo  oxidizer  type  k 

The  complex  frequency  parameter 

The  thermal  conductivity  of  the  propellant 

The  imaginary  part  of  A 

The  real  part  of  A 

Turbulent  thermal  conductivity  of  gas  below  the  AP  flame 


NOMENCLATURE  (Cont'd) 


Turbulent  thermal  conductivity  of  gas  below  the  final 
flame 

Turbulent  thermal  conductivity  of  gas  below  the  primary 
flame 

Stoichiometric  related  variable 
Viscosity  of  combustion  gases 
Density  of  the  core  gas 
Density  of  the  fuel  binder 
Density  of  tne  oxidizer 
Density  of  the  oxidizer  specie  k 
Density  of  the  propellant 
Propellant  density 

Density  of  the  pseudopropellant 

★ 

Natural  log  of  o 

Standard  deviation  of  oxidizer  particle  size  in  the  log 
normal  distribution 

Life  time  of  an  oxidizer  particle  at  the  burning  surface 

Shear  stress  at  the  burning  surface  due  to  a cross  flow 
velocity 

Dimensionless  frequency 

Frequency  of  the  gas  phase  pressure  oscillations 
Non-dimensional  flame  standoff  distance 

Non-dimensional  AP  flame  standoff  distance 

Non-dimensional  primary  flame  standoff  distance 

Non-dimensional  final  flame  standoff  distance 


INTRODUCTION 


Composite  solid  propellants  are  currently  used  in  the  majority  of 
the  solid  propellant  propulsion  systems.  Composite  propellants  are  a 
mixture  of  a plastic  fuel  binder  and  finely  ground  oxidizer  crystals. 

The  fuel  binder  is  initially  a liquid  but  is  polymerized  to  a solid 
after  mixing.  The  oxidizer  crystals  range  in  size  from  less  than  a 
micron  in  diameter  to  several  hundred  microns. 

The  mixing  properties,  combusC'cn  properties,  and  mechanical  pro- 
perties of  composite  propellants  can  be  changed  by  changing  the  mixture 
ratio  of  binder  to  oxidizer  and  by  changing  the  overall  solids  particle 
size  distribution  within  the  propellant  mixture. 

In  the  formulation  of  composite  propellants  to  satisfy  specific 
mission  requirements,  the  propellant  chemist  is  faced  with  the  formidable 
task  of  satisfying  several  stringent  simultaneous  constraints. 

The  basic  variables  the  propellant  chemist  has  at  hand  for  formu- 
lating composite  propellants  arc  the  amount  and  type  of  binder;  the 
amount,  type  and  particle  size  distribution  of  the  oxidizer;  the  amount 
and  type  of  burning  rate  catalyst;  and  the  amount,  type  and  particle 
size  distribution  of  metal  additives. 

In  determining  the  propellant  formulation,  the  propellant  chemist 
must  first  consider  the  desired  ballistic  properties.  That  is,  the 
propellant  must  burn  at  a specified  rate  at  any  given  pressure  and 
initial  propellant  temperature  within  some  specified  range  of  pressures 


and  initial  propellant  temperatures.  The  particle  size  of  the  oxidizer 
has  a profound  effect  upon  the  ballistic  characteristics. 

Currently,  much  attention  is  being  given  to  nozzleless  rockets. 

In  this  configuration  the  propellant  is  subjected  to  very  high  cross 
flow  velocites  which  cause  enhancement  of  the  propellant  burning  rate. 
This  effect  is  termed  erosive  burning.  It  is  suggested  that  the 
oxidizer  particle  size  distribution  has  a pronounced  effect  on  the 
erosive  burning  characteristics  of  composite  propellants.  Thus,  the 
propellant  formulator  may  have  to  meet  constraints  regarding  both  the 
steady  state  and  the  erosive  burning  characteristics  of  the  propellants. 

In  addition  to  the  importance  of  the  steady  state  ballistics,  the 
nonsteacly  burning  characteristics  of  the  propellant  are  also  significant, 
since  the  combustion  processes  may  interact  with  spurious  pressure 
oscillations  present  in  the  combustion  chamber.  That  interaction  can 
be  such  that  the  amplitude  of  the  pressure  oscillations  is  increased 
resulting  in  a condition  of  unstable  combustion  so  dest-uctive  that  the 
rocket  engine  itself  is  destroyed. 

Spurious  cross  flow  velocity  oscillations  can  result  in  a similar 

* 

condition  of  unstable  combustion.  Recent  work  (1-5)  suggests  that 
both  types  of  unstable  combustion  are  strongly  dependent  on  the  oxidizer 
particle  size  distribution.  Thus,  the  propellant  chemist  may  be  forced 
to  meet  constraints  on  combustion  instability. 

The  processability  of  the  propellant  must  also  be  considered.  Vari- 
ations of  the  particle  size  distribution  have  a pronounced  effect  on  the 


* 

Numbers  in  parenthese  indicate  references  listed  in  the  List  of  References. 


after  mix  viscosity  of  tlio  propellant.  The  propellant  must  have  a vis- 
cosity low  cnoiKjh  to  permit  it  to  flow  into  every  corner  of  the  rocket 
chamber  in  order  to  insure  repeatabi 1 i ty  of  performance  of  the  pro- 
pulsion system. 

Having  determined  a propellant  that  meets  the  ballistics  and  pro- 
cessing constraints,  the  propellant  must  now  also  meet  constraints  on 
its  mechanical  properties.  Those  properties  are  also  influenced  heavily 
by  the  particle  size  distribution  of  the  solid  oxidizer  as  wt.-ll  as  th<? 
metal  additive. 

Unfortunately,  definitive  a priori  methodologies  relating  ballistics 
processability  and  meclianical  properties  do  not  exist.  Thus,  the  pro- 
pellant forimilator  is  left  only  with  empiricisms  and  his  own  intuition 
to  guide  in  the  mixing  of  propellants  that  meet  specific  mission  .vpiiro- 
monts. 

It  is  evident,  then,  that  a combustion  model  which  describes  the 
steady  and  nonsteady  burning  processes  of  compost ce  solid  propellants, 
with  particular  emphasis  on  particle  size  effects,  should  result  in  a 
valuable  tool  for  the  propellant  formulator.  The  combustion  model  when 
used  in  conjunction  v/ith  available  empirical  data  (6,7)  on  the  effect 
of  oxidizer  particle  size  and  size  distribution  on  the  propellant's 
viscosity  and  mechanical  properties  should  result  in  a complete  descrip- 
tion of  the  desired  propellant  formulation. 

The  objective  of  this  research  program  is  to  develop  such  a 
combustion  model  for  composite  propellants. 
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KlVItW  OF  THt  LllLRATURE 

GofU’ial  Hi sciiss ion 

Tht’  fcview  of  ttu*  literatui'o  is  concorneO  with  tlioso  poftiiient 
references  pertainimj  to  the  moOels  of  the  steady  state  combustion  of 
solid  composite  propellants  as  well  as  erosive  burniini  and  nonsteady 
state  combustion,  before  reviewiny  ttie  models  a brief  description  will 
bo  given  of  the  composite  solid  profiellant  anil  the  physii.al  and  che.iiical 
processes  that  occur  neai'  tiu'  burning  surface  of  that  composite  propel- 
lant. 

The  Composite  Solid  Fi'opel  laiU 

The  composite  propellant  is  a heterogeneous  mixture  of  finely 
ground  oxidizer  particles,  with  a plastic  fuel  binder  which  is  inittally 
a liquid  but  is  finally  polymerized  to  a solid.  The  o.xidizer  particles 
used  in  composite  propellants  are  usually  yroutui  ammonium  perchlorate 
(AP)  crystals  with  a particle  size  range  of  from  less  than  a micron  in 
diameter  to  diameters  as  large  as  several  hundred  microns.  Other  oxi- 
dizers which  are  used  to  a lessor  extent  include  cyclotetramethy lene- 
(etrani tramine  (HIIX),  hydroxytammoiiium  perchlorate  (liAP),  potassium 
perchlorate  (kP)  and  aninonium  nitrate  (AM). 

The  number  density  of  oxidizer  particles  in  a composite  propellant 
usually  varies  in  a log  normal  distribution  fashion  with  the  size  of 
the  particle.  This  variation  of  particle  number  density  with  particle 
diameter  is  termed  a particle  size  distrihiition.  Figure  1 depicts  a 

IP 


typical  particle  size  distribution  for  a nominal  diameter  of  ten  microns. 
Composite  propellants  usually  contain  several  distribution  modes. 

The  fuel  binder  used  in  the  formulation  of  composite  propellants 
is  typically  hydroxyl -terminated  polybutadiene  acrylic  acid  (HTPB), 
other  binders  include,  dei'ivatives  of  the  polybutadienes,  polysulfides, 
polyurethanes,  and  polyvinyl  chlorides. 

In  many  inst.ances  additives  such  as  burning  rate  catalysts  and  metal 
powders  arc  auded  to  the  propellant  to  modify  the  ballistic  properties 
of  the  propellant. 

Description  of  the  Combustion  of  CojiipqyJe  Solid  _P rope IJants 

The  combustion  of  composite  solid  propellants  is  a very  complex 
process.  This  is  primarily  due  to  the  heterogeneous  nature  of  the  com- 
posite. Consequently , many  aspects  of  the  combustion  process  are 
determined  by,  or  at  least  dependent  upon,  the  propellant’s  physical 
structure.  Therefore,  in  an  attempt  to  describe  the  combustion  of 
composite  propellants  it  is  advantageous  to  first  look  at  the  propellant's 
physical  characteristics. 

The  sketch  in  Figure  2 depicts  the  structure  of  a typical  composite 
propellant  mixture.  Composite  solid  propellants  are  literally  a packed 
bed  of  polydisperse  oxidizer  particles  where  the  interstitial  voids 
between  the  packing  of  the  particles  are  filled  with  a continuous  binder. 

The  oxidizer  particles  are  ground  to  a desired  size  before  the  pro- 
pellant is  mixed.  However,  the  oxidizer  particles  formed  during  the 
grinding  process  are  not  of  one  size  but  rather  are  a distribution  of 
sizes.  This  distribution  can  be  characterized  by  a log  normal  distri- 
bution. 
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Figure  3 shows  the  effect  of  the  distribution  width  parameter  on 
a typical  particle  size  distribution.  In  Figure  3 particle  number 
density  is  plotted  as  a function  of  diameter.  The  weight  median 
(or  50S)  diameter  is  given  by  that  diameter  where  the  number  density 
curve  reaches  a maximum.  The  width  of  the  distribution  is  representa- 
tive of  the  number  of  particles  differing  in  size  from  that  of  the 
weight  median  diameter. 

A distribution  where  all  of  the  particles  are  of  the  same  size  is 
termed  monodisperse.  Such  a distribution  would  plot  as  a vertical 
straight  line  located  at  the  we'ght  median  diameter.  That  is,  no 
particles  would  differ  in  size  from  that  of  the  median  diameter. 

A polydisperse  distribution,  on  the  other  hand,  contains  particle 
sizes  other  than  that  of  the  median  diameter.  Furthermore,  the  wider 
the  distribution,  the  more  particles  there  are  larger  than  and  smaller 
than  the  median  diameter.  A wider  distribution  is  depicted  in  Figure  3 
by  the  dashed  line.  A more  complete  explanation  of  the  particle  size 
distribution  can  be  found  in  Reference  (8). 

Composite  propellants  contain  a high  loading  or  weight  percent  of 
oxidizer  particles  so  that  high  values  of  specific  impulse  can  be 
achieved.  In  order  to  contain  a high  oxidizer  loading,  the  oxidizer 
must  be  made  up  of  several  different  oxidizer  grinds  in  which  the 
smaller  particles  can  fill  the  voids  between  the  packing  of  the  larger 
particles.  Thus,  it  should  be  clear  that  propellants  contain  a poly- 
dispersion of  oxidizer  particle  sizes.  For  those  propellants  it  has 
been  demonstrated  (9,10)  that  above  some  length  scale  characteristic  of 
the  particle  size  the  arrangement  of  the  particles  within  the  propellant 
is  totally  random  due  to  complete  mixing  during  the  propellant  formulation. 
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The  combustion  process  can  best  be  described  by  following  the  evo- 
lution of  one  o;xidizor  particle  from  deep  within  the  propellant  until 
it  emerges  at  the  surface  and  is  consumed  in  the  combustion  process. 

Since  aninonium  perchlorate  is  the  most  extensively  used  oxidizer 
in  composite  propellants,  attention  is  directed  toward  that  oxidizer 
throughout  the  following  discussion. 

A particle  in  the  undisturbed  propellant  has  the  initial,  undis- 
turbed propellant  temperature.  As  the  burning  surface  moves  toward  the 
particle,  the  particle  reaches  the  thermal  wave  penetration  depth  where 
upon  it  begins  to  receive  energy  from  the  burning  surface.  Eventually 
the  oxidizer  particle  reaches  a temperature  of  513°K  undergoing  an  endo- 
thermic change  of  crystal ine  state  from  orthorhombic  to  cubic. 

When  the  particle  emerges  at  the  surface,  it  begins  an  exothermic 
decomposition  to  gaseous  products. 

The  combustion  process  of  the  propellant  in  the  immediate  vicinity 
of  an  exposed  oxidizer  particle  involves  several  different  chemical 
reactions.  Each  of  the  important  steps  in  the  overall  chemical  reaction 
will  be  isolated  for  the  purpose  of  discussion: 

1.  The  exothermic  heterogeneous  reaction  between  the  binder  and 
the  oxidizer  interface. 

2.  The  endothermic  pyrolysis  of  the  binder  or  fuel. 

3.  The  exothermic  decomposition  of  the  crystal. 

4.  The  monopropellant  reaction  between  the  decomposition  products 
of  the  oxidizer. 

5.  The  exothermic  gaseous  reaction  between  the  oxidizer  decomposi- 
tion products  and  the  fuel  pyrolysis  products. 


As  the  oxidizer  particle  emerges  at  the  burning  surface,  some  of 
the  oxidizer  undergoes  a heterogeneous  subsurface  reaction  with  the 
surrounding  fuel  binder  (11).  There  remains,  however,  a great  deal  of 
controversy  in  the  literature  pertaining  not  only  to  the  importance  of 
a subsurface  reaction  but  also  to  the  actual  existence  of  such  a 
reaction. 

The  strongest  case  for  heterogeneous  subsurface  reactions  comes 
from  ignition  studies  of  Anderson  (12,13)  and  co-workers.  In  those 
studies,  oxidizers  such  as  fluorine  and  chlorine  trifluoride  were  used 
to  ignite  composite  propellants  hypergol ical ly.  The  experiments  demon- 
strated that  the  ignition  process  could  proceed  as  a result  of 
spontaneous  heterogeneous  reactions  occurring  at  the  oxidizer-binder 
interface.  However,  extensive  electron  microscope  studies  at  Lockheed 
(14)  on  slices  of  extinguished  propellant  samples  have  shown  no  trace  of 
subsurface  heterogeneous  reaction.  Moreover  calculations  by  Beckstead 
(14)  and  co-workers  have  indicated  that  heterogeneous  subsurface 
reactions  are  not  required  to  sustain  combustion. 

A portion  of  the  fuel  or  binder  surrounding  the  oxidizer  crystal 
can  be  associated  with  that  crystal.  As  the  crystal  nears  the  burning 
surface,  the  surrounding  binder  located  at  the  surface  pyrolyzes  due  to 
the  heat  transfer  to  the  binder  from  the  various  exothermic  reactions 
occurring  nearby.  The  result  is  the  thermal  degradation  of  small  frag- 
ments of  the  polymer  binder  which  serve  as  reactants  in  the  gaseous 
fuel-oxidizer  reaction  zone  above  the  burning  surface. 

At  the  burning  surface  the  oxidizer  initially  undergoes  a decom- 
position step  from  solid  crystal  to  gaseous  ammonia  and  perchloric  acid. 
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The  perchloric  acid  then  further  decomposes  (15)  to  oxygen,  chlorine 
monoxide,  and  hydroxyl.  The  reactions  are  as  follows: 

NH^CIO^  - NHj  + HCIO^  1 

HCIO4  - OH  + CIO  + O2  2 

The  AP  monopropellant  flame  occurs  follow'ng  the  decomposition  of  the 
perchloric  acid.  The  AP  monopropellant  flame  is  a kinetics  controlled 
reaction  between  the  ammonia  (which  resulted  from  the  first  step  in  the 
AP  decomposition)  and  the  chlorine  monoxide  (which  resulted  from  the 
perchloric  acid  decomposition).  The  reaction  produces  oxygen  and  inert 
products. 

The  final  reaction  to  be  considered  is  the  gas  phase  reaction  be- 
tween  the  pyrolized  fuel  binder  and  the  decomposed  oxidizer.  The 

reaction  may  result  from  either  or  both  of  two  paths.  The  first  is  the  ; 

reaction  between  the  binder  pyrolysis  products  and  the  perchloric  acid 

I 

from  the  decomposition  of  the  oxidizer.  This  results  when  the  perchloric  j 

acid  and  binder  pyrolysis  products  diffuse  and  mix  near  the  burning 

I 

surface.  The  second  path  is  the  reaction  between  the  binder  pyrolysis  j 

products  and  the  products  of  the  AP  monopropellant  flame.  This  flame  ^ 

I 

results  when  the  binder  pyrolysis  products  and  the  products  of  the  AP  » 

monopropellant  flame  diffuse  and  mix.  The  heat  generation  from  these  j 

two  reactions  occur  at  varying  distances  above  the  burning  surface  due  i 

to  the  irregular  nature  of  the  propellant  burning  surface  and  the  1 

inhomogeneities  of  the  reaction  zone.  Both  flames  are  strongly  depen- 
dent upon  the  diffusional  mixing  of  the  reactants  and,  trierefore,  are 
very  dependent  upon  the  oxidizer  particle  size. 
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Quite  clearly,  the  combustion  process  of  composite  solid  propellants 
is  a very  complex  situation.  Attempts  to  model  this  process  have  evolved 


from  very  crude  and  simple  models  of  the  late  50' s which  only  attempt 
to  predict  gross  trends,  to  very  sophisticated  models  of  the  70's  which 
attempt  to  predict  quanti tiati ve  results.  Of  the  several  models  which 
have  attempted  to  describe  the  combustion  of  composite  propellants, 
three  models  have  received  the  most  attention  as  well  as  acceptance  in 
the  solid  propellant  propulsion  community.  Those  three  models  are,  in 
chronological  order:  the  Summerfield  Granular  Diffusion  Flame  Model 
(16),  the  Heterogeneous  Reaction  model  of  Hermance  (11,17)  and  the 
Multiple  Flame  Model  of  Beckstead,  Derr  and  Price  (18).  A complete 
description  of  each  of  the  above  models  follows. 

Composite  Propellant  Combustion  Models:  Steady  ^ta^ 

A comprehensive  model  which  describes  the  combustion  of  a solid 
propellant  should  be  based  on  the  key  physical  and  chemical  processes 
involved  in  the  combustion  zone.  Those  processes  should  include  the 
heat  release  and  associated  kinetics  due  to  steps  such  as  fuel  pyrolysis, 
oxidizer  decomposition  and  gas  phase  flame  reactions  including  the  dif- 
fusional  mixing  process.  It  should  also  account  for  the  effect  of 
oxidizer  partical  size  and  concentration  on  propellant  burning  rate. 

The  Granular  Diffusion  Flame  Model 
The  granular  diffusion  flame  (GDF)  model  was  formulated  by 
S'lntnerf ield  (16)  to  describe  the  burning  characteristics  of  anmonium 
perchlorate  (AP)  based  composite  propellants  and  is  of  historical  value 
since  it  was  the  first  attempt  to  model  the  combustion  of  composite 
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propellants.  The  model  was  primarily  for.nulated  in  an  attempt  to  pre- 
dict the  effects  of  pressure,  oxidizer  particle  size,  and  the  ratio  of 
fuel  to  oxidizer  on  the  combustion  process. 

Description  of  the  GDF  Model.  The  physical  structure  upon  which 
the  GDF  model  is  based  is  presented  in  Figure  4a.  The  model  is  a one 
dimensional  model  which  assumes  that  the  decomposition  process  at  the 
propellant  surface  is  controlled  by  conductive  heat  feed  back  from  an 
oxidizer/fuel  flame  occurring  in  the  gas  phase.  In  the  model's  original 
form  in  1960  (16),  a premixed  ammonia/perchloric  acid  reaction  zone  (AP 
monopropellant  flame)  .vas  assumed  to  take  place  at  the  surface  of  the 
propellant  followed  by  a reaction  between  the  oxidizing  products  from 
the  ammonia/perchloric  acid  reaction  and  the  fuel  products  from  the 
pyrolized  binder. 

The  critical  assumption  of  the  model  was  that  vapors  of  fuel  or 
oxidizer  or  both  are  released  in  the  form  of  pockets  of  a certain  mass 
content.  The  average  mass  content  of  a pocket  was  assumed  to  be  much 
smaller  than  that  of  an  average  oxidizer  crystal  and  to  be  independent 
of  the  temperature  and  pressure.  However,  the  size  of  the  pockets  was 
assumed  to  be  related  to  the  average  oxidizer  particle  size.  The  rate 
at  which  the  pockets  are  consumed  in  a chemical  reaction  as  they  pass 
through  the  flame  zone  was  assumed  to  be  controlled  by  diffusional 
mixing  and  the  chemical  kinetics.  The  burning  surface  was  assumed  to  be 
dry  such  that  the  mixing  of  the  oxidant  and  fuel  occurs  only  in  the 
gas  phase. 

Development  of  the  Burning  Rate  Equation.  The  propellant  burning 
rate  relationship  for  the  GDF  model  is  based  on  an  energy  balance  at 
the  burning  propellant  surface.  The  temperature  distribution  through 
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Figure  4.  The  GDF  Physical  Description 
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the  propellant  and  the  gas  phase  reaction  zone  is  presented  in  Figure 
4b.  The  significant  feature  of  Figure  4b  is  that  a linear  temperature 
profile  is  assumed  in  the  gas  phase  reaction  zone.  Thus,  the  temper- 
ature gradient  above  the  burning  surface  is  given  by  (T^  - T^l/L  where 
Tf  is  the  flame  temperature,  is  the  average  temperature  of  the 
burning  surface  und  the  thickness  of  the  gas  phase  reaction  zone  is 
given  by  L. 

- An  energy  balance  at  the  burning  surface,  as  shown  in  Figure  5, 
can  be  written.  Thus 


»,(Tf  - T^)/L  - 


where  X is  the  thermal  conductivity  of  the  hot  combustion  gases,  Cp  is 
the  average  specific  heat  of  the  propellant,  T^  is  the  initial  temper- 
ature of  the  propellant,  m^.  is  the  mass  flux  of  propellant  consumed  per 
unit  of  surface  area  and  is  the  surface  heat  release  per  unit  area 

and  per  pound  of  propellant  consumed  at  the  burning  surface.  The  sur- 
face heat  release  includes  the  heat  required  to  vaporize  the  fuel  binder, 
the  net  heat  released  by  the  AP  decomposition  and  the  net  heat  released 
by  the  heterogeneous  surface  and  subsurface  reactions. 

The  only  remaining  problem  is  to  determine  the  gas  phase  reaction 
zone  thickness,  L.  In  the  determination  of  L,  two  pressure  extremes 
were  considered,  low  and  high. 

Since  the  mass  diffusion  coefficient  for  the  reaction  species  is 
inversely  proportional  to  pressure  while  the  chemical  reaction  rate  is 
directly  proportional  to  pressure,  the  molecular  diffusion  rate  is  much 
faster  than  the  chemical  reaction  rate  at  low  pressure.  Thus,  at  low 
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Figure  5.  The  Sumtierf ield  Surface  Energy  Balance  . | 

I 

I 


( 


! 


uressure,  the  gas  phase  reaction  occurs  in  a premixed  gas  and  the  re- 
action rate  is  the  controlling  mechanism.  Conversely,  at  high  pressure 
the  chemical  reaction  rate  is  much  faster  than  the  intermolecular  dif- 
fusion rate  so  that  the  combustion  mechanism  is  controlled  by  the  mixing 
process. 

In  the  determination  of  L for  the  low  pressure  premixed  flame, 
Sumiiierfield  treated  the  flame  as  a stream  of  velocity  u in  which  a 
second  order  reaction  was  taking  place.  The  mean  gas  velocity  in  the 
flame  zone  is  given  by 


u - mj/Pg  ^ 


where  is  the  density  of  the  gas  in  the  flame  zone. 

The  flame  zone  thickness  for  the  premixed  flame  is  given  by 


L = UT 


whore  t is  the  time  it  takes  the  chemical  reaction  to  occur.  Assuming 
a global  second  order  reaction  for  the  premixed  gases,  the  reaction 
time  is  given  by 

T = [(1  - f)^PgA  exp  (-E/RTg)]  6 

wnere  an  Arrhenius  expression  is  assumed  for  the  rate  of  reaction,  T^ 
is  the  average  gas  phase  reaction  zone  temperature  and  c is  a function 
of  the  products  of  reaction,  i.e.,  c is  equal  to  zero  at  the  propellant 
surface  and  c is  equal  to  unity  at  the  completion  of  the  reaction. 

For  simplicity,  Sunflierfield  assumed  that  the  average  value  for 
(1  - e)  over  the  flame  zone  is  unity.  Physically  this  means  that  most 
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of  the  “overall  reaction"  takes  place  near  the  propellant  surface. 

Combining  equations  4,  5 and  6,  yields  the  followiny  for  the  flame 
zone  thickness  for  the  low  pressure  premixed  flame: 


Lj  = nt^/[A  c\p  (-E/RTg)] 


7 


When  equation  3 is  solved  for  the  total  mass  flux  and  tno  result  sub- 
stituted into  equation  7,  the  following  expression  for  the  thickness  of 
the  low  pressure  flame  zone  is  obtained: 


X (Tj:  - T^) 

p“Ii‘  (T  - T^l  - Q lA'  exp  (-E/rTJ 
g*’  p s 0 s g 
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For  tiij  high  pressure  diffusion  flame,  Summerfield  assumed  that  the 
mass  of  a pocket  could  be  expressed 


where  d is  the  oxidizer  particle  diameter.  The  lifetime  of  a pocket  is 
determined  by  the  rate  of  gas  diffusion  to  the  surrounding  flame.  The 
molecular  di ffusivi t ies  of  the  fuel  and  the  oxidizer  are  assumed  to  be 
the  same  and  are  averaged  over  the  gaseous  reaction  zone.  Under  these 
conditions  the  lifetime  of  a fuel  pocket  is 


where  is  the  average  diffusivity  of  the  oxidizer  and  fuel  vapors. 
The  thicxness  of  the  flame  zone  is  the  product  of  the  velocity  of  the 
gas  stream  and  the  lifetime  of  llie  fuel  pocket.  Thus 


33 


Combining  equations  9and]1  yields  the  diffusion  flame  zone  thickness 

zn- 

u m 

t-g  g 

Solving  equation  3 for  the  total  mass  flux  and  substituting  this  into 
equation  12  yields 


(Tf  - 


Expressions  for  the  flame  zone  thickness  in  the  two  extreme  cases 
can  now  be  written.  For  the  general  intermediate  case,  Summerfield 
assumed  that  the  flame  zone  thickness  varies  with  pressure  partly  as  if 
it  were  reaction  rate  controlled  (L^)  and  partly  as  if  it  were  diffusion 
controlled  (L2).  and  that  the  flame  zone  thickness  could  be  expressed 
as  the  sum  of  the  above  two.  Thus,  the  general  expression  for  the  flame 
zone  thickness  becomes 


[c  (T  - T ) - 

**  p'  S 0* 


I X ^ 1/1 

/ r/nr  xW?  ;7l/2  5/6 
Pg[A  exp  (-E/RFg)]  Dg  Pg 


Substituting  equation  14  for  L into  the  energy  balance,  equation  3,  the 
following  expression  for  the  burning  rate  is  obtained: 
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1/2 


RT 


P[A  exp  (-E/RTg)] 


W 
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or  by  lumping  parameters,  equation  15  can  be  simplified  to 


16 


where  a is  related  to  the  chemical  reaction  time  and  b is  related  to  the 
diffusion  time.  Furthermore,  all  of  the  parameters  which  make  up  the 
two  constants  are  assumed  to  be  independent  of  pressure. 

Discussion  of  the  GDF  Model.  Sunmerfield  used  equation  16  to 
correlate  experimental  burning  rate  data  for  a range  of  styrene  base 
polyester-ammonium  perchlorate  propellant  formulations.  The  equations 
provided  excellent  agreement  between  theory  and  experiment  in  predicting 
the  effect  of  oxidizer  crystal  size,  fuel  to  oxidant  ratio,  and  pressure 
on  burning  rate. 

Although  the  burning  rate  equation  is  empirical  in  nature  and  does 
little  more  than  give  a means  for  curve  fitting  data,  it  was  the  first 
burning  rate  equation  based  on  a theoretical  premise.  Another  signifi- 
cant advantage  of  the  Summerfield  burning  rate  expression  was  that  it 
provided,  in  a compact  equation,  a means  of  estimating  the  effects  of 
varying  the  important  propellant  combustion  parameters  on  burning  rate. 
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The  model  was  later  modified  by  Summerfield  et  al.  (19)  to  take  into 
account  a distended  AP  monopropellant  flame.  This  was  an  attempt  to  make 
the  model  applicable  to  low  pressure  combustion. 

Miller  and  co-workers  (20)  extended  the  GDP  model  to  polydisperse 
propellants  by  applying  a statistical  technique. 

The  Hermance  Heterogeneous  Reaction  Model 

The  first  attempt  to  develop  a model  capable  of  predicting  in- 
fluences on  burning  rate  due  to  propellant  formulation  changes  was  that  of 
Hermance  (11). 

The  model  is  based  on  a detailed  combination  of  the  steady  state 
decomposition  processes  of  the  fuel  and  oxidizer.  The  model  was  first 
formulated  in  1966  and  later  revised  in  1967  (17). 

Hermance  postulates  an  exothermic,  heterogeneous  reaction  between 
the  solid  fuel  and  the  oxidizer  decompositon  products  and  assumes  that 
the  gas  phase  flame  position  is  the  sum  of  the  lengths  associated  with 
diffusional  mixing  and  chemical  reaction.  Incorporation  of  these  pro- 
cesses into  one-dimensional  energy  equation  allows  the  calculation  of  the 
burning  rate  as  a function  of  the  propellant  composition  and  ambient 
pressure.  Hermance’s  heterogeneous  reaction  (HR)  model  was  a significant 
advance  in  the  state  of  the  art  of  combustion  modeling  since  it  was  the 
first  model  to  deal  directly  with  the  heterogeniety  of  the  solid  phase 
by  taking  into  consideration  an  arbitrary  distribution  of  particle  sizes. 

Decription  of  the  HR  Model.  The  model  assumes  that  the  steady 
state  combustion  of  a composite  propellant  occurs  in  three  regions;  the 
propellant  surface,  the  gas  phase  flame  zone,  and  the  region  between 
the  surface  and  the  flame.  The  surface  chemical  processes  are 
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endothermic  fuel  pyrolysis,  exothermic  oxidizer  decomposition  and  an 
exothermic  heterogeneous  chemical  reaction  between  the  fuel  binder  and 
decomposed  oxidizer  in  small  regions  surrounding  individual  oxidizer 
crystals.  Each  of  these  reactions  produces  a mass  flux  from  the  pro-  > 
pellant  surface. 

In  the  gas  phase  'lame  zone,  an  exothermic  chemical  reaction  occurs 
between  the  oxygen  rich  products  from  the  decomposed  oxidizer  and  the 
pyrolized  fuel.  Heat  from  the  flame  zone  is  fed  back  to  the  propellant 
surface  by  conduction. 

In  the  region  between  the  surface  and  the  flame  zone,  the  surface 
decomposition  products  mix  by  diffusion  and  undergo  a chemical  kinetic 
delay  before  ignition  at  the  flame.  In  the  original  model  (11)  only  the 
ignition  delay  was  considered. 

The  processes  depend  on  the  pressure  level  and/or  the  temperature 
at  the  location  of  the  process  and  are  linked  together  by  the  temperature 
distribution  in  the  gas  and  solid  phases. 

The  physical  model  is  depicted  in  Figure  6a.  One  of  the  more 
important  aspects  of  the  model  is  the  calculation  of  the  intersection 
diameter  of  an  oxidizer  crystal  with  the  burning  surface  as  depicted 
in  Figure  6a.  If  the  propellant  burning  surface  is  visualized  as  a 
plane  passing  through  a randomly  packed  bed  of  spherical  oxidizer 
crystals,  it  can  be  shown  (21)  that  the  statistical  average  intersection 
diameter  is 

= /?7J  17 
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Figure  6.  Physical  Description  of  the  Hermance  Heterogeneous 
Reaction  Model 
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Development  of  the  Durning  Rate  Equation.  The  propellant  burning 
rate  can  be  determined  from  the  total  mass  flux  issuing  from  the  pro- 
pellant surface.  It  is  the  sum  of  the  mass  fluxes  of  gaseous  species 
produced  by  each  of  the  surface  decomposition  processes  multiplied  by 
the  fraction  of  the  total  propellant  surface  area  associated  with  each 
of  the  processes.  Thus 

^ ■ 1^  * "srlV«o>’ 

where  pp  is  the  density  of  the  propellant,  m is  the  mass  flux  and  S is 
surface  area.  The  subscripts  f,  ox,  o and  sr  designate  fuel,  oxidizer, 
total  propellant  and  surface  reaction  respectively.  Assuming  a planar 
surface,  the  ratio  of  fuel  and  oxidizer  surface  area  to  the  total  pro- 
pellant surface  area  can  be  expressed  as  a function  of  the  volume 
fraction  of  oxidizer  in  the  propellant.  Thus 


(1  - C) 


19 


and 


20 


The  area  on  which  the  surface  reaction  occurs  is  calculated  by 
postulating  that  an  oxidizer  crystal  decreases  in  size  during  decom- 
position as  shown  in  Figure  6a.  This  size  reduction  produces  a fissure 
of  depth  c between  the  oxidizer  crystal  and  the  fuel  binder  which  is  the 
region  where  the  heterogeneous  reaction  takes  place.  Arrhenius  type 
functions  are  used  to  describe  the  surface  regression  of  both  fuel  and 
oxidizer. 
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The  inodel  divides  the  steady  state  combustion  process  of  composite 
propellants  into  three  regions  as  shown  in  Figure  6b.  The  differential 
form  of  the  energy  equation  is  solved  for  each  region  simultaneously. 

The  required  boundary  conditions  result  by  matching  temperature  and 
heat  flux  at  the  interface  of  each  section.  The  analysis  results  in 
three  equations  with  three  unknowns.  They  are  the  burning  rate  equation, 
the  surface  temperature  equation  and  the  'lame  temperature  equation. 

These  equations  in  condensed  form  ar> 


r = a,  exp  (-E^/EO^)  + ^ Ct  - exp  [ 

Oj  P s 


] 21 


°s  ~ 


* 4^  - -i]  o»p  [— 
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EO. 


^2  '^f  * 

- — exp  (-^)  + a^  exp  (-E  ) 
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0f  = 0$  + a^Cl  - exp(-E  )]  23 

where  a^ , a^.  a^.  b^ , b^,  b^,  d^  are  constants  containing  the  physical 
and  chemical  properties  of  the  propellants.  The  burning  rate,  surface 
temperature  and  flame  temperature  are  calculated  by  numerical  iteration 
using  the  above  thrc'  equations. 

Discussion  of  the  HR  Model.  Hermance  solved  the  above  three  simul- 
taneous equations  with  a set  of  input  data  which  represented  the 
combustion  properties  of  a polysulfide  - AP  composite  propellant.  He 
obtained  a •'airly  reasonable  fit  of  the  theoretical  model  to  experimental 
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burning  rate  data  taken  over  the  pressure  range  from  1 to  400  atmos- 
pheres. The  model  qualitatively  predicts  the  effect  of  oxidizer 
particle  size  on  the  propellant  burning  rate  at  pressures  above  200 
atmospheres  but  the  model  over  predicts  the  effect  of  initial  propellant 
temperature  on  burning  rate. 

At  low  pressure  the  model  predicts  very  little  heterogeneous  reac- 
tion while  at  high  pressure  the  crevice  becomes  huge  and  the 
heterogeneous  reaction  dominates,  as  depicted  in  Figure  7.  Moreover, 
the  calculated  surface  temperature  is  almost  constant  with  increasing 
pressure  indicating  that  the  temperature  dependent  Arrhenius  type 
function  will  be  essentially  constant.  Thus,  the  burning  rate  char- 
acteristics predicted  by  the  HR  model  are  determined  almost  completely 
by  the  crevice  and  its  formation.  Yet  there  is  no  physical  evidence 
to  support  the  formation  of  a crevice  as  postulated  by  Hermance  (14). 

The  Beckstead  Derr  Price  Multiple  Flame  Model 

After  extensive  experimental  observation  of  the  surface  structure 
of  a burning  composite  solid  propellant,  Beckstead,  Derr  and  Price  for- 
mulated a multiple  flame  composite  propellant  combustion  model  based  on 
the  postulated  existence  of  three  flames  in  the  region  of  each  exposed 
oxidizer  particle. 

Description  of  the  Multiple  Flame  Model.  The  model  is  based  on 
the  assumption  that  the  gas-phase  heat  release  can  be  represented  by 
multiple  flames  surrounding  individual  oxidizer  crystals.  The  geometric 
relationship  between  the  oxidizer  crystal  and  the  binder  matrix  is 
evaluated  statistically  as  in  the  Hermance  HR  model  (11).  The  geometric 
relationship  of  the  multiple  flame  is  shown  in  Figure  8.  The  distance 
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Figure  7.  The  Pressure  Dependence  of  the  Hermance  Crevice 
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Figure  8.  The  Flame  Structure  for  an  AP  Composite  Propellant 


from  the  flame  to  the  surface,  the  standoff  distance,  is  an  important 
consideration  in  the  model. 

Three  flames  are  considered:  the  primary  flame,  the  AP  monopropel- 
lant flame  and  the  final  diffusion  flame.  The  AP  flame  occurs  between 
the  decomposition  products  of  the  AP  forming  O2  and  inert  products. 

Thus 


NH^  + HCIO^  Inert  Products  + O2  24 

Since  the  AP  decomposition  products  are  premixed,  the  AP  flame  standoff 
distance  is  dependent  only  upon  chemical  kinetics. 

The  primary  flame  occurs  between  fuel  pyrolysis  products  and  the 
AP  decomposition  products.  Thus 

25 

Fuel  Products  + HCIO*  (decomposition  products)  ->  Primary  Flame 

Combustion  Products 

The  primary  flame  height  is  dependent  on  diffusional  mixing  as  well  as 
kinetics. 

The  final  diffusion  flame  occurs  between  the  binder  decomposition 
products  and  the  oxidizing  products  from  the  AP  flame.  Since  the  re- 
actants of  the  final  diffusion  flame  are  preheated,  the  kinetics  are 
very  fast,  thus,  the  final  flame  height  is  dependent  only  on  diffusional 
mixing. 

The  flame  standoff  distances  are  a function  of  the  combustion 
pressure.  At  low  pressures  (<100  atm)  the  propellant  is  considered  to 
burn  as  a premixed  flame  with  the  oxidizer  and  binder  decomposition 
products  mixing  completely  before  a reaction  occurs.  As  the  pressure 
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increases,  the  nixing  path  length  increases,  and  the  reaction  path 
length  decreases  so  that  the  two  reaction  paths  (ammonia  + oxidizing 
products  and  binder  products  + oxidizing  products)  become  competitive. 
At  higher  pressures  {>100  atm)  the  ammonia  reacts  with  the  oxidizing 
products  before  the  binder  products  can  diffuse  into  the  oxidizer 
stream  and  react.  Due  to  the  fuel  rich  nature  of  composite  solid 
propellants,  a final  diffusion  flame  always  occurs  above  the  AP  flame. 
The  flame  structure  variation  with  pressure  is  depicted  in  Figure  9. 

At  the  propellant  surface  the  initial  decomposition  step  of  both 
the  binder  and  the  oxidizer  is  endothermic.  However,  while  still  ad- 
sorbed on  the  surface,  either  set  of  products  can  undergo  a condensed 
phase  reaction  before  passing  into  the  gas  phase.  Products  from  the 
surface  decompositions  then  pass  into  the  gas  phase  and  begin  mixing 
and  reacting. 

The  burning  surface  of  the  exposed  oxidizer  particle  is  assumed  to 
be  spherical  in  shape.  Due  to  an  ignition  delay,  the  burning  oxidizer 
surface  protrudes  above  the  planar  fuel  surface  at  low  pressure  but  is 
recessed  below  the  planar  fuel  surface  at  high  pressure.  The  oxidizer 
regression  is  assumed  to  be  the  overall  rate  controlling  mechanism,  and 
an  overall  average  temperature  is  defined  for  the  entire  burning  sur- 
face. In  the  model's  original  form,  the  oxidizer  is  assumed  to  be 
monomodal , monodisperse  and  spherical.  However,  the  model  has  since 
been  extended  to  bimodal  and  trimodal  mondisperse  propellants  (23,24). 

In  addition  to  the  above  assumptions,  radiation  heat  transfer  is 
neglected,  an  average  value  for  the  solid  and  the  gas  phase  specific 
heat  is  assumed,  and  the  mass  diffusion  coefficient  and  thermal 
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Figure  9.  Interaction  of  Primary  Flame  and  AP  Monopropellant 
Flame  at  Different  Pressures  (Dashed  Lines  indicate 
hypothetical  position) 
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conductivity  of  the  reacting  gases  is  averaged  over  the  gas  phase 
reaction  zone. 

Development  of  the  Basic  Equations.  In  the  development  of  the 
basic  equations,  the  intersection  diameter  of  the  statistically  average 
oxidizer  particle  at  the  burning  surface  is  (21) 


D.  = /TTS  Dq  20 

From  the  equation  of  continuity  for  a burning  propellant  the  total 
mass  flux  is 


■"t  - '■‘-'p  ■ “flSf/sp  ♦ V( 
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where  m is  the  mass  flux,  r is  burning  rate,  is  propellant  density 
and  S is  the  surface  area.  Subscripts  f,  ox  and  T represent  the  fuel, 
oxidizer  and  total  area,  respectively.  Equation  26  can  be  simplified 
further  by  requiring  that  stoichiometry  be  maintained  over  a long 
period  of  time;  i.e.,  the  ratio  of  oxidizer  consumption  to  fuel  corn- 
sumption  is  equal  to  the  weight  ratio  of  oxidizer  to  fuel  in  the 
propellant.  Thus 


m^  (S  /S,) 
ox^  ox  T'  . g 

ilif(S^/S^)'  ■ i - g 
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The  ratio  of  oxidizer  surface  area  at  the  burning  surface  to  the 
total  burning  surface  area  is  evaluated  from  geometric  relations  and 
from  employing  the  same  ignition  delay  concept  used  by  Hermance.  Thus 
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C[6{f)2  + 1] 

0 

[6;(^)  + 1] 
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where  c is  the  volume  fraction  of  oxidizer  in  the  propellant  and  h/D^ 
has  the  form 


JL  = 1 (1  + — ){i 
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The  ignition  delay  is  assumed  to  have  the  form  derived  by  Shannon  (25). 
Thus 


C.  d“  ^ 
t = .ign  0.- 
ign  pin 
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The  model  utilizes  Arrhenius  expressions  for  the  mass  burning  rate 
of  fuel  and  oxidizer.  Thus 

« 

m^  = A^exp[-E^/RTg]  31 

and 

"ox  ' ''ox'’''’[-Eox'"sJ 

where  A is  the  Arrhenius  frequency  factor,  E is  the  activation  energy 
R is  the  gas  constant  and  T^  is  the  average  surface  temperature. 

Since  the  oxidizer  regression  is  rate  controlling,  the  Arrhenius  ex- 
pression for  the  mass  burning  rate  of  fuel  is  used  only  in  a secondary 
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sense.  That  is,  it  is  used  for  the  determination  of  the  area  ratio 

A one-dimensional  energy  balance  at  the  burning  surface  as  depicted 
in  Figure  10  yields 

'"t^'p^^s  ■ ^0^  ""  ■"'ox^^ox'^^T^\  ■ 

+ (1  - ^ QfF?<p(-^’j)] 

where  Cp  is  the  average  specific  heat  of  the  gases  and  the  solid,  and 
is  the  initial  ;iropellant  temperature,  and  represent  the  energy 
required  to  vaporize  the  oxidizer  and  fuel  respectively,  Qpp,  Q^p,  and 
Qpp  represent  the  energy  released  in  the  final  flame,  and  AP  flame  and 
the  primary  flame  respectively,  Bp  is  the  fraction  of  reactants  that 

★ ★ •k 

react  in  the  primary  flame,  and  f.^p  and  ^pp  are  the  nondimensional 
flame  standoff  distances  for  the  final  flame,  AP  flame  and  the  pri- 
mary flame,  respectively. 

After  rearranging  equation  33  and  combining  with  equation  27,  the 
equation  for  the  average  surface  temperature  can  be  written 

- <'  - 

*(lJff/Cp)oxp(-£y)]  * ?^(Qpp/Cp)«p(-5pf)  3« 

The  AP  flame  heat  release  is  evaluated  from  the  energy  balance 
depicted  in  Figure  11.  Thus 

V - 'p<Vp  - ^o>  ’ 'll 
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where  T^p  Is  the  adiabatic  AP  flame  temperature. 

The  final  flame  heat  release  is  calculated  from  an  etieryy  balance 
as  depicted  in  Fiyure  12.  Thus 

QfF  = (Cp/a)[(T^  - T^)  - ct(T^p  - T^)  + (1  - .OQf/Cp]  36 

rthere  is  the  final  diffusion  flame  temperature. 

The  primary  flame  heat  release  is  evaluated  in  an  analoyoiis  mannei 
to  the  final  flame  heat  release.  Thus 

Finally,  the  nondimensional  yas  phase  flame  standoff  distances  are 
evaluated  by  consideriny  the  reactant  mixiny  lenyth  and  the  time  it 
takes  for  the  chemical  reaction  to  occur  after  mixing  has  taken  place. 

The  mixiny  lenyth  is  evaluated  by  calculating  the  peak  height  for 
a two  dimensional  diffusion  flame  and  then  trar\storminy  the  two  dimen- 
sional flame  into  a one  dimensional  flame  at  a height  above  the  surface 
which  is  some  fraction  of  the  peak  height  of  the  two  dimensional  flame. 

The  height  of  the  one  dimensional  flame  is  such  that  the  conductive 
heat  transfer  to  the  surface  is  approximately  the  same  for  the  one 
dimensional  flame  as  for  the  two  dimensional  flame. 

The  chemical  kinetic  reaction  lenyth  is  evaluated  by  consideriny 
the  reaction  to  take  place  in  a laminar  stream  moving  at  some  velocity 
in  which  a chemical  reaction  of  arbitrary  order  is  taking  place. 

The  AP  flame  is  premixed  and  is  therefore,  dependent  only  on  kinetics. 
The  AP  flame  standoff  distance  is  approximated  by 
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X - UT 

★ 

where  X is  the  standoff  distance,  u is  the  average  gas  stream  vcloc’ty 
and  T is  the  time  it  takes  for  the  reactants  to  react. 

For  an  arbitrary  reaction  order  the  reaction  time  can  be  written 


‘‘ap 


^ [K, 


(>5ap  - 0^-^ 
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where  K is  the  rate  constant  and  6^p  is  the  reaction  order.  Combining 
equations  4,  38  and  39  yields 


"'ox^'‘ap'’ 
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where  k^p  now  includes  the  proportional ity  constant  of  equation  39. 
In  nondimensional  form  the  AP  flame  standoff  distance  is 


^■AP  ’'aP^P%x^^ 
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where  \ is  the  average  thermal  conductivity  of  the  reactive  and  combus- 
tion gases. 

The  reaction  distance  for  the  primary  flai.'e  is  evaluated  in  an 
analogous  manner.  Thus 

* • ^PF  iio 

XpF  = m^/KppP  42 

Burke  and  Schumann  (26,27)  solved  the  differential  form  of  the  mass 
conservation  equation  for  Bunsen  burner  type  of  diffusion  flame.  Be- 
cause of  the  similarity  of  the  Bunsen  burner  to  the  burning  oxidizer 
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particle  in  a pyrolyzing  fuel /oxidizer  matrix,  the  Burke  Schumann 
analysis  was  used  in  the  model  to  determine  the  diffusion  mixing  length. 
The  calculated  diffusion  flame  height  represents  the  maximum  flame 
height  above  the  center  of  the  oxidizer  crystal.  The  actual  flame 
shape  was  not  calculated,  however,  a paraboloid  of  revolution  was 
assumed  for  the  flame  shape. 

To  avoid  the  complications  of  solving  ,=»  two  dimensional  heat  trans- 
fer problem,  an  average  effective  flame  height  was  defined,  The 

average  flame  height  factor,  is  that  fraction  of  the  calculated 

flame  height  which  corresponds  to  the  distance  a one  dimensional  flame 
would  be  from  the  surface  to  provide  the  same  energy  to  the  surface  as 
the  actual  flame.  Using  this  factor  the  effective  heat  transfer  dis- 
tance is 


*0  " *fh^0 
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In  the  primary  flame  the  components  must  mix  together  and  then 
react  so  that  the  standoff  distance  is  a diffusional  mixing  length 
followed  by  the  distance  that  it  takes  the  ingredients  to  react.  Thus 


Y -f  X 
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where  Xpp  is  the  effective  mixing  length  and  ^PR  is  the  reaction  length. 
In  nondimensional  form  the  primary  flame  standoff  distance  is 


Since  a final  diffusion  flame  exists  above  the  AP  flame,  the  final 


flame  standoff  distance  is  composed  of  the  AP  reaction  distance  plus 
the  effective  diffusion  length  for  the  mixing  of  the  AP  flame  products 
and  the  fuel  binder  pyrolysis  products.  Thus 


In  nondimcnsional  form,  the  final  flame  standoff  distance  is 


* 
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The  only  remaining  unknown  is  the  fraction  of  reactants  that 
react  in  the  primary  flame.  It  was  not  possible  . a rigorous  deter- 
mination of  the  diffusion  flame  shape  from  the  Burkt-Schumann  analysis. 
Thus,  a parabolic  shape  was  assumed  for  the  diffusion  flames. 

The  parameter  is  the  determined  by  assuming  that  the  flow  from 
the  oxidizer  suntace  is  perpendicular  to  the  planar  burning  surface. 

The  intersection  of  the  flat  AF  flame  with  the  parabolic  primary 
diffusion  flame  can  be  projected  onto  the  surface  of  the  oxidizer  and 
the  area  corresponding  to  each  flame  determined.  The  parameter  8^  may 
be  considered  as  that  fraction  of  the  AP  surface  whose  products  enter 
tl\e  primary  flame.  This  is  shown  schematical ly  in  Figure  13. 

The  pressure  dependence  of  0^  results  from  the  pressure  dependence 
of  the  primary  and  AP  flame  standoff  distances.  At  low  pressures  the 
combined  primary  flame  mixing  and  reacting  distances  are  less  than  the 
AP  flame  height.  Therefore,  all  of  the  oxidizing  species  react  directly 
with  the  binder  products  rather  than  with  the  ammonia  giving  a value 
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of  unity  for  By  increasing  the  pressure,  the  AP  flame  height  is 
reduced,  the  kinetic  distance  of  the  primary  flame  is  reduced  and  the 
primary  diffusion  distance  is  increased.  At  higher  pressures  the  dif- 
fusion flame  intersects  the  AP  flame.  Near  the  center  of  the  crystal 
where  the  diffusion  distance  is  greatest  the  oxidizing  species  will 
react  preferentially  with  the  ammonia.  Under  these  conditions  is 
between  zero  and  one. 

At  very  high  pressures  the  AP  flame  standoff  distance  and  the  pri- 
mary flame  reacting  distance  become  very  small  while  the  diffusion 
distance  becomes  very  large.  Under  these  conditions  a very  small  amount 
of  the  oxidizing  species  will  react  with  the  binder  products  in  the 
primary  flame  near  the  oxidizer  binder  inte'face  and  approaches 
zero.  The  above  processes  are  illustrated  in  figure  9. 

The  propellant  burning  rate  is  calculated  by  solving  equations 
26,  27,  28,  29,  32  and  34,  using  an  iteration  procedure. 

Discussion  of  the  BDP  Mod^l . The  purpose  of  the  BDP  model  was  to 
describe  the  combustion  process  of  composite  propellants  in  as  realistic 
a manner  as  possible  and  then  relate  the  effects  of  various  combustion 
trends  to  propellant  compositional  variations  and/or  canbustion  mechan- 
isms. Comparisons  were  made  with  a series  of  unimodal  AP-polysul f ide 
propellants.  The  predicted  dependence  of  burning  rate  on  oxidizer  con- 
centration was  essentially  the  same  as  observed  expe'-imental  ly , while 
the  dependence  of  the  burning  rate  on  the  oxidizer  particle  size  was 
somewhat  greater  than  observed  experimentally.  The  dependence  of  the 
burning  rate  on  the  initial  temperature  was  found  to  be  in  qualitative 
agreement  with  the  data. 
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The  original  BDP  model  applied  solely  to  an  additive  free  propel- 
lant with  spherical,  monodisperse  oxidizer.  Cohen,  Derr,  and  Price 
subsequently  extended  the  model  to  propellants  with  aluminum  and  bi- 
modal  oxidizer  (24).  Beckstead  (22,23)  has  since  extended  the  model  to 
trimodal  propellants.  However,  in  all  of  the  above  extensions  of  the 
model  several  limitations  remain.  The  limitations  will  be  discussed 
in  detail  below. 

The  first  limitation  occurs  in  the  calculation  of  the  diffusion 
flame  standoff  distances.  The  original  BOP  model  utilized  the  Burke- 
Schumann  (26)  diffusion  flame  analysis  for  the  calculation  of  the  dif- 
fusion flame  standoff  distances.  The  general  solution  for  the  position 
of  the  diffusion  flame  front  is  given  by  the  following  series  solution: 


2 / \ ^ 
v-c  (v+1 ) _ Y 

2(v+l  )c 


■j](')>jC)J^(>|).0 

4)ij2(^.) 


exp  [(1  -*'1  Mi{i^?)n/2i|/^]  48 


where  v is  a flame  stoichiometry  related  coefficient,  4'^  is  the  i^^ 
zero  of  a Bessel  function  of  the  first  kind,  c is  the  ratio  of  the 
oxidizer  diameter  at  the  surface  to  a diameter  associated  with  the  fuel 
binder,  ip  is  related  to  the  gas  diffusivity,  gas  velocity  and  surface 
geometry,  and  4 and  n are  nondimcnsional  radial  and  axial  coordinates, 
respectivley. 

In  the  BDP  model,  only  the  first  term  of  the  series  solution  is 
considered.  This  truncated  version  of  the  Burke-Schumann  analysis 
yields  the  following  expression  for  the  nondimensional  diffusion  flame 
height  n: 
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for  the  axial  nondimensional  coordinate  ^ equal  to  zero. 

Figure  14  depicts  the  nondimensional  standoff  distance  as  a function 
of  the  oxidizer  volume  fraction.  The  dashed  line  represents  the  exact 
solution  to  equation  48  while  the  solid  line  represents  the  approximate 
solution  (equation  49).  As  can  be  seen,  the  approximate  solution  is 
only  valid  for  over-ventilated  flames  near  stoichiometric  conditions. 

Real  composite  propellants,  on  the  other  hand,  contain  oxidizer 
mass  fractions  which  are  much  lower  than  the  stoichiometric  oxidizer 
mass  fraction.  Consequently,  substantial  errors  in  the  calculated 
diffusion  flame  standoff  distance  may  exist.  In  recent  work  with  the 
BDP  model  by  Beckstead  (22,23),  the  diffusion  flame  analysis  was  re- 
placed by  a one-dimensional  approximate  analysis  similar  to  that  of 
Sunmerfield  (16). 

Another  limitation  of  the  BDP  model  involves  the  assumption  of  the 

packing  structure  for  the  oxidizer  particles.  Figure  15  illustrates 

the  packing  assumed  for  the  BDP  model  where  each  box  contains  an  oxidizer 

particle  and  the  fuel  associated  with  that  oxidizer  particle.  Letting 

1 /3 

n be  the  number  of  particles  per  unit  volume,  then  n is  the  number 
of  particles  per  unit  length  for  the  assumed  packing  structure.  It  is 
clear  that 


6nV3  + = 1 
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OXIDIZER  VOLUME  FRACTION 


Figure  14.  Diffusion  Flame  Standoff  Distance  as  a Function  of 
Oxidizer  Volume  Fraction 
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where  6 is  the  distance  between  particles  at  the  burning  surface  and 
is  the  particle  intersection  diameter  at  the  burning  surface.  There- 


6 = 1/n^/^  - 


The  number  of  particles  per  unit  volume  is 


n - 

VftY 


where  ; is  the  volume  fraction  of  oxidizer  and  is  the  volume  of  an 
oxidizer  crystal . Thus 


n - 6U(rrD^)  53 

Combining  equations  20,  51  and  53  yields  the  following  equation  for 


6 = ' Do  - .'JTJ  Do 


Since 


b = I + -^  55 

Combining  51  and  52  yields 

Since  the  packing  structure  of  composite  propellants  is  not  ordered, 
but  rather  is  a completely  random  structure,  equation  56  for  b is  in- 
correct, for  it  describes  an  ordered  packing  of  particles. 

This  is  further  illustrated  in  Figure  16  where  6/0^  is  plotted  as 
a function  of  oxidizer  volume  fraction.  Clearly,  for  oxidizer  volume 
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5. 


Figure  16.  The  Ratio  of  the  Distance  3etween  Oxidizer  Particles 
to  the  Particle  Diameter  as  a Function  of  Oxidizer 
Volume  Fraction 
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fractions  approaching  unity,  the  ratio  (S/D^  should  approach  zero.  How- 
ever, for  the  BDP  method,  the  ratio  of  5/0^  actually  goes  negative  at 
a volume  fraction  0.78.  This  is  physically  unrealistic. 

The  final  limitation  deals  with  the  determination  of  the  mean  state 
of  an  oxidizer  particle  undergoing  combustion  at  the  burning  surface. 

In  the  BOP  approach,  the  burning  surface  is  generated  mathematical ly  by 
passing  an  imaginary  plane  through  a randomly  packed  monodispersion  of 
spherical  oxidizer  particles.  The  mean  state  in  terms  of  the  diameter 
of  particles  at  the  burning  surface  is  assumed  to  be  Die  average  inter- 
section diameter  of  the  particles  with  the  imaginary  plane.  Since  the 
particles  are  assumed  to  undergo  an  ignition  delay,  the  mean  state  in 
terms  of  the  diameter  for  oxidizer  crystals  described  above  is  incor- 
rect. A more  detailed  explanation  along  with  a description  of  a 
better  determination  of  tlie  mean  state  diameter  of  the  oxidizer  ptir- 
ticle  appears  in  a subsequent  section. 

Hescrjption  of  frosivc  Burning 

Requirements  for  ever  higher  propellant  loading  in  solid  propellant 
rocket  motors  have  led  to  the  development  of  grain  configurations  with 
relatively  low  port-to- throat  area  ratios.  This,  in  turn,  results  in 
high  velocities  of  propellant  combustion  gases  flowing  across  the  burn- 
ing propellant  surfaces  in  the  aft  region  of  the  propellant  grains.  The 
presence  of  a high  cross  flow  velocity  over  a propellant  burning  surface 
.causes  enhancement  of  the  propellant  burning  rate  and  this  enhancement 
of  the  burning  rate  is  termed  erosive  burning,  for  higlt  cross  flow 
velocities  the  burning  rate  can  be  substatitial ly  higher  than  the  burning 
rate  of  the  same  propellants  not  subject  to  a cross  flow  velocity. 
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Moreover,  the  effects  of  erosive  burning  are  critical  in  that  the  ero- 
sive burning  rate  contributions  strongly  influence  performance  level  and 
performance  repeatability  of  propulsion  systems. 

Since  erosive  burning  is  present  in  many  solid  propellant  propulsion 
systems  and  since  there  is  such  a strong  interaction  between  the  local 
flow  environment  and  the  propellant  burning  rate,  it  is  necessary  to  be 
able  to  predict  this  interaction  in  oi^der  to  design  and  calculate  the 
performance  of  a low  port/throat  area  ratio  rocket.  At  present  there 
is  no  unifying  model  which  can  be  used  to  predict  the  propellant  burning 
rates  under  erosive  conditions. 

Experimental  Observations  of  Erosive  Burning 

There  are  several  general  observations  of  importance  which  can  be 
attributed  to  past  experimental  studies  (28-36).  A few  of  the  obser- 
vations follow: 

1.  Threshold  velocities  are  usually  observed,  i.e.,  burning  rate 
is  not  affected  until  the  core  gas  reaches  a certain  velocity. 

2.  Slower  burning  propellants  are  more  strongly  affected  by  cross 
flows  than  higher  burning  rate  formulations. 

3.  At  high  pressure,  the  burning  rate  under  erosive  conditions 
tends  to  approach  the  same  value  for  all  propellants  (at  the  same  flow 
velocity)  regardless  of  the  burning  rate  of  propellants  at  zero  cross 
flow. 

4.  Erosive  burning  rates  do  not  depend  upon  the  gas  temperature  of 
the  cross  flow. 

Most  of  the  above  observations  were  made  in  systems  with  cross  flow 
velocities  having  a Mach  number  less  than  0.3. 


Erosive  Burning  Models 


A large  number  of  erosive  burning  models  have  appeared  in  the  liter- 
ature (37-46).  Of  those  models,  all  fall  into  one  of  four  basic 
categories.  They  are: 

1.  Models  based  on  heat  transfer  from  a core  gas  flowing  parallel 
to  the  burning  surface. 

2.  Models  based  on  the  alteration  of  transport  properties  in  the 
region  of  the  gas  phase  reaction  zone. 

3.  Models  based  on  chemically  reacting  boundary  layer  theory. 

4.  Models  based  on  displacement  of  the  flame  by  the  cross  flow. 

The  only  existing  model  of  the  third  type  (45)  assumes  a homogeneous 

propellant  system.  Therefore,  since  composite  propellants  are  far  from 
homogeneous,  no  further  discussion  will  be  offered  for  this  class  of 
model.  The  model  representative  of  the  1st  type  is  the  Lenoir  and 
Robillard  (37)  model.  Lengell^'s  (43)  erosive  burning  model  is  repre- 
sentative of  the  type  two  models,  and  the  erosive  burning  model  of  King 
(46)  is  the  only  model  of  the  fourth  type  in  existence.  A discussion 
of  those  three  models  follows. 

The  Lenoir  Robillard  Model 

Lenoir  and  Robillard  proposed  an  erosive  burning  rate  law  for  solid 
propellants  based  on  the  postulate  that  the  erosive  component  of  the 
burning  rate  is  proportional  to  the  heat  transfer  rate  from  the  core 
of  the  combustion  gas  stream  to  the  propellant  burning  surface.  A 
transpiration  model  for  the  core  gas  was  used  to  compute  the  heat  trans- 
fer rate. 
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In  the  development  of  the  basic  equations,  it  was  assumed  that  the 
burning  rate  was  dependent  upon  the  amount  of  heat  the  propellant  re- 
ceives. The  main  source  of  heat  was  assumed  to  be  from  the  primary 
burning  zone  in  the  first  few  microns  adjacent  to  the  burning  surface. 
Furthermore,  the  heat  transfer  from  the  primary  combustion  zone  was 
assumed  to  be  independent  of  the  core  gas  velocity.  The  pressure  de- 
pendence of  burning  rate  was  assumed  to  be  controlled  by  the  pressure 
dependence  of  the  location  of  the  primary  combustion  zone  above  the 
burning  surface.  As  pressure  increases,  the  flame  becomes  nearer  to  the 
surface,  heat  flow  to  the  burning  surface  from  the  flame  increases  and 
thus  the  burning  rate  increases.  A second  source  of  heat  was  assumed 
to  be  from  the  core  of  the  hot  combustion  gases.  This  second  source 
of  heat  was  assumed  to  act  through  a convective  mechanism  and  was  there- 
fore dependent  upon  the  cross  flow  gas  flow  rate  . The  total  burning 
rate  was  assumed  to  be  the  sum  of  the  two  effects,  a rate  dependent  upon 
pressure  and  an  erosive  rate  dependent  upon  combustion  gas  cross  flow 
rate.  Thus, 

^ = '■p  * »'e  57 

The  pressure  dependent  part  of  the  burning  rate  was  assumed  to 
correlate  to  the  well  known  linear  burning  rate  law 

r = cP"  58 

P 

where  c and  n are  constants  for  a particular  propellant  and  are  inde- 
pendent of  combustion  pressure. 
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The  combustion  gas  cross  flow  velocity  dependent  portion  of  the 
burning  rate  was  derived  by  assuming  the  erosive  burning  rate  to  be 
proportional  to  the  heat  transfer  coefficient.  Thus 

r.  = kh  58A 

e 

where  h is  the  heat  transfer  coefficient  and  k is  the  proportionality 
constant.  The  relationship  between  the  heat  transfer  coefficient,  h, 
under  conditions  of  transpiration,  to  the  heat  transfer  coefficient  in 
the  absence  of  transpiration  was  expressed  in  exponential  form  as 

h = h^cxp(-3Q/G)  58B 

where  8 is  a dimensionless  constant,  Q is  the  mass  velocity  of  the 
transpiring  gas  and  6 is  the  mass  velocity  of  the  core  gas.  The  heat 
transfer  coefficient  under  conditions  of  zero  transpiration  was  assumed 
to  correlate  by  the  Chilton  Colburn  equation  of  flow  over  a flat  plate. 

^ ^ to. 0288  58C 

where  Cp  is  the  specific  heat  of  the  combustion  gases,  u is  the  visco- 
sity of  the  combustion  gases,  L is  the  length  from  the  head  of  the 
grain,  and  Pj.  is  the  Prandtl  number  of  the  combustion  gases.  The  ex- 
pression in  the  brackets  of  equation  58C  was  assumed  to  be  insensitive 
to  changes  in  propellant  identity  so  that  the  expression  in  the  brackets 
could  be  combined  with  the  proportionality  constant  of  equation  58A  to 
form  a constant, 
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58D 


a = C0.0288C„P°-^P'°'^^^]k 
p r 

Combining  equations  57,  58,  58A,  58B,  and  580,  resulted  in  the 
following  expression  for  burning  rate: 

r = CP"  + exp  (-8Q/G)  58E 

The  Lenoir  Robillard  theory  has  been  compared  with  experimental 
data  and  satisfactory  agreement  indicated.  However,  the  agreement  was 
obtained  by  matching  equation  58E  to  the  experimental  data  in  order  to 
determine  the  constants  a and  f..  Thus,  the  Lenoir  Robillard  theory  is 
empirical  in  nature.  Furthermore,  the  model  predicts  a dependence  of 
burning  rate  on  the  mainstream  gas  temperature.  Later  experimental 
work  by  Marklund  and  Lake  (33)  indicated  no  such  dependence.  The 
Lenoir  Robillard  model  has,  however,  served  the  test  of  time  and  is 
probably  the  most  widely  used  model  in  the  solid  propellant  propulsion 
community. 


The  Flame  Displacement  Model 

The  erosive  burning  model  of  King  (46)  postulates  the  mechanism  for 
erosive  burning  to  bo  an  enhanced  heat  transfer  to  the  burning  surface. 
The  enhanced  heat  transfer  results  from  the  displacement  of  the  diffusion 
flames  by  the  cross  flow  of  combustion  gases. 

The  burning  rate  is  hypothesized  to  be  dependent  upon  the  heat 
feedback  from  two  flame  zones.  They  are  the  AP  monopropellant  flame 
and  a diffusion  flame  between  the  binder  pyrolysis  products  and  the 
products  of  the  AP  monopropel lant  flame.  Thus,  the  burning  rate  is 
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dependent  upon  three  distance  parameters  associated  with  the  two  flames. 
The  three  distance  parameters  considered  are  the  distance  from  the  pro- 
pellant surface  to  the  average  location  of  the  kinetically  controlled 
AP  monopropellant  heat  release,  the  distance  associated  with  mixing 
of  the  oxidizer  and  fuel  for  the  diffusion  flame,  and  the  distance 
associated  with  the  fuel  oxidizer  reaction  time  subsequent  to  mixing. 
Those  distances  are  then  assumed  to  be  modified  in  the  presence  of  a 
cross  flow  velocity  by  the  following  equation: 


L = Lq  A-tn6 
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where  is  the  flame  standoff  distance  in  the  absence  of  a cross  flow 
and  0 represents  the  angle  of  the  average  flow  vector  in  the  mixing 
region. 

This  analysis  was  combined  with  an  analysis  of  the  boundary  layer 
flow  (which  gives  the  cross  flow  velocity  as  a function  of  distance 
from  the  propellant  surface,  mainstream  velocity,  and  propellant  burning 
rate)  to  permit  the  calculation  of  the  angle  of  the  average  flow 
vector  in  the  mixing  region.  The  resulting  equation  for  burning  rate  in 
the  absence  of  cross  flow  is 


r = A3P  1 
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where  P is  the  pressure  and  ^ is  the  average  oxidizer  particle  size  in 
the  propellant.  A regression  analysis  using  no-cross  flow  burning  rate 
data  must  be  performed  to  obtain  best  fit  values  of  the  constants  A3, 

A^  and  Ag.  The  product  of  Ag  with  appearing  in  burning  rate  expres- 
sion may  be  lumped  into  one  parameter  during  the  regression  analysis 
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such  that  an  average  oxidizer  particle  size  need  not  be  defined. 

The  data  of  Mickley  and  Davis  (47)  were  used  to  develop  empirical 

expressions  for  the  local  cross  flow  velocity  as  a function  of  distance 

from  the  propellant  burning  surface,  mainstream  cross  flow  velocity, 

and  transpiration  rate.  The  resulting  analysis  yielded  eight  equations 

and  eight  unknowns  for  the  burning  rate  of  a given  composite  propellant 

at  a given  pressure  and  cross  flow  velocity.  A computer  code  was  used 

to  solve  this  system  of  equations  yielding  a predicted  burning  rate  for 

a given  pressure,  cross  flow  velocity  and  set  of  constant  A^,  A^  and 
—2 

A^D  . The  model  employs  no  empirical  constants  other  than  those  from 
the  analysis  of  the  no-cross  flovr  burning  rate  data. 

Discussion  of  the  Flame  Displacement  Model.  King  has  hypothesized 
that  the  flame  displacement  is  due  to  the  cross  flow  induced  distortion 
of  the  flow  field  in  the  region  of  the  gas  phase  reaction  zone.  In- 
herent in  the  model  is  the  assumption  that  transport  properties  are  not 
affected  by  the  cross  flow.  Thus,  the  characteristic  kinetic  reaction 
times  for  the  gas  phase  reactions  taking  place  above  the  burning  sur- 
face are  the  same  for  both  the  case  of  zero  cross  flow  and  the  case  where 
a strong  cross  flow  is  present.  However,  since  the  flame  is  assumed  to 
be  conical  in  shape,  the  effective  distance  for  diffusional  mixing  is 
decreased  in  the  presence  of  a cross  flow  velocity.  The  reduction  in 
the  diffusional  mixing  length  is  due  to  a bending  of  the  flame  such  that 
its  central  axis  is  pointed  in  a direction  parallel  to  the  average  velocity 
vector  for  the  flow  of  gases  near  the  burning  surface.  Thus,  the  flame  is 
"squeezed"  such  that  the  radial  distance  from  the  central  axis  of  the  cone 
to  the  outsioe  surface  of  the  cone  is  reduced.  Subseguent ly . the  mixing 
length  is  reduced.  Thus,  the  cross  flow  vdocltv  Itidu'  -d  r ^h  nt  of 
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the  burning  rate  is  realized  by  a decrease  in  the  diffusion  flame  stand- 
off distance.  This  effect  is  exactly  the  same  as  the  effect  due  to  the 
enhancement  of  the  diffusion  coefficient  in  the  gas  phase  reaction  zone 
by  the  presence  of  a turbulent  boundary  layer.  However,  in  this  case 
the  effect  of  the  turbulence  on  the  thermal  conductivity  in  the  gas 
phase  reaction  zone  was  assumed  negligible. 

Lengelle's  Erosive  Burning  Model 

The  model  is  based  on  the  combination  of  Summerfield's  (16)  gran- 
ular diffusion  flame  combustion  model  with  a model  for  the  turbulent 
boundary  layer  on  a flat  plate  with  injection.  The  basic  assumption 
adopted  was  that  of  a turbulent  boundary  layer  developing  above  the 
propellant  surface.  The  transport  coefficients  of  diffusion  and  heat 
conduction  are  then  modified  in  such  a way  that  the  heat  flux  to  the 
surface,  and  consequently  its  burning  rate,  can  be  enhanced. 

The  turbulent  boundary  layer  model  used  in  the  analysis  assumes 

incompressible  flow  over  a flat  plate  with  injection.  The  continuity 

9 (i 

and  the  momentum  conservation  equations  were  simplified  by  letting  ^ 
and  u be  negligible  where  u is  the  gas  velocity  parallel  to  the 
burning  surface.  The  resulting  approximate  relationship  for  the  velocity 
profile  is 

Em  (1  + B<i.)  = n"  Em  (1  + B)  61 

where  B is  the  blowing  parameter,  <t'  is  the  ratio  of  downstream  velocity 
at  position  y in  the  boundary  layer  to  free  stream  velocity,  n is  the 
ratio  of  y distance  to  the  boundary  layer  thickness  and  n is  the  well 
known  empirical  velocity  profile  exponent  for  the  case  of  no  injection. 
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The  implicit  relationship  for  the  velocity  profile  given  by  equation 
61  was  awkward  for  obtaining  analytical  relations  for  the  turbulent 
boundary  layer  parameters.  Therefore,  the  implicit  relationship  was 
approximated  by  the  following  explicit  expression 

= . a = In  {}  *B)  62 

The  mixing  length  hypothesis  of  Prandtl  was  combined  with  the 
velocity  profile  given  by  equation  62  resulting  in  the  following  ex- 
pression for  e,  the  coefficient  of  turbulent  diffusivity  of  momentum: 


= .16Re.  an 
U <5 


a+1 
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The  average  value  of  the  coe^'ficient  for  turbulent  diffusivity  of 
momentum  in  the  flame  zone  is  given  by 


(QL)  = 1 

' u 'ave  L 


p £ dy  » .IBRej  (j 


a +1 
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where  L is  the  flame  zone  thickness  and  Re^  is  evaluated  from  an  approx- 
imate expression  which  relates  the  ratio  of  the  boundary  layer  thickness 
in  the  case  of  blowing  and  the  boundary  layer  thickness  in  the  case  of 
no  blowing,  to  the  reduction  of  the  wall  shear  stress.  Thus 
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where 


^ = C'a  (Ha) 
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Lengelle  also  assumed  that  the  turbulent  Lewis  and  Schmidt  numbers 
are  close  to  unity.  Then,  the  thermal  conductivity  and  mass  diffusion 
coefficient  of  the  combustion  gases  can  be  related  to  the  viscosity. 

Thus 

67 
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and 


A = CgpD 


pD  = ud  + pe/u) 
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where  X is  the  thermal  conductivity,  D is  the  mass  diffusion  coefficient 
and  u is  the  viscosity.  Utilizing  equations  67  and  68,  the  transport 
parameters  in  the  GDF  model  were  expressed  in  terms  of  the  coefficient 
for  turbulent  diffusivity  which  in  turn  could  be  evaluated  using 
equation  63.  The  resulting  expression  for  burning  rate  is 
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where 
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Pj  is  the  core  gas  density,  u^  is  the  core  gas  velocity  and  the  velo- 
city profile  exponent  has  been  approximated  by 


n = .52Re 


•0.1 
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The  predicted  results  of  the  model  compared  favorably  to  the 
experimental  results  of  Marklund  and  Lake.  The  appearance  of  a threshold 
velocity,  below  which  the  normal  burning  rate  is  not  effected,  was 
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determined  to  be  due  to  the  blocking  effect  of  the  transpiring  gases. 

The  well  known  observation  that  high  burning  rate  propellants  give  rise 
to  high  threshold  velocities  and  low  burning  rate  propellants  give  rise 
to  low  threshold  velocities,  was  also  explained  on  the  basis  of  the 
blowing  parameter. 

It  should  be  noted  that,  although  the  philosophy  behind  the  ero- 
sive burning  mechanism  of  Lengellf's  model  is  sound,  the  model  is  little 
more  than  a curve  fitting  device  due  to  the  empiricism  of  the  GDF  model. 
It  should  also  be  noted  that  the  model  is  in  serious  error  due  to  the 
assumed  independence  between  burning  rate  and  the  blowing  coefficient. 

In  order  for  the  model  tc  be  correct,  the  value  of  the  burning  rate 
must  be  determined  by  an  iteration  process  accounting  for  the  change  in 
the  blowing  parameter  with  changes  in  the  burning  rate.  Thus,  the  con- 
clusion by  Lengelle  that  the  erosive  burning  effect  is  additive  (i.e., 
the  overall  burning  rate  is  the  sum  of  a nonerosive  rate  and  a contri- 
bution due  to  erosive  burning)  is  incorrect. 

Nonsteady  State  Combustion 

In  evaluating  a propellant's  combustion  stability  characteristics, 
there  are  two  parameters  of  importance:  pressure  coupled  response  and 
velocity  coupled  response.  The  pressure  coupled  response  represents  the 
coupling  mechanism  between  the  gaseous  pressure  oscillations  and  the 
combustion  processes,  while  the  velocity  coupled  response  represents  the 
coupling  mechanism  between  the  cross  flow  velocity  transients  and  the 
combustion  processes.  The  interaction  occurs  principally  in  a relatively 
thin  region  near  the  burning  surface. 
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Following  traditional  acoustics  practice  it  has  usually  been 
assumed  for  the  pressure  coupling  that  the  fluctuation  of  velocity,  v‘, 
normal  to  the  surface  is  proportional  to  the  fluctuation  of  pressure  p' 
at  the  burning  surface.  The  coefficient  of  proportionality,  in  nor- 
malized form,  is  the  admittance  function,  A|j,  for  the  burning  surface. 
Thus 

\ M 

where  M is  the  average  Mach  number  of  the  flow  leaving  the  surface,  p 
is  the  gas  density,  p is  pressure,  y is  the  ratio  of  specific  heats  and 
m is  mass  flux.  The  primes  denote  fluctuating  values  and  the  super- 
script bars  denote  average  values. 

The  admittance  function  is,  in  general,  a complex  quantity  depending 
on  the  properties  of  the  materials  involved  and  the  frequency  of  the 
oscillations.  The  real  part  of  th''  admittance  function  gives  that  por- 
tion of  the  normal  velocity  fluctuation  which  is  in  phase  with  the 
pressure.  Thus  the  attenuation  or  growth  constant  for  steady  waves  has 
a part  proportional  to  the  real  part  of  the  admittance  function  such  that 
the  waves  are  driven  if  the  real  part  of  the  admittance  function  is 
positive.  A larger  value  of  the  real  part  of  the  admittance  function 
implies  a greater  tendency  for  combustion  to  drive  the  waves. 

A similar  interpretation  of  the  admittance  function  arises  in  con- 
nection with  the  velocity  fluctuation  u'  parallel  to  the  burning  surface. 
Due  to  its  possible  erosive  influence  on  the  burning  rate  the  velocity 
fluctuation  parallel  to  the  surface  causes  a fluctuation  of  velocity 
normal  to  the  burning  surface. 


(m'/m) 

(p'/p) 

(p7yF) 

(pVyf) 
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If  one  assumes  that  the  oscillations  are  isentropic  the  second 
term  in  equation  72  is  unity  and  the  admittance  function  becomes  de- 
pendent only  on  the  value  of  the  first  term.  The  first  term,  divided 
by  the  product  (M  y),  is  usually  referred  to  as  the  response  function 
of  the  propellant  and  is  given  the  symbol  R.  The  real  part  of  the  re- 
sponse function  has  the  same  characteristics  as  the  real  part  of  the 
admittance  function.  Thus,  waves  are  driven  if  the  real  part  of  the 
response  function  is  greater  than  1/  y and  a larger  value  of  the  res- 
sponse  function  implies  a greater  tendency  for  the  combustion  process 
to  drive  the  waves. 

The  pressure  coupled  response  is 


R 


P ' TpV  pT 
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and  the  velocity  coupled  response  is 
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Nonsteady  State  Models 

It  is  the  objective  of  the  nonsteady  state  combustion  models  to 
predict  the  effect  of  propellant  composition  on  the  propellant  response 
function.  In  the  past,  models  attempting  to  describe  the  nonsteady  com- 
bustion of  propellants  (48)  have,  in  general,  addressed  only  homogeneous 
propellants.  Although  the  homogeneous  models  have  been  applied  to 
heterogeneous  propellants  with  a limited  degree  of  success,  there  is  a 
formidable  deficiency  in  the  models.  That  is,  the  models  are  entirely 
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unable  to  predict  the  effect  of  oxidizer  particle  size  and  particle 
size  distribution  on  the  pressure  coupled  response. 

Several  of  the  more  pertinent  models  which  have  appeared  in  the 
literature  are:  the  Denison  and  Baum  (DB)  model  (49),  the  Cohen  (C) 
iiwdel  (51),  the  Zel'dovich  Novozhilov  (ZN)  model  (50)  and  the  velocity 
coupling  model  of  Lengelld  (43).  A discussion  of  those  models  follows. 

The  Denison  and  Baum  Model 

The  DB  model  was  an  early  attempt  to  model  the  pressure  coupled 
response  of  propellants.  The  model  is  based  upon  the  solution  of  the 
nonsteady  energy  equation  for  the  solid  phase. 

The  model  assumes  a quasi-steady  gas  phase  with  harmonic  pressure 
oscillations  and  a planar  flame  front. 

The  energy  equation  for  the  solid  phase  is,  for  the  case  of  no  heat 
generation, 

X„(0^/;^x^)  - mC,  (DT/hx)  - P C,  (3T/3t)  « 0 ^5 

P P P P 

The  crucial  boundary  condition  is  the  matching  of  the  heat  flux  at 
the  solid/gas  phase  interface. 

iho  resultirj  expression  for  the  pressure  coupled  response  is 

K %r  ^ ^^^i^  • A f (A/xr  V (T+7iyT  SB 

•here  and  are  respectively  the  real  and  imaginary  parts  of  the 
-.urr  coupled  response  normalized  with  respect  to  the  burning  rate 

irr  exponent.  The  parameter,  X , is  a complex  function  of  the 
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nond imensional  frequency,  fi.  The  real  part  of  X is 


V = 2 ^ (!//?)(!  + t 77 

while  the  imaqinary  part  of  A is 

= [I/?  - 1]’^^  78 

where  the  nond imensional  frequency  can  be  written 

" = Pp  wAp  79 

The  constants  A and  B are  related  to  propellant  properties.  The 
parameter  A is  related  to  the  surface  pyrolysis 

A = E^(l  - Vr'5)/Rr^  80 

where  is  the  activation  enerqy  associated  with  the  surface  pyrolysis, 
T^  is  the  initial  propellant  temperature,  is  the  average  surface 
temperature  and  R is  the  gas  constant. 

The  parameter  8 is  dependent  on  the  temperature  sensitivity  of  the 
propellant.  Thus 

B • Up/tp/jd  . Tp/ip  31 

where  j is  the  initial  temperature  index  in  the  burning  rate  law 

r ■ Cp^T^  82 
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The  model  must  be  classified  empirical,  at  most,  since  the  A and 
B parameters  are  adjusted  for  the  fitting  of  experimental  data.  Fur- 
thermore, the  model  has  not  met  with  a great  deal  of  success  at  fitting 
the  experimental  data  of  composite  propellants.  This  is  probably  due 
to  the  homogeneous  nature  of  the  model . 


The  Cohen  Model 

The  most  complete  attempt  at  accounting  for  the  oxidizer  particle 
size  effects  on  the  pressure  coupled  response  was  developed  by  Cohen 
et  al.  ( 1 ).  Cohen  postulated  that  the  combustion  response  will  peak 
at  a frequency  given  by  some  characteristic  transit  time  for  an  oxidizer 
particle  in  the  condensed  phase.  Thus 


where  fp  is  the  frequency  at  which  the  response  function  for  a parti- 
cular mode  peaks,  r is  the  propellant  burning  rate  and  D is  the  average 
particle  diameter  for  the  mode. 

Furthermo'e,  Cohen  postulated  that  the  magnitude  of  the  peak  re- 
sponse is  proportional  to  the  concentration  of  oxidizer  in  the  propellant. 


n 0 


where  Rpp/n  is  the  normalized  peak  value  of  the  pressure  coupled  re- 
sponse and  a is  the  mass  fraction  of  oxidizer  in  the  particular  mode. 
Cohen  assumed  that  the  response  function  is  only  a function  of 


propellant  type  and  may  be  obtained  from  data  for  a representative 
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propellant  type.  Propellant  data  for  A-13  were  used  to  construct  a 
calibration  response  function  which  could  then  be  used  to  obtain  the 
response  function  for  each  weight  mean  diameter  particle  size  distri- 
bution mode  in  a propellant  of  interest.  The  response  function  of  each 
mode  in  the  propellant  was  then  superimposed  to  obtain  the  total  pro- 
pellant response  function.  This  method  has  potential  but  is  still 
deficient  in  that  it  treats  a multimodal  polydisperse  propellant  only 
as  being  a multimodal  monodisperse  propellant. 

The  Zel'dovich  Novozhilov  Model 

The  ZN  approach  to  nonsteady  modeling  employs  the  same  assumptions 
as  those  of  the  DB  model.  However,  the  ZN  technique  allows  the  pseudo- 
propellant  response  functions  to  be  related  to  steady  state  propellant 
combustion  parameters  in  a more  realistic  manner. 

In  the  derivation  of  the  ZN  pressure  coupled  response  the  propel- 
lant mass  flux,  m,  can  be  written  as  a function  of  burning  surface 
temperature,  T^,  combustion  pressure,  p,  temperature  gradient  at  the 
burning  surface  in  the  condensed  phase,  f,  temperature  gradient  at  the 
condensed/gas  phase  interface  in  the  gas  phase,  f.  and  the  surface  heat 
release,  Q^.  Thus 

m = m(T^,p,f ,f ' ,Q^)  B5 

Similarly  for  the  surface  heat  release. 

Q = Q5(T^.P.f.f') 
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Conservation  of  energy  at  the  condensed/gas  phase  interface  yields 

Kf  = K'f  + ni  87 

where  k and  k'  are  the  thermal  conductivity  of  the  solid  and  gas  repec- 
tively. 

The  mass  flux  in  the  gas  phase,  m' , can  be  characterized  by  the 
surface  temperature,  the  temperature  gradient  in  the  gas  phase,  the 
heat  release  in  the  gas  phase,  Q^,  and  the  combustion  pressure. 

The  heat  release  in  the  gas  phase  can  be  characterized  by  the  surface 
temperature,  the  temperature  gradient  in  the  gas  phase  and  the  combustion 
pressure.  Thus  the  mass  flux  and  heat  release  can  be  written 


m = m'(T^,f‘ .Q’ ,p) 

88 

Q'  = Q’(T^.f',p) 

89 

Since  the  gas  phase  is  assumed  quasi-steady,  mass  continuity  yields 


Thus,  there  are  six  equations  for  seven  unt<nowns:  f.  f,  T^,  Q^, 

Q' , m and  m'  (the  pressure  is  assumed  to  be  a known  function  of  time). 

’ Therefore,  it  is  possible  to  express  any  value  as  a function  of  pressure 

and  the  temperature  gradient  in  the  condensed  phase  at  the  burning 
surface.  Specifically, 

m = m (f,p) 

I ’ Tj  = 
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Thus, 


. /Djflv  df  . ^ 

dt  ” p dt  '3p^^  dt 


93 


and 


dT,  3T^  . 

dt  ' 3f  ^p  dt  ' 9p  dt 
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Nondimensional izing,  and  taking  a first  order  approximation  for  the  per- 
turbations yields 


/ \ A J.  ^ \ M 
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where  and  41  are  the  nondimensional  perturbation  amplitudes 

in  burning  rate,  temperature  gradient,  combustion  pressure,  and  surface 
temperature  respectively. 

Equations  (95)  and  (96)  provide  the  means  for  obtaining  the  desired 
functional  relation  between  the  burning  rate  perturbation  amplitude  and 
the  pressure  perturbation  amplitude. 

With  the  aid  of  Jacobians  and  an  energy  balance  in  the  condensed 
phase  at  the  burning  surface,  the  partial  derivatives  appearing  on  the 
right  hand  side  of  equations  (95)  and  (96)  can  be  expressed  in  terms  of 
the  following  steady  state  parameters: 
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where  the  superscript  o indicates  steady  state,  is  the  surface  temper- 
ature and  Tq  is  the  initial  propellant  temperature. 

The  remaining  needed  relation  between  the  nondimensional  burning 
rate,  the  nondimensional  surface  temperature,  the  nondimensional  temper- 
ature gradient,  and  the  nondimensional  pressure  perturbation  amplitudes 
can  be  determined  by  solving  the  linearized  energy  equation  in  the  solid 
phase.  The  result  is 


\1)^(1+1/Z^)  - + iv^/n  = 0 


where  is  a complex  function  of  nondimensional  frequency  given  by 

1 1 o 1/2  12  1/2  1/2 

Z^  = ^ (16f?Sl)  - 1]  + i (IKl  +1)  -1]}  102 


The  desired  relation  between  the  burning  rate  perturbation  ampli- 


tude and  the  pressure  perturbation  amplitude  can  now  be  written.  Thus 
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n + (nv  - pu)Z^ 

'*1  " 1 - u + (v  - iu/n)Z^ 
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Assuming  an  isentropic  fluctuation  in  the  combustion  pressure,  the 
response  function  is  defined 


R 


0 p, 

r ^1 
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where  r^  and  p^  are  the  fluctuation  amplitude  of  burning  rate  and  pres- 
sure respectively.  Thus 


r^  = v^r 


Pi  = n^P 
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Therefore,  the  pressure  coupled  response  for  the  ZN  technique  is 


R_  = 


n + (nv  - pu)Z^ 
T~-  u + (v  - iu/fl)Zi 


107 


Equation  107  can  be  used  to  determine  a pressure  coupled  response 
function  by  evaluating  the  steady  state  parameters  for  the  propellant. 

The  problems  with  this  approach  are  that  some  of  the  steady  state  para- 
meters can  not  be  evaluated  experimentally  (specifically  those  involving 
the  surface  temperature)  and  the  assumptions  employed  in  the  model  break 
down  for  heterogeneous  propellants  when  the  characteristic  thermal 
wave  penetration  depth  is  not  much  larger  than  the  scale  of  heterogeneity 
in  the  solid  phase. 
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Lengelle's  Model 

Lengelle  formulated  a nonsteady  combustion  model  for  velocity 
coupling  based  on  the  linear  perturbation  of  an  erosive  burning  model 
based  upon  the  GDF  combustion  model.  The  erosive  burning  model  was 
discussed  in  a prior  section  of  this  review  and  will  not  be  discussed 
here. 

The  gas  phase  reaction  zone  was  assumed  to  respond  quasi-steadily 
to  external  oscillations  so  that  only  the  condensed  phase  introduced 
inertia  into  the  process. 

The  energy  equation  for  the  condensed  phase  with  harnic  lic  pertur- 
bations was  solved  for  the  heat  flux  penetrating  into  the  solid  phase. 
Thus 


ds  = Cpm  [(T^ 


T ) + T'  + -*(T  - 

o'  s m ' s 


T„)/s] 


108 


where  is  the  average  surface  temperature,  T^  is  the  initial  propel- 
lant temperature,  T^  is  the  amplitude  of  perturbation  of  the  surface 
temperature  and  Cp  is  the  propellant  specific  heat.  The  parameter  s is 


s = [1  + (1  + 4in)^/^]/2 
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where  fi  is  the  nondimensional  frequency  and  i designates  an  imaginary 
number. 

From  the  GDF  model,  with  transport  properties  modified  to  account 
for  turbulence,  the  following  may  be  written: 
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q+  = + pe)/L 

s y ' 

L = + pe)  in 

0"  ■ ‘>“a  I-  <f 
0 

where  the  terms  were  defined  in  the  previous  section  on  erosive  burning. 

Equations  110,  111  and  112  were  then  perturbed.  The  Arrhenius  law 
was  also  perturbed  to  relate  transient  mass  flux  to  transient  surface 
temperature.  Thus 


m' 

RT 
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Combining  the  perturbed  versions  of  equation  110,  111  and  112  with 
equations  108,  110,  111  and  112  results  in  the  following  relation  be- 
tween pressure  and  velocity  coupling: 


r " 0:WB  A(Vs-l)  + ABC]  - ( 1/3+0. 9B)AB/Rp)  114 


where 


_ _ _ _ 1 

B = [T^/(T^-T^)][(T^/T^)C(0.5  + B)  + ~ - 0.98]  115 

C = 1 + 8 + C(0.5  + 6)RT^/E5  116 

c = V(Vf5)  + C^yo  117 
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and  Q is  the  energy  required  to  heat  the  propellant  to  the  surface 
temperature  and  to  transform  it  into  gases.  The  parameter  A is  defined 
the  same  as  for  the  Denison  and  Baum  model,  i.e.,  it  is  related  to  the 
surface  pyrolysis. 

Thus  the  linearized  response  of  the  burning  rate  of  propellants 
to  small  pressure  and  velocity  oscillations  around  a mean  turbulent 
boundary  layer  flow  was  obtained.  The  results  of  calculations  by 
Lengelld"  showed  that  the  pressure  coupled  response  is  rather  strongly 
amplified  when  the  erosive  effect  becomes  more  and  more  pronounced. 

The  velocity  coupled  response  follows  the  same  trend. 

Lengelle's  nonsteady  state  combustion  model  represents  the  only 
attempt  to  model  theoretically  the  phenomenon  of  velocity  coupling.  An 
important  implication  derived  from  Lengelle's  work  is  that  erosive  con- 
ditions must  exist  before  velocity  coupling  can  manifest  itself.  There 
exists  some  disagreement  in  the  solid  propellant  propulsion  community 
regarding  whether  or  not  erosive  conditions  must  be  present  in  order  for 
velocity  coupling  to  occur.  This  argument  is  based  on  claimed  evidence 
of  velocity  coupling  in  T-burners  operating  in  the  velocity  coupled 
mode  in  the  absence  of  erosive  conditions.  The  T-burner  results  are, 
however,  subject  to  interpretation. 

Lengelle's  technique  offers  the  most  logical  and  physically  plaus- 
ible mechanism  for  velocity  coupling.  The  nonsteady  model  formulated 
by  Lengelle  is,  however,  based  on  a burning  rate  model  which  is  primarily 
empirical  in  nature  and  thus  leads  to  a nonsteady  model  which  is  empirical 
in  nature  and  has  little  predictive  capability. 
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THE  PETITE  ENSEMBLE  MODEL 


Introduction 

Conventional  analyses  of  nonmetal ized  composite  solid  propellant 
combustion  are  restricted  to  either  one  flame  or  one  spatial  dimension. 
Thus,  the  details  of  the  combustion  process  that  are  related  to  the 
physical  structure  of  composite  propellants  are  smeared  or  lost  in  the 
analyses  and  the  calculated  results  are  qualitative  at  best.  Statis- 
tical analysis  of  composite  propellants  combustion,  on  the  other  hand, 
offers  the  possibility  of  quantitative  results.  The  history  of  this 
approach  began  when  C.  E.  Hermance  (11)  combined  statistical  concepts 
with  a detailed  quasi-one  dimensional  physiochemical  model  in  an 
analysis  of  the  steady  state  combustion  of  composite  solid  propellants. 

In  this  work,  statistical  concepts  were  employed  to  relate  the  char- 
acteristic dimension  in  the  physiochemical  model  to  the  oxidizer  particle 
size. 

In  1970  Beckstead,  Derr,  and  Price  (18)  continued  the  same  "statis- 
tical plus  physiochemical  model"  approach  adopting  Hermance' s statistics 
but  advancing  an  improved  physiochemical  model. 

Unfortunately,  the  Hermance  and  Beckstead,  Derr  and  Price  approaches 
are  questionable  on  two  accounts,  the  physical  validity  of  replacing 
the  behavior  of  an  ensemble  of  different  flames  by  that  of  a single 
characteristic  flame  and  the  statistical  process  employed  to  select  the 
characteristic  dimension. 
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In  1970  Miller,  Hartman,  and  Myers  (20)  overcame  the  aforementioned 
problems  by  assuming  that  the  burning  surface  was  an  ensemble  of  quasi- 
one-dimensional  flamelets.  Miller,  Hartman  and  Myers  approached  the 
statistical  aspect  by  showing  that^  for  a planar  burning  surface,  the 
population  of  panicles  at  the  burning  surface  was  the  same  as  in  the 
propellant.  Thus,  it  was  assumed  that  the  area  associated  with  each 
particle  size  at  the  surface  was  related  to  the  distribution  of  par- 
ticle sizes  in  the  propellant.  The  Summerfield  GDF  model  was  imbedded 
in  this  statistical  frame  work  resulting  in  a model  capable  of  cor- 
relating the  effect  of  oxidizer  particle  size  distribution.  The  ana- 
lytical results  of  this  model  agree  well  with  experimental  data.  How- 
ever, the  model  is  based  on  an  empirical  burning  rate  model.  Thus,  it 
is  itself  empirical . 

The  statistical  formalism  presented  in  this  paper  was  originally 
conceived  by  Glick  (51).  The  combining  of  Click's  statistical  frame 
work  with  a physiochemical  model  similar  to  that  of  the  BDP  model  was 
first  reported  by  Glick  and  Condon  (10)  in  1974.  In  that  model  the 
macro  flame  structure  above  the  burning  surface  of  a composite  propellant 
was  treated  as  a collection  of  different,  non-interacting,  quasi-one- 
dimensional  flames  and  the  burning  surface  was  treated  as  a fuel  plane 
dotted  with  concave  and/or  convex  oxidizer  surfaces.  The  model  has  been 
termed  the  Petite  Ensemble  Model  (PEM).  A complete  description  of  the 
model  follows. 


The  Physical  Model 

The  burning  surface  of  a composite  propellant  can  be  visualized  as 
a random  arrangement  of  polydisperse  oxidizer  particles/fuel  surface 
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pairs.  A planar  section  along  with  an  elevation  depicting  the  burning 
surface  of  a polydisperse  propellant  appears  in  Figure  17.  The  poly- 
disperse  propellant  burning  surface  is  comprised  of  individual  oxidizer 
crystals  protruding  above  or  recessed  below  the  planar  burning  surface 
depending  on  the  combustion  pressure.  There  exists  a distribution  of 
particle  sizes  at  the  burning  surface  due  to  the  random  mixing  of  the 
various  oxidizer  particles  during  the  propellant  mixing  process. 

Each  oxidizer  particle  at  the  burning  surface  has  associated  with 
it  some  portion  of  the  available  fuel  at  the  burning  surface.  The  dif- 
fusion flame  above  the  surface  consumes  the  decomposed  gaseous  AP  and 
fuel  binder.  As  a result,  each  oxidizer/fuel  surface  pair  will  produce 
a unit  flame. 

If  it  is  assumed  that  all  oxidizer/fuel  surface  pairs  or  unit 
flames  burn  independently  of  each  other,  the  propellant  surface  can  be 
rearranged  into  imaginary  families  of  monodisperse  propellants  (pseudo- 
propellants) containing  one  oxidizer  type.  This  rearrangement  into 
subareas  is  depicted  in  Figure  18. 

Since  the  burn^'ng  rate  of  the  propellant  is  a necessary  parameter 
for  both  the  steady  and  unsteady  models,  the  burning  rate  of  a polydis- 
perse propellant  is  calculated  from  the  sum  of  the  monodisperse  subarea 
mass  fluxes  as  determined  from  an  appropriate  unit  flame  combustion 
model.  The  unit  flame  combustion  model  selected  should  be  capable  of 
predicting  the  effects  of  oxidizer  size,  oxidizer  mass  fraction,  oxidizer 
volume  fraction  and  psuedopropel lant  density  on  the  combustion  processes 
of  a unit  flame.  These  desired  characteristics  are  satisfied  by  a highly 
modified  version  of  the  original  Beckstead  Derr  Price  (BDP) 
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Figure  17.  The  Packing  Structure  and  Burning  Surface 
of  a Polydisperse  Propellant 


93 


TYPE  A,  SIZE  I 


Konodi sperse  Pseudopropellants 


iltiple  flame  combustion  model  (18). 


In  order  to  treat  the  distribution  of  oxidizer  particle  sizes  on 
the  surface  correctly,  the  statistical  formalism  for  describing  the 
heterogeneous  surface  structure  of  composite  propellants  derived  by 
Click  (51)  will  be  used  in  deriving  the  burning  rate  for  the  steady 


state  model. 


The  Theoretical  Steady  State  Model 


The  theoretical  model  for  the  steady  state  combustion  of  a poly- 
disperse  composite  propellant  is  concerned  primarily  with  specifying  a 
relationship  for  the  burning  rate  of  the  propellant  as  a function  of  the 
requisite  physical  and  chemical  parameters  outlined  in  the  physical 


model . 


Derivation  of  the  Statistical  Formalism 


The  derivation  is  initiated  by  applying  the  conservation  of  mass 
equation  at  the  burning  propellant  surface.  Thus 


dS  - 

s--  -^^p 


where  m^  is  the  average  mass  flux  per  unit  area  from  the  burning  surface 
m"  is  the  mass  flux  per  unit  surface  area,  S denotes  surface  area,  r 
is  the  average  burning  rate,  is  the  propellant  density,  and  the  sub- 
scripts b and  p designate  burning  and  planar  areas  respectively. 

The  oxidizer  particle/fuel  surface  pairs  with  common  particle 
diameters  and  species  can  be  rearranged  into  pseudopropel lants  so  that 
those  with  common  particle  diameters  and  species  are  neighbors.  Each 
oxidizer  particle/fuel  surface  pair  is  then  assumed  to  burn  independently 
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of  dll  Other  oxidizer  particle/fuel  surface  pairs.  By  rearranging  the 
burning  surface  area  into  Q monodisperse  pseudopropellants  with  s 
oxidizer  types,  equation  118  becomes 


^ f=l  k*l 


AS 


b,d,k 


119 


where  is  the  average  total  mass  flux  per  unit  area,  AS|j  is  the 
portion  of  the  burning  surface  occupied  by  the  pseudopropellant  and 
^ is  the  mass  flux  from  pseudopropellants  possessing  oxidizer  par- 
ticles with  diameters  between  and  + dD^  and  oxidizer  specie  k. 

Application  of  the  mean  value  theorem  for  integrals  to  the  right 
side  of  equation  119  yields 


0 s 
A ( I 

1=1  k=l 


, AS.  . . )/S 
d,k  b,d,k'  p 
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The  mass  flux  term  m^  is  based  on  the  burning  surface  area.  The  mass 
flux  based  on  the  planar  area,  mp  ^ can  be  written 


m"  . AS,  ^ K = n'"  ^ L.  AS„  . . 

Q,K  b,d,k  p,d,k  p,d,k 
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where  AS_  j i,  is  the  projection  of  the  monodisperse  pseudopropellant 

P »Q  » K 

burning  subarea  AS^j  (j  the  propellant  planar  surface. 

If  an  j 1,  is  the  number  of  oxidizer  particles  on  Sp,  the  planar 

p )U  9 K r 

surface,  having  diameters  between  Op  and  Dp  + dDp  and  oxidizer  specie 
k per  unit  planar  area,  then 


^ k/S„  = AS^  . . AN„  . ^ 

pyUyK  p P^Clyk  pfCXfK 
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where  ASp  ^ represents  the  average  value. 

A distribution  function,  Fp  ^ can  be  defined  such  that 


dN„  . . = NF„  . . dD^ 

p,d,K  p,d,K  0 
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where  N is  the  total  number  of  particles  on  the  burning  surface. 
Combining  equations  120  through  123  yields 


■"t  = y,  -"j.d.k  ^Sp.d.k  ^.d.k^^o^ 
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Evaluating  the  sum  over  all  f as  Q -*  “>  yields 


m"  = N I 
k=l 


'"p.d.k  ‘^^p.d.k  ^p,d,k  *^^0 
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Equation  125  is  the  statistical  formalism  that  enables  the  mean  burning 
rate  of  a propellant  with  several  different  polydisperse  oxidizers  to 
be  computed  from  the  mean  burning  rates  of  a sequence  of  monodisperse 
pseudopropellants  assuming  the  combustion  of  each  unit  flame  to  burn 
independently  of  adjacent  ones.  The  above  development  permits  consider- 
ation of  the  real  propellant  having  an  oxidizer  size  distribution  as  well 
as  having  several  different  oxidizers.  It  represents  an  improvement  over 
the  conventional  requirement  that  a single  particle  be  selected  to  re- 
present a polydispersion  of  oxidizer  particle  sizes. 

Equation  125  is  not,  however,  in  a useful  form  since  the  average 
planar  surface  area  of  a pseudopropellant  subarea,  AS  . . , and  the 

P » K 


distribution  function,  F . . , have  not  been  related  to  propellant  for- 

P tO  » K 


mulation  variables.  Moreover,  relations  for  the  pseudopropellant 
properties:  oxidizer  mass  fraction,  oxidizer  volume  fraction  and 
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pseudopropellant  density  must  be  formulated  since  they  along  with 
oxidizer  particle  size,  determine  the  mass  flux,  mj*|^,  emanating  from 
a pseudopropellant  subarea-. 

In  order  to  derive  the  above  relations  the  mean  statistical  char- 
acteristics of  the  fuel  surface  in  each  monodisperse  subarea  must  be 
determined.  This  requires  that  a statistical  determination  be  made  of 
the  arrangement  of  the  smaller  particles  inside  the  packing  of  the  larger 
particles.  Such  information  is  unavailable.  However,  an  examination 
of  the  characteristics  of  regular  geometric  packings  (52)  suggests  that 
the  mean  volume  of  fuel  associated  with  a particle  should  be  roughly 
proportional  to  its  surface  or 

AV-  ^ aoj  126 

f ,d  0 

where  AV^:  ^ is  the  average  volume  of  fuel  associated  with  an  oxidizer 
particle  of  diameter  between  Dq  and  +,dDQ.  To  permit  some  variation 
of  fuel  volume  for  a sphere,  it  will  be  assumed  that 

• CdJ  127 

where  C is  the  proportionality  constant  and  n is  an  empirical  constant. 
Once  a value  is  assumed  for  n,  C can  be  determined  in  terms  of  propel- 
lant formulation  variables  and  pseudopropellant  properties  by  considering 
conservation  of  mass. 

Since  particle  diameter  alone  is  important  in  a packing,  the  average 
volume  of  fuel  associated  with  an  oxidizer  particle  is  independent  of 
the  composition.  Thus 
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Finally,  the  volume  fraction  of  fuel  associated  with  particles  having 
diameters  between  and  + dDQ  and  oxidizer  type  k Is 

■"f.d  • ‘^f.d  ■‘"v.d 

where  dN^  ^ is  the  number  of  oxidizer  particles  per  unit  volume  with 
diameters  between  Dq  and  Dq  + dD^. 

Similarly,  the  volume  f'^action  of  oxidizer  particles  with  diameters 
between  and  + dDg  is 

= ‘'o.d  “"v.d 

and  the  mean  volume  of  oxidizer  particles  with  diameters  between  and 
0„  * d0„  Is 

‘Vo.d  ■ "“d/®  '31 

Combining  equations  130  and  131,  the  number  of  oxidizer  particles  per 
unit  volume  with  diameters  between  and  + dD^  can  be  written 

dNy^^  = 6dc/iTD^  132 

If  the  fraction  of  particles  with  diameters  between  and  + dO^  and 
oxidizer  type  k in  the  propellant  is  |^,  the  number  of  these  particles 
per  unit  planar  surface  area  is 


‘"'p.d.k  ■ ■'d.kO^d  ■ ("'’“o'  ‘"'d.k 
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The  total  volume  fraction  of  fuel  Is  (l-C)  where  c,  Is  the  total 
volume  fraction  of  oxidizer.  The  total  volume  fraction  of  fuel  Is  also 
given  by  the  Integral  of  equation  129  over  all  oxidizer  particle 
diameters.  Thus 


Employing  equations  127,  129,  131,  132,  and  134  and  rearranging  yields 
an  expression  for  the  proportionality  constant  C. 


C « {iT(l-q/6}/  [ 0"‘^d!:.  135 

With  both  particle  diameter  and  mean  volume  of  fuel  associated  with 
that  particle  diameter  known,  the  volume  and  mass  fractions  and  the 
density  of  a pseudopropellant  formed  from  an  oxidizer  particle/fuel 
surface  pair  can  now  be  computed.  The  pseudopropellant  volume  fraction 
of  oxidizer  Is 


^d,k  “ '^^o,d,k'^^*^o,d,k  * ^''f,d,k^ 


'^d,k  * ^ 6C0"‘V)“’ 


The  mass  fraction  of  oxidizer  Is 


“d,k  ” ^^o.d,k  f'o,k/<'^''o,d,k  Po,k  ^ *''f,d,k 
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Thus 


{1  + 


and  the  pseudopropellant  density  Is 


^d,k  * ^p.k  ^d,k^°‘d,k 
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* * * 

The  above  pseudopropellant  properties  (a^  1^,  r,^  and 
along  with  the  particle  diameter,  are  used  in  the  unit  flame  com- 
bustion model  to  determine  the  mass  flux  emanating  from  the 

pseudopropellant  subareas. 

In  order  to  relate  the  average  planar  surface  area  of  a pseudo- 
propellant subarea,  1^,  to  propellant  properties,  it  is  first 

necessary  to  consider  the  statistical  characteristics  of  the  intersection 
of  a plane  with  a randomly  packed  polydispersion  of  spheres.  It  can  be 
shown  that  the  mean  statistical  characteristics  of  the  intersection  of 
a plane  with  a randomly  packed  polydispersion  of  spheres  is  the  same 
as  that  for  a randomly  packed  monodispersion  of  spheres.  Namely,  if 
one  considers  all  particles  being  intersected  that  have  diameters  be- 
tween Dq  and  Dq  + dD^,  the  Blum  and  Wilhelm  (21)  mean  intersection 
diameter  applies.  Thus 


D “ Dp 


141 


where  D is  that  mean  intersection  diameter. 

The  average  oxidizer  particle  intersection  diameter  can  now  be 
determined.  Thus, 
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Since  the  volume  fraction  of  oxidizer  in  the  pseudopropellant  subarea  is 


* ♦ 

'^d.k  * ^o.d.k''^p.d.k  “ ‘''^o.d.k^^^p.d.k 


equation  142  becomes 


•'^^p.d.k  * 

Equation  144  provides  the  necessary  relationship  for  determining 
the  average  planar  surface  area  of  a pseudopropel lart  subarea  with 
oxidizer  particles  having  diameters  between  D^j  and  0^  + dD^  and  for 
oxidizer  species  k. 

During  the  derivation  of  the  expressions  for  the  pseudopropellant 
properties  a new  unknown,  d^^j  |^,  the  volume  fraction  of  oxidizer  having 
particle  sizes  between  and  + dD^,  was  introduced.  It  now  becomes 
necessary  to  relate  d;^  ^ and  the  only  other  remaining  unknowns 
(the  distribution  function)  and  (the  propellant  density)  to  the  in- 
dependent propellant  formulation  variables  the  mass  fraction  of 

oxidizer  of  specie  k and  mode  j.  El,  ^ the  distribution  function  of 
oxidizer  particle  size  of  oxidizer  specie  k and  mode  j.  the  number  of 
modes  for  specie  k and  s the  number  of  oxidizer  species. 

The  volume  fraction  of  oxidizer  specie  k having  diameters  between 
Dq  and  Dp  + dD^  is  defined 


■‘‘a.k  ■ 
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Since 


dV  . , = dm  j l/p_  L 
Ofdfk  o»d»k  o»k 


Combining  equations  145,  146,  and  147  yields 


The  distribution  function  can  be  defined 


fk.j  ■ ’'"o,().k,J''<'"o.k.J  “o*  ’ 

where  the  element  of  mass  of  oxidizer  of  specie  k and  mode  j and  of 

sizes  between  D and  D + dD  is 
0 0 0 

dm  j I ~ y ^ Lf  i 

o,d,k  o,d,k,j 

j ' 

Combining  equations  148,  149,  and  150  yields 


o^d.k  * <'’t^‘’o.k>  T,  “k.j  ^k.j 


Since  the  volume  of  oxidizer  with  diameters  between  Dq  and  Dq  + dD^  is 
given  by  the  summation  of  volumes  associated  with  each  specie,  then 


■“d' J,  V^.k’  ^k 


J 


where 


\ 

^ ■ J,  'kj  “k-j 


The  volume  fraction  of  oxidizer  is 


^0  = 


but 


s \ 


Jl  j = l 


Combining  equations  147,  154,  and  155  yields 

^0  ■ J,  / ‘^,ky"t' 


Since 


c.  . - m , ./m. 
k,j  o,k,j  t 


and 


P,  , " ni.,  ,,  ;/V, 
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155 


156 


157 

\ 
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o.k  * o,k,j^  o,k,j 
Combining  equations  155,  157,  and  158,  the  volume  fraction  of  oxidizer 
can  be  written 


Ml 
s k 


'•0  ■ "t  j,  T,  '''k,j'‘’o,k> 
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The  propellant  density  can  be  written  as  follows: 


Pt  = m^/V^ 


160 


where 


s Mk 


't  = I,  ''o.kj  ^ ''f 
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Since 


s \ 


'’t  ■ ‘k?i  ?i  * (Vf/iijHmf/VfPf)) 


-1 


162 


Combining  equation  162  with  equations  157  and  158  and  noting  that 


yields 


(1  -a)  = m^/m^ 


s "k 
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Pt  = “k.j/Po.k,o  " ^ " “V 


-1 
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but 


s «k 


“ " ^ ^ “k  j 

k=l  j=l 
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Therefore 


s \ 


^t  " ^^''Pf  “k.j^ 


-1 


166 


Employing  equations  123,  133,  144,  and  152,  equation  125  can  be  re- 
written 
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■*  ■ "o'.k  fn  ^ 


Therefore 


^o.k  ,p  ^"’J.d.k^^d.k^  ’'k  ‘*°o 


Equation  168  expresses  the  polydisperse  propellant  burning  rate  in  terms 
of  propellant  formulation  variables  and  the  mean  state  subarea  mass 


fluxes. 


In  order  to  implement  the  computation  Indicated  by  equation  168, 
it  is  necessary  to  specify  a mathematical  relationship  for  the  distri- 
bution function,  representing  the  distribution  of  oxidizer  particle 


sizes. 


A log  normal  distribution  was  selected  as  representative  of  the 
particle  size  distribution.  The  assumption  of  a log  normal  distribu- 
tion permits  each  particle  size  distribution  to  be  characterized  by 
two  parameters:  the  weight  median  diameter,  D,  of  particles  in  the 
distribution  mode  and  a parameter,  o,  which  reflects  the  width  of  the 
distribution  of  particle  sizes  about  the  weight  median  diameter  of 
particles  in  the  mode. 

The  log  noriikil  distribution  may  be  represented  by  the  following 
relationship: 


y » 1/(0  2ti)  exp  [-1/2  (~-)  3 


where 


X * ^ 0 

170 

m = In  5 
(n 

171 

* 

a = In  a 

172 

» • ”o 

173 

, the  weight 

percent  of  particles  less  than  a given  diameter 

1 1 / ( 2vtna ) ) 

[ ^ exp  [-  ^ ( (£KD-fnD)/£no)^]dD  174 

J 

0 

This  two  parameter  method  of  representing  the  particle  size  distribution 
of  a mode  is  in  excellent  agreement  with  experimental  particle  size 
data  for  most  oxidizer  grinds. 

Figure  19  depicts  a comparison  of  the  log  normal  distribution  to 
particle  size  measurements  made  by  Miller  (6,7)  for  several  different 
particle  size  distribution  modes.  The  standard  deviations  of  the  experi- 
mental data  from  the  log  normal  distributions  which  result  when  the 
weight  percent  of  particles  less  than  a diameter  is  plotted  as  a function 
of  particle  diameter  was  never  more  than  3.1%  and  on  the  average  was 
only  1.7%.  Thus,  it  is  apparent  that  the  log  normal  approximation  is 
in  good  agreement  with  the  actual  distribution. 

The  formalism  has  now  been  written  for  computing  the  burning  rate 
of  a polydisperse  composite  propellant  given  the  average  mass  flux  from 
each  pseudopropellant. 
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WEIGHT  PERCENT  LESS  THAN  DIAMETER 


Derivation  of  the  Physiochemical  Model 


The  BDP  multiple  flame  model  was  used  as  a foundation  upon  which 
to  build  a more  realistic  physiochemical  model.  A complete  description 
of  the  BDP  model  can  be  found  in  the  literature  review.  The  major 
similarity  between  the  BOP  and  PEM  models  is  the  assumed  flame  structure. 
However,  the  assumed  surface  geometry  and  statistical  approach  differ 
extensively.  The  diffusion  flame  standoff  distance  and  the  AP  flame 
rate  constant  calculations  also  differ  from  the  BDP  approach.  A 
complete  derivation  of  the  physiochemical  model  will  follow. 

The  Mean  Deflagrating  State.  A critical  problem  in  statistical 
combustion  modeling  is  to  define  an  oxidizer  particle's  mean  deflagrating 
state.  The  method  employed  is  to  calculate  the  burning  rate  of  the  oxi- 
dizer particle  as  a function  of  time  during  the  entire  lifetime  of  the 
I particle  from  its  birth  at  the  propellant  burning  surface  until  the 

particle  is  fully  consumed.  The  regression  rate  is  calculated  at 
specific  increments  during  the  regression  of  the  particle  in  order  that 
a rate  versus  time  curve  can  be  generated.  The  mean  mass  flux  can  then 
be  determined  by  the  numerical  integration  of  the  curve. 

- 1 f" 

m = - mdt  175 

■^0 

where  t is  the  lifetime  of  the  particle  at  the  burning  surface.  This 
process  is  illustrated  in  Figure  20.  The  calculation  of  the  pseudopro- 
pellant mass  flux  at  the  indicated  increments  is  based  on  the  assumption 
that  the  mixture  ratio  of  oxidizer  to  fuel  remains  constant  during  the 
combustion  of  the  oxidizer  particle.  A complete  description  of  the 

j 


Is 
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equations  used  to  calculate  the  mass  flux  as  well  as  the  solution 
methodoloay  will  follow. 

The  Flame  Structure.  The  flame  structure  is  depicted  in  Figure  8. 

A complete  description  of  the  structure  was  given  in  the  BOP  section  of 
the  literature  review  and  will  not  be  discussed  here. 

The  Conservation  of  Hass.  For  a burning  propellant,  the  equation 
of  continuity  may  be  written 

= rpp  = 1/6 

where  m represents  the  mass  flux,  S the  surface  area,  and  the  subscripts 
f,  ox  and  T stand  for  fuel,  oxidizer  and  total  surface  area  respectively. 
This  equation  can  be  simplified  further  since  it  is  assumed  that  the 
ratio  of  oxidizer  consumption  to  fuel  consumption  remains  constant 
during  the  consumption  of  the  oxidizer  particle. 


177 


where  n*  is  the  mass  fraction  of  oxidizer  in  the  pseudopropellant.  Thus 


or 


(S^/S.j.)m^/(l-ot  ) 


178 


179 
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Thus,  the  burning  rate  c,in  be  calculated  from  either  equation  178  or 
equation  179.  Since  it  is  assumed  that  the  oxidizer  regression  controls 
the  overall  burning  rate,  equation  178  is  chosen  as  the  continuity 
equation. 

The  Pseudopropellant  Surface  Geometry.  To  evaluate  the  area  ratio 
in  equation  178,  the  geometric  relationship  between  the  oxidizer  and  the 
binder  must  be  established.  The  oxidizer  surface  area  at  any  point 
during  its  consumption  is  calculated  by  assuming  that  the  regressing 
oxidizer  surface  is  spherical  in  shape  as  depicted  in  Figure  21,  where 
R is  the  radius  of  curvature  of  the  surface.  It  is  assumed  that  the 
fuel  binder  regresses  continually,  but  the  oxidizer  is  assumed  to 
undergo  an  ignition  delay  upon  reaching  the  burning  propellant  surface. 
The  ignition  delay  time  is  assumed  to  have  the  form  derived  by  Shannon 
(25).  Thus 


t.  = C. 
ign  ign  o ^ 
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where  C^gp,  U and  U are  constants.  is  the  oxidizer  particle  diameter 
and  P is  the  pressure.  Therefore,  the  distance  d^^^  in  Figure  21  is 


d = 1/p 
ox  ox 


'ign. 

m^  dt 
ox 
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the  distance  d^  is 


d^  = 1/P^ 


m^  dt 
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and  h in  Figure  21  is 


" • ■’ox  ■ “f 

The  radius  of  curvature,  R,  can  be  determined  by  triangle  I in  Figure 
21 . Thus 


=«  (R  - h)^  + (072)^  18^ 

where  O'  is  the  intersection  diameter  of  the  oxidizer  particle  with  the 
fuel  plane  at  any  instant.  Therefore, 

R - h/2  + (D72)^/2h  135 

The  surface  area  of  the  spherical  segment  representing  the  burning 
AP  surface  is 


Combining  equations  185  and  186  yields 


186 


Sqx  = ^ (072)^]  187 

The  intersection  diameter.  O',  can  be  determined  from  geometric  consid- 
erations 


D' 


188 
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The  total  burniny  surface  area  for  the  pseudopropel lant  is  yiven  by 
the  sum  of  the  planar  fuel  area  with  the  oxidizer  surface  area.  Since 
the  mixture  ratio  of  oxidizer  to  fuel  is  assumed  to  be  constant  through- 
out the  combustion  of  the  oxidizer  particle,  the  planar  surface  area  of 
the  fuel  binder  associated  with  the  pseudopropellant  can  easily  be 
determined.  From  Figure  22,  it  is  clear  that 


P 

ox 
+ S 


P 

ox 
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where  is  the  volume  fraction  of  oxidizer  in  the  pseiidopropel lant 

and  is  the  projection  of  the  oxidizer  surface  on  the  fuel  plane, 
ox 

The  planar  oxidizer  area  is 


190 


and 


SP  - nb^  191 

ox  f 

Combining  equations  189,  190,  and  191  yields 
b » 07(2x'^) 

The  planar  fuel  surface  area  can  now  be  determined.  Thus 

sf  * WD72)^[1/;*  - 1]  193 
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and  the  total  surface  area  Is 


Sf  * Sf  ^ Sjx  * ♦ (0V2)"J  ‘ {L>'/2)^[1/-/  - U 


' 194 


The  ratio  of  oxidizer  surface  area  to  total  pseudopropel lant  surface 
area  is 

(S^jx/Sj)  * [h^  » (P72)^]/{[h^  ^ (DV2)^]  ^ {D72)^[l/:/  - 1]}  1 

Surface  Pyrolysis.  The  mass  flux  of  the  oxidizer,  is  cal- 
culated by  assuming  that  global  kinetics  adequately  describes  the 
decomposition  of  the  oxidizer  and  the  fuel.  Thus,  utilizing  an 
Arrhenius  expression,  the  mass  flux  of  oxidizer  can  be  expressed 


m =>  A i'xp  (-E  /RT  ) 

ox  OK  ^ OX  S 


Similarly  for  the  fuel  pyrolysis  the  mass  flux  can  be  expressed 


If  • cxy  (-tj:/RT^) 


where  A is  the  pre-exponent ial  constant,  E is  the  effective  activation 
energy,  R is  the  gas  constant  and  T^  is  the  surface  temperature.  It 
should  be  noted  that  equation  31  is  used  only  in  the  determination  of 
the  surface  geometry  and  is  not  used  directly  for  the  burning  rate 
calculation.  The  surface  temperature  used  in  equations  J1  and  32  is 
assumed  to  he  the  average  surface  temperature  for  the  entire  pseudo- 
propellant  surface. 
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The  Conservation  of  Energy  at  the  Surface.  It  is  apparent  from 
equation  32  that  the  surface  temperature  plays  an  important  role  in 
determining  the  pseudopropellant  mass  flux.  The  surface  temperature 
is  determined  by  an  energy  balance  at  the  burning  surface.  The  energy 
balance  is  depicted  in  Figure  23.  The  energy  required  to  raise  the  pro 
pellant  from  the  initial  propellant  temperature,  T^,  to  the  surface 
temperature,  T^,  plus  the  energy  required  to  vaporize  the  oxidizer, 
plus  the  energy  required  to  vaporize  the  fuel  binder  is  equal  to  the 
energy  arriving  at  the  surface  from  the  gas  phase  flame  zone.  Assuming 
constant  specific  heat,  the  energy  balance  at  the  surface  can  be 
written 

mfCpCTs-To)  + am-^Qj_  + {l-a)m.j.Q^  = ENERGY  ARRIVING  FROM  FLAME  ZONE  196 

where  is  the  heat  of  pyrolysis  of  the  fuel  binder  and  is  the  heat 
of  gasification  of  the  oxidizer.  The  heat  of  gasification  of  the 
oxidizer  includes  the  endothermic  decomposition  step  plus  an  exothermic 
decomposition  step  plus  the  exothermic  reactions  which  occur  on  the  de- 
composing surface. 

In  order  to  evaluate  the  amount  of  energy  received  by  the  surface 
from  the  flame  zone,  it  is  assumed  that  the  two  dimensional  flame  shape 
can  be  represented  by  quasi  one  dimensional  flames  for  the  purpose  of  a 
heat  transfer  analysis.  This  is  illustrated  in  Figure  24. 

If  it  is  assumed  that  the  gas  phase  energy  release  is  concentrated 
only  at  the  flame  location,  the  1-D  energy  equation  for  the  region  be- 
tween the  final  flame  and  the  AP  flame  becomes 
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Figure  23.  The  Pseudopropellant  Surface  Energy  Balance 
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The  boundary  conditions  are 

0 X = X^p  T = T(X^p)  198 

0 X = Xpp  T = Tpp  199 

where  X^p  and  Xpp  are  defined  in  Figure  24,  Tpp  is  the  adiabatic  final 
flame  temperature  and  T(Xy^p)  is  the  temperature  at  the  location  of  the 
AP  flame.  With  the  above  boundary  conditions,  the  solution  to  equation 
197  becomes 

T = Tpp  + {Tpp-T(X^p))[cxp(m^^CpX/k) 

- cxp(mg^CpXpp/k)]/[cxp(m^^Cp  Xfp/k) 

An  energy  balance  at  the  final  flame  yields 


m Qrn 
ox^FF 


201 


where  Qpp  is  the  heat  release  in  the  final  flame  per  unit  of  mass  flux 
of  oxidizer. 

Evaluating  the  derivative  of  equation  200  with  respect  to  X and 
letting 

E » mCpX/k  202 
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yields 


%x^FF  ^ %x^p^^FF  ■ ' ‘^’‘p(^*AP^^ 

Equation  203  can  be  rearranged  to  yield  an  expression  for  the  AP  flame 
temperature. 

T(Xap)  = Tpp  - Qppd  - cxp(Eap  - 

The  energy  equatioh  (197)  can  also  be  solved  for  the  region  between  the 
AP  flame  and  the  burning  surface  with  the  following  boundary  conditions 


P X = Xap 

T • T(*,p) 

205 

0 X = 0^ 

^ 

206 

The  solution  is 

T = Ts  + {T(Xap)  - Tg)(exp(rJ  - 1 )/(e.xp(i;Ap)  - D 207 

The  energy  received  by  the  surface  from  the  final  flame  and  the  AP  flame 
may  now  be  determined 


m 


ox 


AP.FF 

surf 
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Evaluating  the  derivative  of  equation  207  and  combining  the  result 
with  equation  T fields 


The  surface  temperature  in  equation  209  can  be  eliminated  by  performing 
an  energy  balance  at  the  AP  flame 


9X  . * "OX^AP  " 3X  |y  y 


where  Q^p  is  the  heat  released  in  the  AP  flame  per  unit  of  mass  flux  of 
oxidizer.  The  derivative  on  the  left  hand  side  is  evaluated  with 
equation  200  while  the  derivative  on  the  right  hand  side  is  evaluated 
with  equation  207.  Equation  210  becomes 


Combining  equation  211  with  equation  209  the  heat  flux  to  the  surface 
from  the  final  flame  and  the  AP  flame  becomes 


• Ap  Fp  • . • 

V^surf  " '"ox^AP'^’^P^'^AP^  %x^FF‘^’^P^'^FF^ 


In  a similar  fashion,  the  heat  flux  to  the  surface  from  the  primary 
flame  is 


"T^surf  " "’t^PF‘^’‘P^‘^FF^ 
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The  surface  energy  balance  can  now  be  written 

m^Cp(T^-T^)  = - - (l-a)m^O^  + B|-Qppm^Pxp( -Cpp) 

+ (l-8p)mja[Q^pCxp(-C^p)Qppexp(-Cpp)]  2H 

where  the  parameter,  appearing  in  equation  d'i  represents  the 
fraction  of  the  reactants  that  react  in  the  primary  flame.  Equation 
214  can  be  rearranged  so  that  the  average  pseudopropellant  surface 
temperature  becomes 


Ts  = - aQ^/Cp  - (l-a)Qp/Cp  + (1 -Bp)a/Cp[Q^pCxp( -^^p) 


^ Qpp^^p ( ”^PP ) ^ ^pQpp^^p ) ~^’pp 
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The  energy  balances  for  the  gas  phase  heat  releases  (Q^p.  Qpp.  and  Qpp) 


were  presented  in  the  review  of  the  BDP  model.  Thus 


Qap  = (T,p  - TJ  ^ Q, 


AP 
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and 


'IpF  " S ■ ^o>  * *(!-»)  Qf 
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where  T is  the  adiabatic  flame  temperature  oi  the  AP  flame.  Thus,  the 
AP 


energy  equation  at  the  surface  yields  a value  for  the  average  pseudo- 
propellant surface  temperature  which  is  a necessary  parameter  for  the 
burning  rate  calculations. 

Flame  Standoff  Distances.  The  nondimensional  flame  standoff  dis- 
tances appearing  in  the  surface  energy  balance  are  defined 
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* '"ox  ^’^AP  ^ 

CpF  ■ Cp  (Xpp  + Xpp)/k 

where  X^p  and  Xpp  are  the  kinetics  distances  for  the  AP  and  primary 
flame  respectively.  The  diffusion  distances  for  the  primary  ami  final 
flame  are,  respectively,  Xp^  and  Xpp. 

The  kinetics  distances  were  evaluated  by  integration  of  the  re- 
action rate  for  a reaction  of  arbitrary  order.  Thus 


* . 

X = m/KP° 
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where  K represents  a pseudo  rate  constant  and  is  evaluated  by  an 
Arrhenius  expression.  The  primary  flame  rate  constant  may  be  written 


*^PF  ~ ^PF  ^pp/^^pp) 
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where  Tpp  is  the  adiabatic  flame  temperature.  Similarly,  the  AP  flame 
rate  constant  is 


^AP  ~ ^AP  ^ 


AP' 


"AP' 
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where  T(X/\p)  is  the  temperature  of  the  AP  flame  given  by  equation  205 
(not  necessarily  the  adiabatic  flame  temperature). 

The  diffusion  distances  are  evaluated  by  a Burke  Schumann  Bunsen 
burner  type  of  analysis  (26).  The  geometry  for  the  analysis  is  depicted 
in  Figure  25  where  the  dimensions  6,  O'  and  b were  described  previously 
in  the  section  on  surface  geometry.  Burke  and  Schumann  solved  the  mass 
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conservation  equation  for  a Bunsen  burner  type  of  flame. 
The  diffusion  equation  can  be  written 


^ r 6r  [ 6r  . 


u 6 B 
D 61 


225 


where  B is  a concentration  term.  In  the  original  solution  to  equation 
225  by  Burke  and  Schumann  the  second  derivative  of  B with  respect  to  Z 
was  assumed  negligible.  This  assumption  is  valid  for  tall  flames  where 
the  peak  flame  height  is  large  compared  to  the  burner  dimensions.  The 
PiirkL'-Schuiiut)n  solution  of  equation  225  is 


y^U±vlc__  = I J ) J (T-  f.)/[T'.(J  oxp(f^n)]  226 

2(l+v)c“  i=l  ' ^ ° ’ 101 


where  J represents  a Bessel  function,  v is  a stoichiometric  ratio,  c is 
related  to  the  burner  geometry,  represents  roots  of  the  first  order 
Bessel  function,  F.  is  the  non-dimensional  radial  distance,  n is  the 
nondimensional  axial  distance  and  the  burner  geometry  parameter,  c,  may 
be  defined 


c = D72b 
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The  solution  to  equation  226  was  modified  by  Beckstcad  et  al.  (18) 
to  include  the  second  derivative  of  B with  respect  to  the  axial  direction, 
thus  yielding  a solution  that  is  valid  for  short  flames. 
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V - 


= I ' 
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exp 


([1  + (2^'4'.)2]’^^  - l)n 
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which  is  identical  to  equation  226  with  the  exception  of  the  term  in 
the  exponent,  and  where 


229 


Equation  228  could  be  used  to  calculate  the  actual  flame  shape, 
however  it  is  not  necessary  since  the  heat  transfer  to  the  burning 
surface  from  the  flame  is  assumed  to  be  one  dimensional.  Also,  a cal- 
culation of  the  actual  flame  shape  would  be  prohibitive  in  terms  of 
computational  time.  Thus,  a parabolic  flame  shape  is  assumed  and  equa- 
tion 228  is  used  only  to  calculate  the  peak  height  of  the  diffusion 
flame.  The  value  of  the  peak  height  calculated  by  equation  228  is  then 
multiplied  by  a constant  factor  so  that  the  result  is  the  approximate 
height  of  a one  dimensional  flame  above  the  surface  in  order  to  transfer 
the  same  amount  of  heat  to  the  surface  as  the  parabolic  flame. 

The  equation  for  the  flame  height  at  the  center  of  the  flame 
becomes 


V - (l+v)c^ 
2(l+v)c 
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exp 
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If,  however,  the  flame  is  underventilated,  the  equation  for  the  flame 
height  becomes 


V - (l-v)c^ 
2{l+\))c 


([1  + (2#.)"] 
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The  diffusion  mixing  distance  then  becomes 

xj  * bn/<*>  232 

Determination  of  the  Fraction  of  Reactants  that  React  in  the  Primary 
Flame.  The  only  remaining  unknown  in  the  surface  energy  balance  is 
the  fraction  of  reactants  that  react  in  the  primary  flame.  The  Bf 
parameter  is  determined  in  the  same  manner  in  the  present  model  as  in 
the  BDP  model.  This  method  was  discussed  in  the  review  of  the  BOP  model 
and  will  not  be  repeated  here. 

Summary  of  the  Basic  Equations.  The  equations  of  the  physiochemical 
model  have  now  been  derived  and  the  main  equations  are  summarized  below. 

Total  Mass  Flux  from  Surface  of  Pseudopropellant 

178 

195 

196 

215 


^ "ox'“* 

Surface  Geometry 

(Sp^/Sf)  = [h^  + (D72)^]/{[l/  + (072)^] 

+ (D‘/2)^  [1/5*  - 1]1 

Oxidizer  Mass  Flux 

m = A exp  (-E /RT^) 

ox  ox  ' ox  s 

Surface  Temperature 

% ■ ■'o  - "VS  - 

[Q^p  cxp(-5^p)  + QppCxp(-  pp)]  ^pQpp<^*p( “^pp )/Cp 
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Solution  of  the  Basic  Equations.  The  simultaneous  set  of  equations 


which  make  up  the  model  are  solved  numerically  by  a secant  iteration  on 
the  surface  temperature.  Embedded  within  the  iteration  on  the  surface 
temperature  are  two  inner  iteration  loops.  One  is  a Newton  iteration 
on  the  solution  of  equation  228  for  the  diffusion  flame  standoff  dis- 
tance. The  other  is  an  iteration  on  the  AP  flame  standoff  distance. 
Once  the  surface  temperature  has  been  determined,  the  mass  flux  from  a 
pseudopropellant  at  a given  instant  in  time  can  be  determined  using 
equation  178. 

This  process  is  repeated  at  time  increments  over  the  lifetime  of  a 
particle  in  a particular  pseudopropellant.  The  mean  state  or  average 
mass  flux  for  the  pseudopropellant  can  then  be  evaluated  by  equation 
175 


m"  = y I in^  dt  175 

This  entire  process  is  then  repeated  at  increments  in  oxidizer 
particle  diameter  such  that  the  total  oxidizer  particle  size  distribu- 
tion is  represented  by  a series  of  pseudopropellants. 

Having  obtained  the  mean  states  for  each  of  the  pseudopropellants, 
as  well  as  the  physical  properties  of  the  pseudopropellants,  the  over 
all  composite  propellant  burning  rate  is  determined  b>  equation  168. 
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Modification  of  the  Model  to  Include 
Burning  Rate  Catalysts 

The  equations  of  the  PEM  were  modified  so  that  the  effect  of  burn- 
ing rate  catalysts  could  be  evaluated. 

There  are  two  flame  zones  in  the  PEM  model  in  which  the  catalyst 
can  manifest  itself.  One  of  those  is  the  primary  flame  which  is  the 
reaction  of  the  oxidizer  decomposition  products  with  the  fuel  binder 
products.  The  other  is  the  AP  monopropellant  flame  which  is  the  re- 
action between  the  AP  decomposition  products  to  form  oxygen  and  inert 
products.  An  Arrhenius  expression  was  assumed  for  both  the  primary 
flame  and  the  AP  flame  rate  constants.  Thus 

KpF  = ApF  exp  (Epp/RTpp)  233 

•^AP  " \p  ^^AP^'^^AP^ 

where  A is  the  Arrhenius  frequency  factor,  E is  the  activation  energy, 

R is  the  gas  constant  and  T is  the  flame  temperature. 

The  action  of  a catalyst  on  solid  propellant  combustion  was  assumed 
to  be  through  the  decrease  of  the  activation  energies  of  either  or  both 
of  the  AP  and  primary  flames. 

The  catalyst  is  assumed  not  to  affect  the  final  flame  kinetics 
(the  reaction  between  the  binder  decomposition  products  and  the  products 
of  the  AP  monopropellant  flame).  The  reason  being  that  the  reactants 
for  the  final  flame  are  preheated  to  at  least  the  AP  flame  temperature, 
thus,  the  kinetics  are  instantaneous  compared  to  the  mixing  time. 

A primary  problem  in  including  the  effects  of  catalysts  in  a 
propellant  combustion  model  is  a lack  of  knowledge  of  where  a particular 
catalyst  acts. 


131 


r ^ 


Since  the  primary  flame  is  more  predominant  at  low  pressures  and 
for  small  oxidizer  particles,  one  would  expect  that  a catalyst  which 
acts  on  the  primary  flame  kinetics  would  have  a more  pronounced  effect 
on  the  burning  rate  at  low  pressure  than  at  high  pressure.  Moreover, 
a shift  in  particle  size  distribution  to  finer  oxidizer  in  a catalyzed 
propellant  would  be  expected  to  produce  a more  pronounced  effect  on  the 
burning  rate  than  the  same  shift  in  particle  size  in  a comparable  un- 
catalyzed propellant.  The  AP  flame,  on  the  other  hand,  is  more  predom- 
inant at  high  pressure  and  for  large  oxidizer  particles.  Thus,  one 
would  expect  that  a catalyst  which  acts  on  the  AP  flame  kinetics  would 
have  a more  pronounced  effect  on  the  burning  rate  at  high  pressure  than 
at  low  pressure.  Moreover,  a shift  in  particle  size  distribution  to 
coarser  oxidizer  in  a catalyzed  propellant  would  bo  expected  to  pro- 
duce a more  pronounced  effect  on  burning  rate  than  the  same  shift  in 
particle  size  distribution  in  a comparable  uncatalyzed  propellant. 

Thus,  a small  amount  of  experimental  data  could  be  employed  to 
determine  the  action  site  of  various  burning  rate  catalysts. 

Experimental  evidence  (53)  suggests  that  there  are  two  parameters 
of  importance,  other  than  the  type  of  catalyst  employed,  in  determin- 
ing the  effect  of  catalysts  on  propellant  burning  rate.  They  are  the 
mass  fraction  of  catalyst  present  in  the  propellant  and  the  specific 
surface  area  of  the  catalyst.  Furthermore,  continued  increases  in 
either  of  the  above  two  parameters  in  a specific  propellant  produces 
diminishing  returns  in  terms  of  an  increase  in  burning  rate.  A point 
is  finally  reached  where  further  increases  in  the  mass  fraction  of 
catalyst  produce  reduced  burning  rates. 
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The  actual  effect  of  burning  rate  catalyst  on  the  chemical  reaction  5 

kinetics  is  not  completely  understood  and,  at  the  present  time,  is  be- 
yond the  scope  of  this  modeling  effort.  However,  once  the  site  of  the 
catalyst  action  is  determined,  the  effect  of  catalyst  on  the  reaction 
kinetics  can  be  handled  in  a semi-empirical  fashion. 

The  catalysts  are  accounted  for  in  the  present  combustion  model  by 
modifying  the  rate  constants  through  polynomial  expressions  for  the 
activation  energies.  The  constants  in  the  polynomial  expressions  are 
determined  by  a regression  analysis  in  which  the  parameters  of  impor- 
tance are  the  mass  fraction  and  the  specific  surface  area  of  the  catalyst. 

Since  catalysts  comprise  a distribution  of  particles,  the  addition 
of  catalyst  particles  to  a propellant  shifts  the  particle  size  distri- 
bution of  the  oxidizer  particles  such  that  there  is  a redistribution  of 
the  fuel  binder  associated  with  each  particle.  This  particular  aspect 
has  not  been  accounted  for  in  the  present  model.  However,  since  the 
amount  of  catalyst  usually  present  in  a propellant  is  quite  small  (on 
the  order  of  1%  or  less  by  weight)  this  effect  is  expected  to  be 
negligible. 

The  only  other  modification  to  the  steady  state  model  to  account 
for  catalysts  is  in  the  energy  balance  at  the  burning  surface. 

The  Erosive  Burning  Model 

The  erosive  burning  model  is  based  on  the  hypothesis  that  the 
erosive  burning  effect  is  due  to  enhancement  of  the  transport  properties 
in  the  region  of  the  gas  phase  reaction  zone  as  suggested  by  Lengelle 
(43).  The  enhancement  of  the  transport  properties  is  assumed  to  due  to 
the  presence  of  a turbulent  boundary  layer. 
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The  turbulent  boundary  layer  is  modeled  by  an  approximate  solution 
of  the  Navier-Stokes  and  continuity  equations  in  the  region  near  the 
burning  surface  for  flow  over  a flat  plate  with  injection.  In  the 
analysis  it  is  assumed  that  a turbulent  boundary  layer  develops  above 
the  propellant  surface  and  that  the  transport  coefficients  of  diffusion 
and  heat  conduction  are  modified  in  such  a way  that  the  heat  flux  to  the 
surface  and,  consequently,  the  propellant  burning  rate  are  enhanced. 

The  model  for  the  turbulent  boundary  layer  is  coupled  with  the  PEM 
to  form  the  erosive  burning  model.  The  basic  assumptions  and  the 
development  of  the  PEM  have  been  discussed  in  a prior  section.  A dis- 
cussion of  the  turbulent  boundary  layer  model  and  the  development  of  the 
equations  which  relate  transport  properties  to  the  cross  flow  velocity 
will  follow. 

The  Turbulent  Boundary  Layer.  The  turbulent  boundary  layer  model 
is  based  on  an  approximate  analysis  of  the  flow  over  a flat  plate  with 
injection.  The  basic  equations  used  in  this  analysis  are  the  Navier- 
Stokes  equations  representing  conservation  of  momentum  and  continuity. 
These  equations  are  applied  to  turbulent  flow  by  representing  the  de- 
pendent quantities  in  the  equations  as  the  sum  of  a time  averaged 
quantity  plus  an  instantaneous  quantity.  When  the  resulting  equations 
are  time  averaged  and  the  usual  boundary  layer  assumptions  are  imposed 
the  resulting  boundary  layer  equations  are  (54): 

235 

236 
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where  the  barred  quantities  represent  time  averages  and  the  primed 
quantities  represent  instantaneous  values  of  the  fluctuating  quantities. 
The  coordinate  system  and  a physical  description  of  the  model  are  de- 
picted in  Figure  26. 

For  the  present  analysis  it  is  assumed  that  the  fluid  properties 
of  viscosity  and  density  remain  constant  throughout  the  flow  field.  In 
order  to  achieve  an  engineering  solution  to  the  problem  it  is  also 
assumed  that  the  variations  of  the  dependent  variables  with  respect  to 
X are  negligible  compared  to  their  variations  with  respect  to  y.  The 
latter  assumption  allows  the  partial  differential  equations  for  con- 
servation of  momentum  and  mass  to  be  replaced  by  ordinary  differential 
equations.  It  has  been  found  (54)  that  this  assumption  combined  with 
assumptions  concerning  the  mechanism  of  turbulence  can  lead  to  solutions 
for  the  skin  friction  coefficients  which  agree  well  with  experimental 
results.  Moreover,  the  present  analysis  is  concerned  primarily  with 
the  inner  region  of  the  boundary  layer  where  gradients  in  the  y direction 
are  much  larger  than  gradients  in  the  x direction.  Thus,  the  assumptions 
should  be  valid  even  in  the  case  of  an  imposed  velocity  and  pressure 
gradient  parallel  to  the  surface,  as  is  the  case  within  a rocket  motor 
chamber. 

With  the  above  assumptions  equations  235  and  236  can  be  simplified 
to 


= 0 237 

dy 


135 


/jTjrrrrrj  rjTT  / ////////;/////////  f //,;/// • 

X 


Figure  26.  The  Physical  Model  for  the  Turbulent  Boundary  Layer 
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— du  _ d^u  du 'v ' 
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Integration  of  equation  yields 


pv  = constant  = (pv) 


where  {pv)^^  is  the  mass  flux  from  the  burning  propellant  surface.  Thus, 
equation  238  becomes 


, ..X  du  _ d u d u'v' 

If’-'w  37  ■ ^ Tp— 


The  fluctuating  terms  in  equation  240  can  be  dealt  with  by  employing 
the  Boussinesq  relation  for  the  eddy  viscosity. 


- PU'V  = pe  ly 


where  c is  the  eddy  diffusivity  of  momentum. 
Combining  equations  240  and  241  yields 


(pv).^=  (p*pc)0 


Note  that  the  bars  representing  time-averaged  quantities  have  been 
dropped  as  all  of  the  terms  are  mean  values  and  the  fluctuating  terms 
no  longer  appear. 

If  it  is  assumed  that  the  Prandtl  mixing  length  hypothesis  is  valid 
for  the  case  of  blowing,  the  eddy  momentum  diffusivity  can  be  written 
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Substituting  equation  243  into  equation  242  results  ina  non-linear 
second  order  differential  equation.  It  was  anticipated  that  a solution 
of  this  non-linear  boundary  value  problem  would  be  too  prohibitive  in 
terms  of  computer  time  to  be  of  value  in  the  erosive  burning  rate 
iteration.  Thus,  an  approximate  solution  of  equations  242  and  243  was 
deemed  appropriate  for  the  present  analysis 

The  governing  equation  can  be  approximated  (43,66) 

(pv)^  u = (P  + PE)  ^ 244 

for  the  case  of  constant  properties  where  the  no  slip  condition  at  the 
wall 

0 y = 0 u = 0 245 

has  been  used  and  is  the  shear  stress  at  the  wall. 

Substituting  equation  243  into  equation  244  results  in  the  follow- 
ing non-linear  first  order  ordinary  differential  equation  relating  the 
parallel  and  normal  components  of  velocity: 

PC'y^(^)^  + P ^ + (pv)^^u)  = 0 246 
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Equation  246  can  be  linearized  by  use  of  the  quadratic  equation 


du  _ 
dy  " 


2pc'y  2pc-y^ 


T + 
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(pv)^^u)/(pc'y^)] 
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Since  adverse  pressure  gradients  or  flow  reversal  are  never  con- 
sidered, negative  velocity  gradients  are  physically  unrealistic.  Thus, 
only  the  plus  sign  of  the  "plus  or  minus"  in  equation  247  needs  to  be 
considered.  Thus 


.2  , 1/2 

* pv)„u)/(pc'y2)j 
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2pc'y 


2pc’y" 


The  boundary  condition  for  equation  248  is 


0 y = 6 u = u 

or> 
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where  5 is  the  boundary  layer  thickness  and  is  the  free  stream 
velocity.  Obviously,  a numerical  solution  is  still  required.  However, 
by  assuming  a value  for  t.,,  the  above  system  can  bo  treated  as  an 
initial  value  problem  with  the  initial  condition 


0 y = 0 


^ = T /U 

dy  ” 
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A numerical  solution  of  equation  248  can  be  determined  b>  employing  a 
marching  scheme.  The  solution  must  be  detennined  at  several  down 
stream  positions  in  the  vicinity  of  the  position  of  interest  in  order 
that  a relation  between  the  shear  stress  at  the  wall  and  the  x position 
can  be  described.  The  resulting  velocity  profiles  can  then  be  used  to 
check  on  th6  assumed  value  for  the  surface  shear  stress  by  use  of  the 
von  Karman  momentum  integral  equation.  For  the  case  of  surface  blowing 
at  constant  pressure,  the  von  Karman  momentum  integral  can  be  written 


T + (cv)  u = j— 
w ' 'w  “ dx 


pu(u^  - u)dy 
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The  derivative  in  equation  251  can  be  determined  numerically  after 
performing  the  indicated  integration  at  several  locations  about  the 
position  of  interest.  If  the  calculated  value  for  the  wall  shear  does 
not  agree  with  the  assumed  value,  the  wall  shear  stress  must  be  deter- 
mined by  an  iteration  process. 

A less  rigorous  but  acceptable  alternative  to  the  iteration  process 
is  to  use  an  empirical  relation  for  the  wall  shear. 

At  this  point  it  is  convenient  to  define  the  skin  friction  coeffi- 
cient 
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For  the  case  of  no  blowing  the  skin  friction  coefficient  is  well 
represented  by  the  empirical  relation  (55), 
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ICf  = 0.0296  Re'°-^ 
0 


where  Re^^  is  the  Reynolds  number  based  on  down  stream  position.  The 
relation  for  skin  friction  given  by  equation  253  is  valid  for  Reynolds 
numbers  in  the  following  range 

5 X 10^  < Re  < 1 X 10^  25^ 

X 


Another  empirical  relation  for  the  skin  friction  coefficient  is  the 


Schul tz-Grunow  equation  (56) 


-2.584 


^ = 0.185  (fog^Q  ReJ  2 

' 0 

g 

which  is  valid  to  Reynolds  numbers  based  on  length  to  as  high  as  10  . 

The  reduction  in  shear  stress  at  the  wall  due  to  blowing  can  be 
expressed  by  the  empirical  relation  (43,57,58) 


C^/C^  = frt(l  + B)/B 

' 0 


where  is  the  local  skin  friction  coefficient  with  blowing  and  B is 
the  blowing  parameter  defined 


Figure  27  compares  the  prediction  of  equation  256  to  the  experimental 


data  of  Simpson  and  McQwaid  as  indicated  by  Mills  (59). 
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Combining  equations  256  and  257,  the  blowing  parameter  can  be 


written 

B = exp  C(pv)y(}c^  P„  u„)]  - 1 258 

0 

It  is  now  possible  to  obtain  an  empirical  relation  for  the  wall 
shear  stress  Rearranging  equation  252  and  combining  with  equations 
255,  256,  and  258  yields 

-2.584 

= i^u>(Pv)y[exp((pv)^/(0.185(Cog^QRe^)  u^P»))  * H 259 

An  approximate  solution  to  equation  248  can  now  be  formulated  in 
terms  of  an  empirical  expression  for  the  shear  stress  at  the  burning 
surface.  The  boundary  condition  given  by  equation  245  is  changed  to 

0 y = 0 xjv  260 

After  evaluating  the  wall  shear,  with  equation  259,  equation 
248  can  be  solved  by  a numerical  marching  scheme.  In  the  present 
analysis  a modified  Euler  predictor  corrector  (60)  scheme  is  employed. 

Figure  28  depicts  several  velocity  profiles  calculated  at  various 
blowing  rates  and  free  stream  velocities.  It  can  be  seen  that  as  the 
blowing  rate  is  increased  or  the  free  stream  velocity  is  decreased,  the 
velocity  gradient  in  the  region  near  the  burning  surface  decreases  (as 
would  be  expected).  It  is  apparent  that  the  behavior  of  the  predicted 
velocity  profile  is  at  least  qualitatively  correct. 

Having  obtained  a solution  for  the  velocity  profile,  the  eddy 
momentum  diffusivity  in  the  region  near  the  wall  is  computed  by 
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rhe  Boundary  Layer  Velocity  Profile 


employing  the  Prandtl  mixing  length  hypothesis 


e 


c'y 


du 

dy 


243 


Thus  the  eddy  diffusivity  varies  with  ordinate  in  the  flame  zone. 
An  average  value  for  the  eddy  diffusivity  in  the  flow  region  between 
the  propellant  surface  and  the  flame  of  interest  can  be  determined  by 
the  following  expression: 


(c)  = 1/L 

' 'av 


e dy 
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where  L is  the  pseudo  one  dimensional  flame  height.  Figure  29  depicts 
the  average  eddy  diffusivity  as  a function  of  distance  above  the  surface 
for  the  velocity  profiles  presented  in  Figure  28.  It  can  be  seen  that 
the  effect  of  blowing  at  the  surface  has  a dramatic  effect  on  the  cal- 
culated values  for  the  eddy  diffusivity.  Increased  blowing  by  means  of 
an  increase  in  propellant  burning  rate  or  by  means  of  a decrease  in  the 
cross  flow  velocity  results  in  a drastic  decrease  in  the  average  eddy 
diffusivity  in  the  region  near  the  propellant  surface. 

I The  effect  of  cross  flow  velocity  on  the  propellant  burning  rate 

[ can  now  be  determined.  It  is  assumed  that  in  the  normal,  nonerosive, 

combustion  regime,  the  transport  coefficients  of  thermal  conductivity 
and  mass  diffusion  are  the  molecular  ones.  However,  under  erosive 
conditions  the  transport  coefficients  are  modified  by  turbulent  com- 
ponents which  are  dependent  upon  the  nature  of  the  flow.  Moreover,  in 
the  present  analysis  it  is  assumed  that  the  ratio  of  the  molecular 

I 

transport  property  to  the  transport  property  in  the  presence  of 
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turbulence  is  proportional  to  the  ratio  o1  tin;  laminar  viseosit>  to  the 
average  value  of  the  total  viscosity  in  the  presence  of  turbulence. 

Thus 


W 

er~  a 

t (p  + pc) 
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Since  there  are  three  flames  considered  in  the  propellant  combustion 
model,  and,  in  general,  all  are  of  differing  heights  above  the  surface, 
separate  average  turbulent  transport  properties  must  be  calculated  for 
each  flame.  Thus,  there  are  turbulent  thermal  conductivities  associated 
with  the  AP  inonopropel lant  flame,  the  primary  flame  and  the  final  diffu- 
sion flame. 
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and  there  are  mass  diffusion  coeff iciei\ts  associated  with  the  (irimary 
flame  and  the  final  diffusion  flame 


Ppp  • kjj(l  + pc/u)P  266 

I^PF  ■ kp(l  *■  pf/p)P  2o7 


where  and  k^j  are  the  proportionality  constants  for  equation  262,  c 
is  the  average  eddy  diffusivity  of  momentum  evaluated  by  equation  262, 
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and  X and  Pare  the  molecular  thermal  conductivity  and  mass  diffusivity 
respectively. 

Since  the  PEM  burning  rate  calculation  is  an  iterative  solution 
of  several  simultaneous  equations,  a closed  form  solution  for  the  ero- 
sive burning  rate  is  not  possible.  Furthermore,  the  transport  properties 
used  in  the  burning  rate  calculation  are  a function  of  the  propellant 
burning  rate  with  a strong  coupling  between  the  propellant  blowing 
effect,  the  burning  rate  and  the  transport  property  calculation. 

An  iterative  solution  for  the  erosive  burning  rate  in  which  the 
transport  properties  were  averaged  between  the  values  for  consecutive 
iterations  was  found  to  be  successful.  The  scheme  employed  for  the 
solution  was  as  follows: 

1.  Calculate  the  burning  rate  using  the  molecular  transport 
properties. 

2.  With  this  burning  rate,  calculate  the  blowing  parameter 
and  the  velocity  profile  in  the  turbulent  boundary 
layer. 

3.  Using  the  calculated  velocity  profile,  calculate  the  transport 
properties  for  the  three  flames. 

4.  Averaging  the  latest  calculated  values  of  the  transport 
properties  with  the  values  computed  for  the  prior  iteration, 
recycle  to  step  1. 

The  process  described  above  is  repeated  until  the  desired  conver- 
gence is  reached.  Step  4 of  the  above  method  has  proved  to  be  extremely 
significant  in  obtaining  a solution  since  convergence  is  often  speeded 
up  by  a factor  of  ten  or  more  (61). 
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The  erosive  burning  model  described  above  was  derived  for  the 
purpose  of  velocity  coupling  calculations  and  was  not  expected  to  pre- 
dict quantitative  erosive  burning  results.  However,  as  will  be  pointed 
out  in  the  results  section,  the  model  has  several  unique  and  desirable 
features.  For  example,  the  model  predicts  a threshold  velocity. 
Furthermore,  it  predicts  the  correct  dependence  of  the  threshold 
velocity  on  the  combustion  pressure  and  burning  rate.  In  addition, 
the  model  is  able  to  predict  the  effect  of  oxidizer  particle  size  on 
the  erosive  burning  characteristics. 

The  PEM  Nonsteady  State  Combustion  Model 
As  discussed  in  the  introduction,  the  two  parameters  of  importance 
in  the  nonsteady  state  combustion  of  composite  propellants  are  the 
pressure  coupled  response  and  the  velocity  coupled  response.  Those 
parameters  were  defined  and  discussed  in  the  nonsteady  state  combustion 
section  of  the  literature  review. 

In  order  to  evaluate  the  nonsteady  state  burning  characteristics 
of  a composite  propellant,  a nonsteady  state  statistical  combustion 
model  was  derived  from  the  steady  state  PEM. 

Derivation  of  the  Statistical  Nonsteady  State  Model 
The  total  propellant  mass  flux  from  a burning  composite  propellant 
is 
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The  linear  perturbation  of  the  above  equation  with  respect  to  time 
yields 


■"t  " ^t  pQ.k  Ijj 


dD 
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where  the  primes  denote  perturbed  quantities. 
The  pressure  coupled  response  function  Is 
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Combining  equation  73  and  268  yields 
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where  is  the  pressure  coupled  response  of  the  pseudopropellant 
with  diameter  between  0^  and  + dD^  and  of  oxidizer  type  k. 
Similarly  for  the  velocity  coupled  response 


^ " Pt"*  ^ ^o\k  Ip  ^,d.k"'d,k  '"k^^’d.k 

^0 
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Equations  269  and  270  provide  the  formalism  for  computing  the  over- 
all composite  propellant  response  functions.  All  that  remains  is  to 
determine  the  pseudopropellant  response  functions.  Note  that  the  use 
of  equations  269  and  270  provides  the  contingency  for  predicting  oxidizer 
particle  size  effects  on  the  nonsteady  state  combustion.  Furthermore, 
the  equations  reduce  the  large  scale  heterogeneity  of  the  burning 
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surface  such  that,  to  a good  approximation,  a homogeneous  nonsteady 
model  can  be  employed  to  calculate  the  pscudopropel lant  pressure 
coupled  responses. 

The  Pseudopropellant  Pressure  Coupled  Response  Function 

Several  methods  are  available  for  calculating  the  pressure  coupled 
response  of  each  pseudopropellant:  (a)  a small  perturbation  technique 
can  be  applied  directly  to  the  steady  state  PEM  to  determine  the  pseudo- 
propellant response;  (b)  the  Denison  and  Baum  model  for  nonsteady 
burning  can  be  used  by  calculating  the  "A"  and  "B"  parameters  directly 
from  the  steady  state  combustion  model;  or  (c)  the  "A"  and  "B"  parameters 
of  the  Denison  and  Baum  model  can  be  calculated  from  the  Cohen  postulates; 
and  (d),  the  pseudopropellant  response  functions  can  be  calculated  using 
the  technique  of  Zeldovich  and  Novozhilov. 

With  the  exception  of  the  first  approach,  the  above  models  were 
discussed  in  the  literature  review.  The  first  approach,  that  of  per- 
turbing the  steady  state  PEM,  was  not  employed  in  the  present  research 
program.  Method  (b)  reduces  to  method  (d)  and,  therefore,  was  also  not 
included  in  the  present  program. 

The  Cohen  Postulate  Method.  For  calculations  utilizing  method  (c), 
the  Cohen  postulates  were  applied  to  the  Denison  and  Baum  pressure 
coupled  response  function  to  obtain  the  pseudopropellant  response  func- 
tions. The  expression  for  the  Denison  and  Baum  pressure  coupled  response 
is 


I (R^  + iR^)  = AB/[X  + (AA)  - (1  + A)  + AB] 


76 
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where  R^/n  and  R^/n  are  the  real  and  imaginary  parts  of  the  pressure 
coupled  response  normalized  with  respect  to  the  burning  rate  pressure 
exponent,  and  the  parameter  X is  a complex  function  of  the  nondimen- 
sional  frequency,  The  real  part  of  the  complex  frequency  is 

X = i {1  + (1/‘''2)[(1  + 7 

r c 


The  imaginary  part  is 


o 1/2  1/2 

X.  = (1/2.7)[(1  + 1K2^)  - U 


where 


78 


= kp^w/rti^Cp  79 

k is  the  thermal  conductivity  of  the  propellant,  is  the  total  pro- 
pellant density,  Cp  is  the  specific  heat  of  the  propellant  and  m is  the 
mass  flux  from  the  burning  surface. 

Separating  the  real  and  imaginary  parts  of  equation  76,  the  re- 
sulting relation  for  the  real  part  of  the  pressure  coupled  response  is 

= AB(xJ  + X^){(xJ  + A^)[X^  - (1  + A)  + AB]  + AX^}/ 

{[(xj  + X^)[X^  . (1  + A)  + AB]  + AX^]^  + 

2/1 

[(xj  + x5)[(Xr  - (1  + A)  ^ AB]2  + 2(X^  - x2)AX^ 

[(X^  - (1  + A))  + AB]]^} 
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The  frequency  at  which  the  response  function  peaks  and  the  magnitude  of 
the  peak  of  the  response  function  are  given  by  the  Cohen  postulates 
Thus 


I 


fp  = H r/Do  272 

R^/n  = Ga/DQ  2/3 

where  H and  6 are  proportionality  constants  and  at  is  the  mass  fraction 
of  oxidizer  in  a particular  particle  size  distribution  mode. 

Differentiating  equation  271  with  respect  to  nondimensional  fre- 
quency and  setting  the  result  equal  to  zero,  results  in  a fourth 
equation  valid  at  the  frequency  were  the  response  peaks.  Thus 

i ■ 0 

Thus,  the  system  reduces  to  a set  of  two  implicit  equations  with 
the  two  unknowns  A and  B.  The  system  can  be  solved  by  a numerical 
secant  iteration  technique  to  determine  A and  B for  each  pseudopropellant 
in  the  polydisperse  propellant.  The  resulting  pressure  coupled  response 
function  for  each  pseudopropellant  is  then  integrated  over  all  particle 
diameters  using  equation  269  to  determine  the  pressure  coupled  response 
of  a polydisperse  propellant  at  a specific  frequency.  This  integration 
is  repeated  at  finite  frequency  increments  to  determine  the  response  as 
a function  of  frequency. 

The  ZN  Technique.  Method  (d),  the  ZN  technique,  allows  the  pseudo- 
propellant response  functions  to  be  related  to  steady  state  propellant 
combustion  parameters.  The  ZN  technique  was  described  in  the  literature 
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review.  The  expression  for  the  pressure  coupled  response  given  by  the 
ZN  technique  is 


n + (nu  - pu)Z^ 

*^p  ' 1 - u + {v  - iu/n)Z^ 
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The  superscript  o indicates  steady  state,  T^  is  the  surface  tem- 
perature and  Tp  is  the  initial  propellant  temperature. 

Equation  107  can  now  be  used  to  determine  the  pseudopropellant 
pressure  couped  response  functions  by  evaluating  the  steady  state 
parameters  of  each  pseudopropellant  with  the  steady-state  polydisperse 
petite  ensemble  combustion  model.  The  polydisperse  propellant  pressure 
coupled  response  can  then  be  determined  using  equation  269. 


The  Pseudopropellant  Velocity 
Coupled  Response  Function 

The  subject  of  velocity  coupled  response,  although  not  new,  has 
received  very  little  theoretical  treatment  in  the  open  literature.  The 
present  state  of  the  art  of  understanding  the  mechanism  for  velocity 
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coupled  response  is  further  exemplified  by  the  following  statement 
taken  from  the  JANNAF  Erosive  Burning/Velocity  Coupling  Workshop  (62): 

Velocity  coupling  could  not  be  discussed  due  to  a lack  of  knowledge 
on  the  subject  and/or  the  unwillingness  of  researchers  working  on  the 

subject  to  present  their  results " 

The  work  of  Lengelle  (43)  is  one  of  the  few  comprehensive  treat- 
ments of  the  subject;  see  the  literature  review  section.  The  mechanism 
of  velocity  coupling  proposed  by  Lengelle  represents  a logical  and  phy- 
sically realistic  approach  to  the  velocity  coupling  problem.  Therefore 
a similar  approach  was  taken  in  the  present  research  effort.  The  only 
difference  between  the  present  approach  and  the  approach  taken  by 
Lengelle  is  that  the  present  model  is  based  on  a more  physically  real- 
istic and  more  complete  erosive  burning  model. 

Several  methods  are  available  for  determining  the  pseudopropellant 
velocity  coupled  response  by  the  Lengelle  philosophy.  The  most  rigorous 
approach  is  to  perturb  the  steady  state  PEM  erosive  burning  model,  and 
thus  follow  the  same  formalism  Lengelle  used  (},n  the  Summerfield 
(modified  for  erosive  burning  by  Lengelle)  combustion  model.  This 
approach  involves  a very  lengthy  analysis  and  was  not  included  in  the 
present  research  program  for  that  reason. 

Another  method  involves  the  utilization  of  the  ZN  technique.  This 
method  can  be  implemented  by  replacing  all  of  the  pressure  terms  in 
the  ZN  formalism  by  cross  flow  velocity.  The  velocity  coupled  response 
is  then  relate  to  the  following  steady  state  erosive  burning  parameters: 


The  present  approach  for  velocity  coupling  was  to  employ  the 
closed  form  relation  between  pressure  coupled  response  and  velocity 
coupled  response  in  the  presence  of  erosive  conditions  derived  by 
Lengelle.  It  is  realized  that  this  does  not  represent  an  advancement 
inthestate  of  the  art  of  the  modeling  of  the  velocity  coupled  response. 
However,  this  method  does  of^er  the  capability  to  estimate  the  oxidizer 
particle  size  effects  on  the  velocity  coupled  component  of  the  response 
function  (an  effect  which  has  not  been  calculated  previously). 

The  expression  wtiich  relates  pressure  and  velocity  coupled  response 
can  be  written 
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= 0.9i'.gAB/{[\  - 1 + A(l/A  - 1)  + ABC] 

- (1/3  + 0.9B)  AB/Rp)  ?80 

where  is  the  erosive  strength,  A is  the  Denison  and  Baum  constant 

which  is  dependent  on  the  solid  phase,  and  \ is  the  complex  frequency. 
The  parameter  C arose  in  the  Lengelle  formalism  due  to  the  inclusion 
of  surface  coupled  reactions.  Thus 

C = 1 + hg  + C(l/2  + Bg)RT5/Es  281 

where  R is  the  gas  constant,  is  the  surface  temperature  and  is 
the  surface  decomposition  activation  energy.  The  parameter  C is  re- 
lated to  the  surface  heat  of  pyrolysis  and  typically  has  a value  less 
than  unity. 

In  the  present  analysis  it  is  assumed  that  the  term  on  the  right 
in  equation  can  be  neglected  since  the  parameter  C is  less  than  unity 
and  the  surface  activation  energy  is  of  the  order  10^.  Thus,  equation 
i.'81  becomes 


C = 1 + Be 

The  above  assumption  is  necessary  since  the  Denison  and  Baum 
analysis  used  to  compute  the  pressure  coupled  response  did  not  include 
surface  coupled  reactions.  With  the  above  assumption,  the  Lengelle 
pressure  coupled  response  for  the  case  of  no  cross  flow  velocity  re- 
duces  to  the  Denison  and  Baum  pressure  coupled  response.  Thus,  by 
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utilizing  equation  282  for  the  parameter  C , equation  280  remains  com- 
patible with  the  pressure  coupled  response  calculation  in  the  limiting 
case  of  zero  cross  flow  velocity. 

The  Velocity  Coupled  Response  Solution  Methodology.  The  erosive 
burning  rate  is  first  calculated  by  the  erosive  burning  model  described 
in  a previous  section.  The  A and  B constants  in  equation  280  are  then 
determined  for  each  pseudopropellant  by  the  method  used  for  determining 
the  pressure  coupled  response.  The  pressure  coupled  response,  Rp,  of 
each  pseudopropel lant  is  computed  by  employing  the  Cohen  postulates  as 
in  the  case  of  no  erosive  burning.  The  erosive  strength  of  each  pseudo- 
propellant is  calculated  with  the  following  expression 


= (r  - r )/r 
e ' e n'  n 
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where  r^  is  the  pseudopropellant  burning  rate  in  the  presence  of  erosive 
conditions  and  r^  is  the  pseudopropellant  burning  rate  in  the  absence 
of  a cross  flow  velocity.  The  parameter  C is  also  calculated  for  each 
pseudopropellant. 

The  velocity  coupled  response  of  each  pseudopropel lant  can  then 
be  calculated  by  equation  280.  The  overall  polydisperse  propellant 
velocity  coupled  response  is  computed  by  integrating  the  pseudopropel lant 
velocity  coupled  response  over  all  pseudopropellants  by  equations  270. 
Thus 
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The  integration  indicated  above  is  repeated  at  finite  frequency  incre- 
ments to  determine  the  velocity  coupled  response  as  a function  of 
frequency. 
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ANALYTICAL  RESULTS 


Introduction 

The  results  of  calculations  with  the  theoretical  model  for  the 
combustion  of  composite  solid  propellants  are  divided  into  four  parts. 
They  are  the  steady  state  burning  rate,  the  erosive  burning  rate,  the 
pressure  coupled  response  and  the  velocity  coupled  response  results. 

The  results  of  calculations  with  the  steady  state  burning  rate  and  the 
pressure  coupled  response  model  were  compared  to  available  experimental 
data.  However,  experimental  data  containing  an  adequate  description  of 
the  propellant's  composition  in  terms  of  particle  size  and  size  distri- 
bution were  not  available  for  comparison  with  the  results  of  calculations 
for  the  erosive  burning  or  velocity  coupling  models. 

Steady  State  PEH  Burning  Rate  Results 

The  PEM,  not  unlike  the  earlier  BDP  model  and  other  burning  rate 
•models,  relies  on  the  numerical  value  of  several  constants  whose 
magnitude  are  known  only  to  a small  degree  of  precision.  The  constants 
include  activation  energies,  pre-exponential  factors,  reaction  orders, 
heats  of  reaction,  specific  heats,  thermal  conductivities,  etc.,  for 
the  various  flames,  gas  and  solid  phase  reactions. 

In  order  to  determine  the  appropriate  values  for  these  input  con- 
stants, the  experimental  burning  rate  data  presented  by  Miller  (6,7) 
was  used.  Miller  determined  the  burning  rate  and  pressure  exponent  for 
a series  of  AP/HTPB  non-metalized  propellants  at  a pressure  of  68.03  atm 
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and  an  initial  temperature  of  294®K.  These  composite  propellants  were 
formulated  from  eight  individual  oxidizer  grinds  with  the  total  oxidi- 
zer mass  fraction  held  constant  at  87.4  percent  for  each  blend.  Figure 
19  presents  the  particle  size  distribution  as  weight  percent  less 
than  diameter  versus  particle  diameter  for  some  of  the  oxidizer  grinds 
used  in  the  Miller  formulations.  The  solid  lines  through  the  data 
points  represent  a least  squares  fit  of  the  measured  particle  size 
distribution  to  a log  normal  distribution  function.  The  values  of  the 
weight  median  diameters,  IT,  and  the  distribution  width  parameter,  o, 
are  presented  in  Table  1 for  each  of  the  oxidizer  grinds.  Also  appear- 
ing in  Table  1 are  the  standard  deviations  of  error  for  the  log  noniial 
fit  of  each  of  the  grinds.  Thus,  the  particle  size  distributions  of  the 
eight  grinds  are  accurately  represented  by  the  log  normal  approximations. 


Table  1 

Correlation  for  the  Log  Normal  Distributions 


Nominal 

Diameter  (p) 

D(p) 

o(u) 

Stnd. 

Dev.  (%) 

400 

448.0 

1.222 

.775 

200 

195.0 

1.628 

3.079 

90 

71.0 

1.370 

2.382 

50 

44.2 

1.445 

.979 

20 

22.6 

1.676 

2.075 

10 

10.8 

1.308 

2.830 

6 

5.23 

1.878 

1.041 

2 

1.89 

1.305 

1.134 

0.7 

.686 

2.716 

2.204 

The  twonty-five  Miller  propel  l.uilb  jre  es‘jent.ii\l  ly  five  similar 
groups  of  five  families  of  propellants  (in  terms  of  the  oxidizer  par- 
ticle size  distribution).  The  families  range  from  very  broad  overall 
distributions  to  very  narrow  distributions,  all  with  the  same  percentage 
of  oxidizer.  The  compositional  grid  for  the  twenty-five  propellants 
are  presented  tn  Table  2.  The  total  percentage  of  oxidizer  for  each 
propellant  equals  87.4,  a constant  for  all  f the  propellants. 

The  burning  rate  and  the  burning  rate  pressure  exponents  deter- 
mined by  Miller  for  the  twenty- five  pro[iellants  are  presented  in 
Table  3.  Those  propellants  for  which  no  data  are  presented  are  pro- 
pellants that  presented  problems  during  the  propellant  mixing  operation. 

The  results  of  calculations  with  the  PF.M  are  presented  in  Figures 
30  and  31.  Figure  30  depicts  the  theoretical  burning  rate  versus  the 
experimental  burning  rate  for  the  twenty-one  propellant  configurations 
considered,  and  Figure  31  depicts  the  theoretical  pressure  exponent 
versus  the  experimental  exponent  for  the  same  propellants.  In  both 
cases,  the  PEM  prediction  matched  the  experimental  dat.i  quite  well  with 
over  seventj'-f i ve  percent  of  the  data  being  within  ten  percent  of  the 
predicted  values.  The  PEM  input  parameters  used  for  the  above 
calculations  are  presented  in  Table  4. 

In  addition  to  the  parameters  listed  in  Table  4,  the  flame  temper- 
atures for  the  primary  and  final  flames  were  calculati'd  with  the  NASA 
equilibrium  thermochemistry  program  (63).  The  flame  temperatures  are 
stored  in  the  PEM  computer  code  in  tabular  form  as  a fuKtion  of  combus- 
tion pressure  and  propellant  composition.  Quadratic  interpolating 
polynomials  are  used  to  retrieve  data  at  pressures  and  compositions 
between  the  tabular  grid  points. 
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Table  2 


Miller  Propellant  Compositional  Grid 


Propellant 

Designation 


Per  Cent  Nominal  Particle  Size 
400  200  90  50  20  6 2 0.7 


SD 

1 1 1-88-1 

45, 

.26 

31 

58 

10. 

53 

SO 

III-88-2 

31 

00 

13 

.68 

- 

42. 

11 

SD 

III-88-3 

- 

55, 

.79 

- 

31. 

58 

SD 

1 1 1-88-4 

31 

58 

- 

24, 

.21 

- 

31. 

58 

SD 

1 1 1-88-5 

42 

11 

- 

13 

,68 

“ 

31. 

58 

SD 

1 1 1-88-6 

31 

58 

13, 

.68 

31 

58 

10. 

53 

« 

SD 

1 1 1-88-7 

31 

.58 

13, 

.63 

42. 

11 

- 

SD 

1 1 1-88-8 

- 

31 

58 

24, 

,21 

31. 

58 

- 

SD 

1 1 1-88-9 

31 

58 

- 

24, 

.21 

31. 

58 

- 

SD 

in-88-10 

42 

21 

- 

13, 

,68 

31. 

58 

- 

SD 

III-88-11 

- 

45, 

,26 

42 

11 

- 

- 

SD 

III-88-12 

31 

.58 

13, 

,68 

42 

11 

- 

- 

SD 

III-88-13 

- 

55, 

,79 

31 

58 

- 

- 

SD 

III-88-14 

31 

58 

- 

24, 

.21 

31 

58 

- 

- 

SD 

III-88-15 

42 

11 

- 

13, 

.68 

31 

58 

- 

- 

SD 

rri-88-16 

31 

58 

31 

58 

24. 

.21 

- 

- 

SD 

III-88-17 

31 

.58 

00. 

,79 

- 

- 

SD 

III-88-18 

42 

.11 

45. 

,26 

- 

- 

SD 

III-88-19 

31 

58 

- 

55. 

,79 

- 

- 

SO 

1 1 1-88-20 

42 

11 

- 

45, 

.26 

- 

- 

SD 

III-88-21 

31 

58 

31 

58 

10 

53 

13. 

,68 

- 

SD 

1 1 1-88-22 

31 

58 

- 

42 

11 

13. 

,68 

- 

- 

SD 

1 1 1-88-23 

42 

,11 

- 

31 

58 

13. 

,68 

- 

- 

SD 

1 1 1-88-24 

31 

.58 

- 

42 

11 

13. 

,68 

- 

- 

SD 

III-88-25 

42 

11 

31 

58 

13. 

68 

* 
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Table  3 


Strand  Burning  Rate  of  the  Miller  Propellant  Compositions  at  1000  psia 


Propellant 

Designation 

(c/sac) 

Pressure 

Exponent 

SD  III-88-1 

SD  1 1 1-88-2 

2.977 

0.916 

SD  III-88-3 

3.636 

. 0.689 

SD  III-8S-4 

2.847 

0.797 

SD  1 1 1-88-5 

2.253 

0.928 

SD  111-88-6 

2.903 

0.621 

SD  III-88-7 

- 

- 

SD  1 1 1-88-8 

2.786 

0.692 

SD  1 1 1-88-9 

2.743 

0.771 

SD  1 1 1-88-10 

2.278 

0,841 

SD  1 1 1-88-11 

• 

- 

SD  111-88-12 

2.626 

0.617 

SD  III-88-13 

- 

- 

SD  III-88-14 

2.477 

0.613 

SD  III-88-15 

1.824 

0.690 

SD  III-88-16 

1.417 

0.451 

SD  III-88-17 

2.118 

0.474 

SD  III-88-18 

1.803 

0.437 

SD  III-88-19 

1.974 

0.529 

SD  1 1 1-88-20 

1,405 

0.610 

SD  1 1 1-88-21 

.828 

0.430 

SD  1 1 1-88-22 

1.316 

0.458 

SD  1 1 1-88-23 

1.171 

0.463 

SD  1 1 1-88-24 

1.364 

0.449 

SD  1 1 1-88-25 

1.120 

0.528 
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PEM  BURNING  RATE  (CM /SEC) 


EXPERIMENTAL  BURNING  RATE  (CM/SEC) 


Figure  30.  A Comparison  of  PEM  Burning  Rate  with  Experimental 
Burning  Rate:  Miller  Data 
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Table  4 


?' 


PEM  Input  Parameters  for  the  Miller  Propellants 


Description  of  Parameter  Value  Units 

Number  of  oxidizer  types  1 


Start  Pressure 

34.0136 

atm. 

Stop  pressure 

136.0544 

atm. 

Number  of  pressure  points 

3 

- 

Maximum  number  of  terms  in 
series  solution 

50 

Maximum  error  in  series  solution 

.0000001 

- 

Number  of  integration  steps  for  rate 
calculation 

40  - 50 

. 

Fuel  volume  parameter  from  Vf=CDQ 

3.0 

- 

Diameter  to  start  integration 

.1 

microns 

Diameter  to  stop  integration 

1100.0 

microns 

Fuel  binder 

HTPB 

- 

Initial  propellant  temperature 

294.15 

°K 

Mass  fraction  of  oxidizer 

.6737 

- 

Heat  of  pyrolysis  of  fuel 

433.0 

cal/gm 

Density  of  fuel 

.9030 

gm/cm’ 

Pre-exponential  for  fuel  surface 
decomposition 

4000 

gm/cm^-sec 

Activation  energy  for  fuel  surface 
decomposition 

16,900.0 

cal /mole 

Oxidizer  type 

AP 

- 

Molecular  weight  of  final 
flame  products 

20.78 

gm/mole 

Final  flame  stoichiometry  variable 

4.18 

- 

Primary  flame  stoichiometry  variable 

3.0 

- 

Molecular  weight  of  primary 
flame  products 

28.0 

gm/mole 

Oxidizer  heat  of  pyrolysis 

-120.0 

cal/gm 

Density  of  oxidizer 

1.95 

gm/cm’ 

Arrhenius  factor  for  oxidizer 
decomposition 

4 X 10^ 

gm/cm^-sec 

Activation  energy  for  oxidizer 
decomposition 

23000.0 

cal/mole 

AP  Flame  temperature  at  100  psia 

1388.0 

AP  Flame  temperature  at  2000  psia 

1413.0 
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Table  4 Continued 


I 


■ 

Description  of  Parameter 

Value 

Units 

0xidi2er  ignition  delay 
proportionality  factor 

190.0 

sec  atir/"/ 

Oxidizer  ignition  delay 
pressure  exponent  (m) 

.721 

cm^ 

Oxidizer  ignition  delay 
diameter  exponent  (x) 

.80 

* 

Activation  energy  for 
primary  flame 

15,000 

cal/mole 

Activation  energy  for  AP  flame 

29,000 

cal /mole 

Arrhenius  frequency  factor  for 
primary  flame 

188.0 

sec“^ 

Arrhenius  frequency  factor  for 

AP  flame 

70,000 

sec"' 

Primary  flame  reaction  order 

2.0 

- 

AP  flame  reaction  order 

1.8 

Specific  heat  of  solid  and  gases 

.3 

cal/gm  ^ 

Thermal  conductivity  of  gases 

.003 

cal/cm-sec- 

Diffusion  parameter 

7.6  X 10*^ 

Average  Diffusion  flame  standoff 
parameter 

.3 

- 
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The  effect  of  oxidizer  particle  size  on  the  polydisperso  propellant 
burning  rate  can  be  demonstrated  by  presenting  the  pseudopropellant  mean 
state  burning  rate  and  pressure  exponent  as  a function  of  oxidizer 
particle  size.  Figure  32  presents  pseudopropellant  burning  rate  versus 
oxidizer  diameter  and  Figure  33  presents  the  pseudopropellant  pressure 
exponent  versus  oxidizer  diameter. 

Erosive  Burning  Results 

The  PEM  erosive  burning  model,  not  unlike  the  PEM  steady  state 
burning  rate  model,  requires  a very  complete  description  of  the  propel- 
lant composition  in  order  to  make  quantitative  predictions.  Especially 
important  is  a complete  description  of  the  oxidizer  particle  size 
distribution. 

In  the  past,  experimental  results  have  not  attempted  to  systemati- 
cally correlate  erosive  burning  characteristics  with  propellant 
composition.  Furthermore,  current  erosive  burning  data  is  incomplete 
in  regard  to  a complete  description  of  the  propellant  composition. 
Therefore,  it  was  not  practical  to  attempt  a quantitative  comparison  of 
the  theoretical  erosive  burning  results  with  experimental  results.  How- 
ever, in  order  to  illustrate  that  the  PEM  erosive  burnTng  model  does 
indeed  predict  the  experimentally  observed  trends,  the  model  was  employed 
to  predict  the  erosive  burning  characteristics  of  a "computer  propellant". 

The  computer  propellant  considered  was  a composite  containing  87.4 
percent  ammonium  perchlorate  with  a HTPB  fuel  binder.  The  constants 
used  for  the  Miller  propellants  as  listed  in  Table  4 were  also  employed 
in  the  computer  propellant. 


Figure  32.  The  Pseudopropellant  Burning  Rate  Versus  Oxidizer  Particle  Diameter 


Figure  33.  The  Pseudopropellant  Burning  Rate  Pressure  Exponent  Versus 
Oxidizer  Particle  Diameter 


The  effect  of  free  stream  velocity  on  the  propellant  burnitig  rate 
is  depicted  in  Figures  34  and  35,  where  the  propellant  burning  rate  is 
plotted  as  a function  of  the  free  stream  velocity.  In  Figure  34  the 
combustion  pressure  was  68.03  atmospheres  while  in  Figure  35  the  pres- 
sure was  34.02  atmospheres. 

The  effect  of  burning  rate  on  the  threshold  velocity  is  presented 
in  Figure  36. 

The  effect  of  combustion  pressure  on  the  threshold  velocity  (the 
velocity  below  which  the  burning  rate  is  independent  of  the  velocity) 
is  illustiated  in  Figure  37  where  threshold  velocity  versus  combustion 
pressure  is  plotted. 

The  effect  of  oxidizer  particle  size  on  the  erosive  burning  rate 
is  illustrated  in  Figure  38  where  the  pseudopropel lant  burning  rate  is 
plotted  as  a function  of  oxidizer  diameter  for  several  free  stream 
velocities. 

Figure  39  presents  the  critical  particle  diameter  for  erosive 
burning  as  a function  of  free  stream  velocity. 

Pressure  Coupled  Response  Results 

Two  techniques  were  employed  to  calculate  the  pressure  coupled 
response,  the  Zeldovich  Novozhilov  (ZN)  method  and  the  Cohcn-Denison 
and  Baum  method.  Both  methods  were  compared  to  experimental  data  for 
several  propellants. 

One  of  the  most  extensively  characterized  propellants,  in  terms  of 
nonsteady  combustion,  is  A-13  propellant.  That  propellant  is  composed 
of  20.4  percent  PBAN,  76.0  percent  ammonium  perchlorate  and  3.6  percent 
epor  828  resin.  The  ammf'niurr  perchlorate  consists  of  monomodal 
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Figure  34.  Propellant  Burning  Rate  Versus  Cross  Flow  Velocity  at  1000  psia 


threshold  velocity  (M/SEC) 
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iqure  37.  The  Effect  of  Pressure  on  the  Threshold 


CRITICAL  PARTICLE  DIAMETER  (MICRONS) 


particle  size  distribution  with  a 90  micron  nominal  weight  median 


diameter.  Experimental  measurements  of  the  actual  oxidizer  particle 
size  distribution  have  not  been  reported.  However,  normal  90  micron 
grinds  typically  have  a distribution  width  parameter  of  1.4  and  a 
weight  median  diameter  of  71  microns. 

The  above  particle  size  distribution  parameters  were  used  for  the 
theoretical  prediction  of  the  pressure  coupled  response  of  the  A-13 
propellant. 

Figure  40  compares  the  results  of  calculations  of  the  pressure 
coupled  response  of  A-13  propellant  utilizing  the  Cohen-Denison  and 
Baum  method  to  the  experimental  results  (48)  for  that  propellant.  The 
proportionality  constants  in  equations  83  and  84  used  for  this  calcu- 
lation were  those  used  by  Cohen  in  reference  (1). 

Figure  41  presents  a comparison  of  the  results  of  calculations 
using  the  ZN  technique  to  the  experimental  data  for  A-13  propellant. 

Another  well  characterized  propellant  is  NWC-SP-520.  Pressure 
coupled  response  data  over  a wide  frequency  range  as  well  as  complete 
formulation  data  were  available  (64)  for  this  propellant. 

NWC-SP-520  propellant  is  composed  of  13.1  percent  HTPB  binder  and 
86.9  percent  ammonium  perchlorate.  The  ammonium  perchlorate  consisted 
of  a trimodal  blend  of  oxidizer  grinds.  The  oxidizer  particle  size 
distribution  data  is  in  Table  5. 
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Figure  40.  The  Real  Part  of  the  Normalized  Response  Fu 
Comparison  of  the  Cohen  - Denison  and  Baum 


EXPERIMENTAL 


Figure  41.  The  Real  Part  of  the  Normalized  Response  Function  Versus  Nondimensional  Frequency 
Comparison  of  ZN  Method  to  the  Experimental  Results  of  A-13  Propellant 


Table  5 


Oxidizer  Particle  Size  Distribution  Data  for  NWC-SP-520 


Nominal 

Diameter 

Weight 

Percent  in 
Propellant 

Weight 

Median 

Diameter 

Distribution 

Width 

Parameter 

200  31.9  170  1.426 

45  25.0  44  1.481 

6 30.0  4.85  1.85 


Figures  42  and  43  present  a comparison  of  the  theoretical  results 
of  the  Cohen-Denison  and  Baum  and  the  ZN  techniques  respectively  to  the 
experimental  data  for  the  NWC-SP-520  propellant. 

A limited  amount  of  pressure  coupled  response  data  are  available 
for  the  four  HTPB/AP  propellants  of  Miller  (6).  The  propellants  con- 
tained 87.37  percent  ammonium  perchlorate  and  12.63  percent  HTPB  binder. 
All  of  the  propellants  contained  a trimodal  blend  of  oxidizer  grinds. 

The  oxidizer  distribution  grid  for  the  four  propellants  is  listed  in 
Table  6. 


Table  6 

Oxidizer  Composition  of  Miller  Acoustic  Propellants 


Nominal 

Diameter 

50DC-870 

50DC-871 

50DC-872 

50DC-873 

400 

51.55 

47.62 

43.68 

200 

63.61 

90 

31.45 

50 

11.36 

12.93 

20 

21.84 

28.39 

10 

10.83 

6 

12.24 

2 

13.98 
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The  particle  size  data  for  the  oxidizer  grinds  in  the  above  pro- 
pellants are  presented  in  Table  1. 

Figure  44  presents  the  results  of  pressure  coupled  response  calcu- 
lations for  the  four  propellants  using  the  Cohen-Denison  and  Baum 
technique  while  Figure  45  presents  the  results  of  the  same  calculations 
using  the  ZN  technique.  Only  one  experimental  data  point  was  available 
for  each  propellant. 

In  an  attempt  to  illustrate  the  dramatic  effect  the  oxidizer  par- 
ticle size  distribution  can  have  on  the  pressure  coupled  response, 
theoretical  calculations  were  performed  for  two  "computer  propellants". 

The  propellants  were  representative  of  candidate  propellants  for  the 
Maverick  missile  propulsion  system.  Both  of  the  propellants  were  com- 
posed of  14  percent  HTPB  binder  and  86  percent  ammonium  perchlorate  having 
a bimodal  blend;  70  percent  of  the  blend  being  of  192  micron  weight 
median  diameter  and  30  percent  being  of  16.6  micron  weight  median  dia- 
meter. The  only  difference  between  the  two  propellants  is  in  the  width 
of  the  particle  size  distributions  in  the  two  distribution  modes  con- 
tained in  the  overall  distribution.  One  propellant  had  distribution 
width  parameters  of  2.0  for  the  16.6  micron  grind  and  2.0  for  the  192 
micron  grind  while  the  other  propellant  had  distribution  width  parameters 
of  4.0  for  the  16.6  micron  grind  and  1.6  for  the  192  micron  grind. 

The  results  of  the  theoretical  calculations  for  the  above  two  pro- 
pellants are  presented  in  Figure  46. 

The  effect  of  particle  size  on  pressure  coupled  response  is  further 
exemplified  in  Figure  47  where  the  real  part  of  the  pseudopropol lant 
response  is  presented  as  a function  of  oxidizer  particle  size  for 
several  frequencies. 
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Figure  42.  Comparison  of  the  Theoretical  Results  with  the  Cohen  - Denison  and  Baum  Method  to 
Experimental  Pressure  Coupled  Response  Data  of  NWC-SP-520  Propellant 


870 


to  the  Experimental  Results  for  the  Miller  Propellants 


Figure  46.  The  Effect  of  Oxidizer  Particle  Size  Distribution 
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Velocity  Coupled  Response  Results 

At  present  there  have  been  no  experimental  results  for  the  velocity 
coupled  response  of  composite  propellants  reported  In  the  open  liter- 
ature. Therefore.  It  was  not  possible  to  attempt  a quantitative 
comparison  of  the  theoretical  velocity  coupled  response  results  with 
experimental  results. 

Figure  48  presents  the  results  of  velocity  coupled  response  cal- 
culations for  the  four  previously  described  Miller  propellants. 

The  effect  of  the  oxidizer  particle  size  on  the  velocity  coupled 
response  Is  Illustrated  In  Figure  49  where  the  pseudopropellant  velocity 
coupled  response  Is  plotted  as  a function  of  oxidizer  particle  diameter 
and  frequency. 

The  effect  of  cross  flow  velocity  on  the  velocity  coupled  response 
of  propellant  50DC-870  Is  illustrated  in  Figure  50. 
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Figure  48.  The  Velocity  Coupled  Response  for  the  Killer  Acou 


Figure  49.  The  Effect  of  Oxidizer  Particle  Size  on  the  Velocity  Coupled  Response 
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DISCUSSION  OF  RESULTS 


Introduction 

The  discussion  of  the  previously  presented  results  is  divided  in- 
to four  parts.  They  include  a discussion  of  the  results  for  the  steady 
state  burning  rate,  the  erosive  burning  rate,  the  pressure  coupled 
response,  and  the  velocity  coupled  response. 

The  Steady  State  Burning  Rate 

Comparison  with  Experimental  Results 

The  comparison  of  the  results  of  calculations  with  the  steady  state 
PEM  to  the  experimental  results  of  Miller  (6,7)  represents  a very  com- 
prehensive test  of  the  PEM's  ability  to  predict  steady  state  burning 
characteristics  of  polydisperse  composite  propellants.  The  Miller  pro- 
pellants spanned  a wide  range  of  burning  rates  and  pressure  exponents 
along  with  a wide  range  of  oxidizer  particle  size  distributions. 

As  is  indicated  Figures  30  and  31,  the  PEM  did  a remarkable  job 
of  predicting  the  oxidizer  particle  size  effects  on  the  steady  state 
burning  rate  and  pressure  exponent. 

The  propellants  for  which  the  theoretical  results  were  in  the 
poorest  agreenwnt  with  the  experimental  results  were,  in  general,  those 
propellants  which  contained  a very  wide  overall  distribution  (for  example, 
propellants  SO  III-88-5,  SO  IlI-88-lO  and  SO  111-88-15).  Also,  the 
(h'^oretical  results  for  those  prope’ lants  cnntainini  100  micron  AP 


were.  In  general,  in  poorer  agreement  with  the  experimental  results. 

The  PEM  exponent  prediction,  on  the  other  hand,  was  in  good  agreement 
with  the  experimental  results  for  nearly  every  propellant. 

One  possible  cause  for  the  minor  discrepancies  in  the  burning  rate 
prediction  Is  that  the  geometric  considerations  In  the  PEM  do  not 
account  for  the  micro  scale  surface  roughness  on  the  burning  surface. 
That  Is,  for  the  case  of  a wide  overall  particle  size  distribution,  one 
would  expect  the  burning  surface  geometric  structure  (on  a micro 
scale)  to  be  much  rougher  than  for  the  case  of  a narrow  overall  distri- 
bution. This  is  due  to  the  wide  variation  of  burning  rate  with 
oxidizer  particle  size  and  the  higher  probability  that  a large  particle 
will  be  surrounded  by  a large  popjlation  of  small  particles  (In  the 
case  of  a wide  overall  distribution).  The  rougher  micro  surface  struc- 
ture would  result  in  an  increase  In  mass  flux  due  to  the  increase  in 
the  ratio  of  the  total  burning  surface  area  to  the  planar  surface  area. 

Since  the  PEM  does  not  presently  account  for  the  roughness  of  the 
burning  surface,  one  would  expect  the  PEM  to  underpredict  the  burning 
rate  of  propellants  which  contain  a wide  overall  particle  size  distri- 
bution. That,  however,  was  not  the  case.  The  PEM  actually  over- 
predicted the  burning  rate  for  the  two  propellants  with  the  widest 
overall  particle  size  distribution  (SD  III-88-5  and  SD  III-88-10). 

Thus,  it  may  be  concluded  that  the  surface  roughness  effects  are  of  a 
secondary  nature  In  the  burning  rate  determination. 

Another  possible  explanation  for  the  discrepancies  Is  due  to  the 
assumption  in  the  PEM  that  the  pseudopropellants  or  unit  flames,  which 
compose  the  burning  surface,  burn  Independently  of  the  surrounding 


pseudopropellants.  It  Is  obvious  that  for  the  case  of  a wide  overall 
distribution,  there  is  a higher  probability  that  the  large  particles 
In  the  distribution  will  be  surrounded  by  a large  population  of  small 
particles.  This  situation  would  Increase  the  probability  for  flame 
Interactions  between  adjacent  pseudopropellants.  Since  flame  inter- 
actions are  not  accounted  for  In  the  PEM  one  would  expect  the 
theoretical  prediction  to  be  In  poorer  agreement  with  the  experimental 
results  In  the  case  of  propellants  with  wide  overall  distributions. 

This  was  Indeed  the  case. 

It  should  be  noted  that  the  values  for  the  constants  used  in  the 
PEM,  as  presented  In  Table  4,  all  represent  physically  realistic 
values.  Furthermore,  most  of  the  values  have  previously  been  reported 
In  the  literature.  One  constant  used  In  the  PEM  for  the  comparison  with 
the  Miller  data,  however,  was  not  a realistic  one.  Its  value  was 
selected  to  achieve  a good  agreement  of  the  theoretical  prediction  with 
the  experimental  results.  That  constant  Is  the  particle  diameter  ex- 
ponent for  the  amount  of  fuel  associated  with  a particular  oxidizer 
particle  (see  equation  127). 

For  an  ordered  packing  of  monodlsperse  particles  it  can  be 
shown  (21)  that  the  value  for  the  diameter  exponent  is  two.  In  the 
case  of  a randomly  packed  polydispersion  of  spheres  (the  case  for  a 
composite  propellant)  the  value  of  the  diameter  exponent  is  expected 
to  be  near  two  but  certainly  between  values  of  one  and  three.  The 
effect  of  the  diameter  exponent  on  the  pseudopropel lant  properties  of 
density,  oxidizer  mass  fraction  and  oxidizer  volume  fraction  are 
depicted  in  Figures  51,  52.  and  53  respectively.  The  effect 


Figure  51.  The  Pseudopropellant  Density  Versus  Oxidizer  Particle  Diameter 


Figure  52.  The  Pseudopropellant  Oxidizer  Mass  Fraction  Versus  Particle  Size 


of  the  diameter  exponent  on  burning  rate  is  dramatic.  A value  of  two 
for  the  diameter  exponent  (the  expected  physically  realistic  value) 
ca.ses  the  small  particles  in  the  overall  distribution  to  burn  very  fuel 
rich  while  the  large  particles  burn  fuel  lean.  This  in  turn  causes  a 
dramatic  change  in  the  shape  of  the  curve  which  represents  pseudopropel- 
lant burning  rate  versus  oxidizer  particle  diameter.  Instead  of  the 
S-shaped  curve,  as  presented  in  Figure  32,  the  burning  rate  of  the  small 
particles  starts  at  a lower  value,  increases  with  particle  diameter  and 
reaches  a maximum  at  some  intermediate  diameter  and  then  decreases  with 
increasing  particle  size  to  a value  slightly  higher  (for  the  largest 
particle  in  the  distribution)  than  the  value  shown  in  Figure  32. 

In  order  to  fit  the  experimental  results  of  Miller  it  was  necessary 
for  the  diameter  exponent  to  be  a value  of  three.  Thus  all  of  the 
pseudopropellants  in  each  of  the  Miller  propellant's  particle  size  dis- 
tributions took  on  the  same  density,  oxidizer  mass  fraction  and  oxidizer 
volume  fraction  as  the  polydisperse  propellant. 

Physically  this  indicates  that  there  are  indeed  interactions  between 
the  flames  of  adjacent  pseudopropellants.  Thus,  although  on  a geometric 
basis  the  small  particles  iiave  proportionally  more  fuel  associated  with 
them  than  the  large  particles  in  an  oxidizer  particle  distribution,  much 
of  that  fuel  is  given  up  to  the  larger  particles  in  the  distribution 
during  the  gas  phase  mixing  process. 

Oxidizer  Particle  Size  Effects 

Figure  32  depicts  the  pseudopropellant  burning  rate  as  a function 
of  particle  diameters.  The  general  S-shaped  curve  as  shown  in  Figure 
32  is  representative  of  both  experimental  and  theoretical  burning  rate 
-r^Ml ts. 
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The  small  particles  yield  high  burning  rates.  As  particle  size 
decreases  the  burning  rate  approaches  a limiting  value  determined  by 
the  premixed  kinetic  reaction  standoff  distance  of  a kinetically 
controlled  primary  flame  between  the  oxidizer  and  the  fuel  binder. 
Conversely,  the  large  particles  yield  low  burning  rates.  As  particle 
size  increases,  the  burning  rate  approaches  a limiting  value  determined 
by  the  kinetically  controlled  AP  monopropellant  flame.  The  burning 
rate  for  the  intermediate  diameters  is  determin^'d  by  competing  processes 
within  each  of  the  three  flames  above  an  individual  oxidizer  particle. 

Figure  33  depicts  the  pseudopropellant  burning  rate  pressure 
exponent  as  a function  of  oxidizer  particle  diameter.  The  exponent 
decreases  from  a value  of  near  unity  for  the  submicron  size  oxidizer 
to  a minimum  of  about  .34  for  intermediate  oxidizer  particles  of 
diameters  between  10  and  100  microns.  Above  iOO  microns  the  pressure 
exponent  begins  to  increase  approaching  near  unity  for  very  large 
oxidizer  particles. 

To  explain  the  pressure  exponent  phenomena  exhibited  by  the  pseudo- 
propellants,  it  is  necessary  to  observe  the  controlling  combustion 
mechanism  for  the  given  oxidizer  particle  size.  For  small  particles 
the  controlling  flame  is  the  primary  flame.  The  limiting  condition  for 
small  particles  is  a kinetically  controlled  primary  flame.  The  physio- 
chemical  model  then  reduces  to  the  following  four  equations: 

"ox  ■ ‘ox  “P  <-Eox^'*^> 
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The  above  four  equations  can  be  differentiated  with  respect  to 
pressure  to  yield  a closed  form  solution  for  the  pressure  exponent,  n. 
The  result  after  simplification  is: 
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Thus,  the  limiting  value  of  the  pressure  exponent  for  small 
particles  is  determined  by  the  value  of  the  term  inside  the  brackets 
times  the  primary  flame  reaction  order.  The  value  of  that  term  is 
determined  by  the  heat  released  in  the  primary  flame,  the  specific 
heat  of  the  combustion  gases,  the  activation  energy  for  the  oxidizer 
decomposition  process  and  the  primary  flame  standoff  distance.  For 
typical  values  of  the  above  parameters,  the  term  inside  the  brackets 
has  a value  of  about  .5  and  the  pressure  exponent  is  seen  to  be  a little 
less  than  half  of  the  primary  flame  gas  phase  reaction  order,  or  of  the 
order  unity. 
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The  same  reasoning  can  be  applied  to  the  limiting  condition  of 
large  oxidizer  particles.  For  large  particles  the  controlling  mechanism 
is  the  kinetically  controlled  AP  monopropellant  flame,  as  the  character- 
istic times  for  diffusion  are  so  great  that  the  heat  feedback  from 
the  final  diffusion  flame  is  negligible.  The  physiochemical  model  then 
reduces  to  the  following  three  equations: 

V ■ *ox  '■‘P 

«Ip- 

The  above  three  equations  can  be  differentiated  with  respect  to 
pressure  to  yield  a closed  form  solution  for  the  pressure  exponent,  n. 
The  result  after  simplification  is: 

n ■ V[12  * 25° 

Thus,  the  limiting  value  of  the  pressure  exponent  for  large 
particles  is  analogous  to  the  limiting  condition  for  small  particles. 

The  limiting  value  of  the  exponent  for  large  particles  is  determined  by 
the  value  of  the  term  inside  the  brackets  times  the  ammonium  perchlorate 
monopropellant  flame  reaction  order.  For  typical  values  of  the  para- 
meters in  equation  290  the  term  inside  the  brackets  has  a value  of 
about  .5.  Thus,  the  pressure  exponent  is  again  seen  to  be  a little 
less  than  half  of  the  AP  flame  monopropellant  reaction  order,  or  of 
order  unity. 
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The  pressure  exponent  for  the  intermediate  particle  sizes  is  not 
so  easy  to  characterize  since  the  controlling  mechanism  for  the  com- 


bustion of  the  intermediate  particles  is  a complicated  competitive 
process  between  the  three  flame  zones.  The  surface  geometry  also 
changes  with  pressure  for  the  intermediate  particles  to  further  compli- 
cate the  process.  The  value  of  the  exponent  for  the  case  of  competing 
flames  is  partially  diffusion  controlled  and  the  diffusion  standoff 
distance  is  only  a very  weak  function  of  pressure.  Thus,  the  pressure 
exponent  is  lower  for  the  case  where  diffusion  takes  a part  in  the 
controlling  mechanism.  The  higher  the  dominance  of  diffusion  in  the 
controlling  combustion  process  the  lower  will  be  the  value  of  the 

pressure  exponent.  | 

Thus,  it  is  obvious  that  the  oxidizer  particle  size  has  a pro- 
nounced effect  on  the  steady  state  burning  characteristics.  Furthernwre, 
the  burning  characteristics  can  be  tailored  by  appropriate  adjustments 
in  the  oxidizer  particle  size  distribution. 

The  Erosive  Burning  Rate 

Comparison  with  Experimental  Observations 
As  stated  previously,  the  primary  objective  of  performing  calcula- 
tions with  the  erosive  PEM  was  to  illustrate  that  the  model  does  indeed 
predict  the  experimentally  observed  trends. 

One  of  the  more  predominant  experimental  observations  is  that,  as 
the  flow  velocity  across  a burning  propellant  surface  is  increased,  no 
effect  on  the  burning  rate  is  observed  until  the  cross  flow  velocity 
reaches  a certain  value  (the  threshold  velocity).  Figures  34  and  35 
illustrate  that  the  model  does  indeed  predict  a thresnold  velocity. 
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Physically,  the  threshold  velocity  is  due  to  the  interaction 
between  blowing  at  the  burning  surface  and  the  cross  flow  induced 
turbulent  boundary  layer.  At  low  cross  flow  velocities  the  blowing 
parameter  is  very  high.  Consequently,  the  shear  stress  at  the  surface 
is  very  low  and  the  velocity  gradient  in  the  flame  zone  region  directly 
above  the  surface  is  very  low.  Thus,  the  contribution  of  turbulence  to 
the  transport  properties  in  the  gas  phase  reaction  zone  is  negligible. 

The  above  situation  is  depicted  in  Figure  54a. 

As  the  cross  flow  velocity  is  increased,  a point  is  finally  reached 
where  the  effect  of  turbulence  is  finally  felt  in  the  gas  phase  reaction 
zone  as  illustrated  in  Figure  54b. 

Further  increases  in  the  cross  flow  velocity  result  in  increases 
in  the  transport  properties  within  the  gas  phase  reaction  zone  due  to 
the  increased  level  of  turbulence  in  that  region. 

The  effect  of  burning  rate  on  the  threshold  velocity  has  been  widely 
observed.  Slower  burning  propellants  have  lower  threshold  velocities 
than  higher  burning  rate  propellants.  In  order  to  illustrate  that  the 
erosive  PEM  predicts  the  correct  effect  of  burning  rate  on  the  threshold 
velocity,  it  was  necessary  to  compare  the  threshold  velocities  of 
propellants  which  contain  identical  oxidizer  particle  size  distribution 
as  well  as  Identical  composition;  thereby  eliminating  the  possible 
effects  of  flame  temperature  and  oxidizer  particle  size  on  the  thres- 
hold velocity.  The  computer  propellants  chosen  for  the  illustration 
were  all  identical  in  composition  and  particle  size  distribution. 

Changes  in  the  burning  rate  of  the  base  line  propellant  were  accom- 
plished by  simulating  catalysis  of  the  primary  and  AP  flames.  The 
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Figure  54. 


Gas  Phase  Reaction  Zone  in  the  Turbulent  Boundary  Layer 


206 


r 


r ^ 


catalysis  was  simulated  by  adjusting  the  activation  energies  in  the 
Arrhenius  expressions  for  the  primary  and  AP  flame  reaction  rate  con- 
stants. The  results  presented  in  Figure  36  illustrate  that  the  model 
predicts  the  correct  dependence  of  threshold  velocity  on  the  propellant 
burning  rate. 

The  physical  effect  of  burning  rate  on  the  threshold  velocity  can 
be  explained  solely  by  the  effect  of  blowing  on  the  turbulent  boundary 
layer.  For  low  burning  rates  the  blowing  rate  is  such  that  the  turbu- 
lence penetrates  into  the  gas  phase  reaction  zone  at  lower  cross  flow 

velocities  than  in  the  case  of  high  burning  rate.  As  the  burning  rate  is  j 

increased,  the  blowing  parameter  increases.  This,  in  turn,  increases  j 

the  difficulty  for  turbu''ence  penetration  into  the  flame  zone. 

The  effect  of  pressure  on  the  threshold  velocity  would,  at  first 
thought,  be  expected  to  be  analogous  to  the  effect  of  burning  rate  on 
the  threshold  velocity.  However,  experimental  evidence  (65)  suggests 
that  this  is  not  the  case.  Thus,  the  threshold  velocity  has  been 
observed  to  decrease  with  increasing  pressure.  The  same  effect  was  also 
exhibited  by  the  erosive  PEM  as  illustrated  in  Figure  37.  Thus,  even 

though  the  blowing  increases  with  increased  pressure,  the  effect  of 

pressure  on  the  turbulent  boundary  layer  is  such  that  the  shear  stress 

(or  velocity  gradient)  in  the  gas  phase  reaction  zone  increases  with 
pressure  at  a faster  rate  than  the  increased  blowing  can  decrease  the 
shear  stress. 

Oxidizer  Particle  Size  Effects 

The  effect  of  oxidizer  particle  size  on  the  erosive  burning  char- 
acteristics was  illustrated  in  Figure  38. 
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Physically,  as  the  free  stream  velocity  increases,  the  turbulence 
extends  closer  to  the  burning  surface.  Now,  the  combustion  of  the 
large  particles  is  normally  dominated  by  the  kinetically  controlled  AP 
monopropellant  flame  (since  the  characteristic  diffusion  times  are  so 
large).  However,  as  the  free  stream  cross  flow  velocity  is  increased 
and  the  turbulence  finally  extends  into  the  final  flame  zone,  the 
average  value  of  the  transport  properties  for  the  final  flame  are 
enhanced.  This  decreases  the  characteristic  diffusion  time  and  in- 
creases the  heat  transfer  from  the  final  flame  zone.  The  kinetically 
controlled  AP  monopropellant  flame,  on  the  other  hand,  is  driven  away 
from  the  surface  due  to  the  increased  mass  flow.  Thus,  as  the  free 
stream  cross  flow  velocity  is  increased,  the  kinetically  controlled  AP 
flame  becomes  less  dominant  in  the  controlling  mechanism  and  the  final 
flame  becomes  more  dominant.  A point  is  finally  reached  where  the  AP 
flame  has  extended  into  the  turbulence  region.  At  this  point  the 
controlling  mechanism  is  divided  between  the  AP  flame  and  the  final 
flame. 

As  particle  size  is  decreased,  the  characteristic  time  for  diffusion 
becomes  smaller.  Thus,  the  degree  of  extension  of  the  final  flame  into 
the  turbulent  region  decreases  with  particle  size. 

The  critical  particle  diameter  for  erosive  burning  is  presented  as 
a function  of  free  stream  velocity  in  Figure  39  (for  a nominal  burning 
rate  of  1.3  centimeters  per  second  and  a pressure  of  68  atmospheres). 
Thus,  for  a free  stream  cross  flow  velocity  of  300  meters  per  second 
and  a burning  rate  of  1.3  centimeters  per  second,  there  would  be  no 
erosive  burning  at  a pressure  of  68  atmospheres  if  all  of  the  oxidizer 
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particles  In  the  propellants  particle  size  distribution  were  smaller 
than  100  microns  in  diameter. 

The  Pressure  Coupled  Response 

Comparison  with  Experimental  Results 

The  two  methods  employed  to  calculate  the  pressure  coupled  response 
were  the  ZN  technique  and  the  Cohen-Denison  and  Baum  technique.  The 
latter  technique  utilizes  two  empirical  constants  (the  proportionality 
constants  for  equations  83  and  84).  The  ZN  technique,  on  the  other  hand, 
enables  the  pressure  coupled  response  to  be  calculated  from  first 
principles.  Both  methods  predicted  the  correct  order  of  magnitude  of  the 
pressure  coupled  response  for  each  of  the  propellants.  See  Figures  40 
through  45  where  the  comparisons  with  the  experimental  data  are  presented. 

In  those  propellants  that  contained  a multimodal  distribution  of 
particle  sizes,  both  methods  predicted  multiple  peaks  in  the  response 
function  (Figures  42  through  45).  The  Cohen-Denison  and  Baum  technique 
predicted  the  correct  magnitude  of  the  pressure  coupled  response,  at 
least  up  to  a frequency  of  4000  hertz.  The  ZN  technique,  on  the  other 
hand,  appears  to  underpredict  the  magnitude  of  the  response.  The  latter 
technique  did,  however,  correctly  predict  the  frequency  at  which  the  ex- 
perimental data  indicated  a maximum  in  the  response  function  for  both 
A- 13  and  NWC  SP-520  propellants. 

It  is  interesting  to  note  that  the  Cohen-Denison  and  Baum  technique 
predicted  the  pressure  coupled  response  of  the  four  Miller  acoustic  pro- 
pellants (see  Figure  44)  to  well  within  the  experimental  error  in  the 
data.  However,  all  four  of  the  propellants  contained  400  micron  AP, 
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which  proved  to  be  troublesome  during  the  steady  state  burning  rate 
predictions. 

Oxidizer  Particle  Size  Effects 

The  possible  dramatic  effect  of  the  oxidizer  particle  size  distri- 
bution on  the  pressure  coupled  response  is  illustrated  in  Figure  46. 

The  particle  size  effect  can  be  explained  by  looking  at  the  response 
of  each  pseudopropellant  in  the  oxidizer  particle  size  distribution. 

Figure  47  presents  the  pseudopropellant  pressure  coupled  response  as  a 
function  of  oxidizer  particle  diameter  for  several  perturbation  fre- 
quencies. It  is  seen  that  the  large  particles  resonate  and  therefore 
dominate  the  response  at  low  perturbation  frequencies.  As  the  frequency 
is  increased,  the  resonance  moves  to  smaller  particle  diameters.  Further- 
more, the  magnitude  of  the  response  at  the  resonant  particle  diameter 
increases  as  the  frequency  increases.  Thus,  it  is  the  small  oxidizer 
particles  in  the  distribution  which  drive  up  the  magnitude  of  the  pres- 
sure coupled  response  and  the  effect  of  the  small  particles  is  felt  to  a 
greater  degree  as  the  perturbation  frequency  is  increased. 

The  Velocity  Coupled  Response 

As  stated  previously,  there  have  been  no  experimental  results  for 
the  velocity  coupled  response  of  composite  propellants  reported  in  the 
open  literature.  Thus,  it  was  not  possible  to  attempt  a quantitative 
comparison  of  the  theoretical  velocity  coupled  response  results  with 
experimental  results. 
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The  Oxidizer  Partic'e  Size  Effect 


The  velocity  coupled  response  of  each  pseudopropellant  is  a com- 
plicated function  of  the  pressure  coupled  response  and  the  erosive 
strength  of  the  pseudopropellant.  Thus,  it  is  much  more  difficult  to 
determine  the  effect  of  oxidizer  particle  size  on  the  velocity  coupled 
response  than  in  the  case  of  the  pressure  coupled  response  alone.  The 
effect  of  particle  size  on  the  erosive  burning  characteristics  as  well 
as  the  pressure  coupled  response  of  the  pseudopropellants  were  described 
previously.  Since  the  erosive  effect  is  more  dominant  with  the  large 
oxidizer  particles  and  the  large  oxidizer  particles  are  only  resonant 
at  low  frequencies,  the  velocity  coupled  response  has  the  potential  to 
be  more  dominant  in  propellants  that  contain  large  oxidizer  particles. 
However,  due  to  the  complicated  interaction  between  the  pressure  coupled 
response  and  the  erosive  strength,  a^d  the  fact  that  the  magnitude  of 
the  pressure  coupled  response  of  the  pseudopropellants  is  higher  for  the 
smaller  oxidizer  particles,  the  velocity  coupled  response  is  positive 
only  for  the  intermediate  particle  diameters.  Figure  49  illustrates 
that  as  the  perturbation  frequency  is  increased,  the  magnitude  of  the 
velocity  coupled  response  decreases.  This  is  due  to  the  fact  that  the 
resonant  oxidizer  particle  diameter  decreases  to  diameters  where  the 
erosive  effect  is  negligible. 
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CONCLUSIONS 


Introduction 

The  conclusions  for  the  combustion  of  polydisperse  composite  pro- 
pellants are  divided  into  four  parts.  They  are  the  steady  state  burning 
rate,  the  erosive  burning  rate,  the  pressure  coupled  response  and  the 
velocity  coupled  response. 

The  Steady  State  Burning  Rate 

The  results  of  calculations  with  the  steady  state  PEM  indicate  that 
the  model  can  accurately  predict  the  burning  rate  and  pressure  exponent 
of  HTPB/AP  class  propellants  with  a very  wide  range  of  oxidizer  particle 
sizes. 

Propellants  containing  a very  wide  overall  distribution  of  particle 
sizes  along  with  large  particle  diameters  were  in  the  poorest  agreement 
with  the  experimental  results.  This  was  attributed  to  the  PEM's  lack 
of  consideration  of  interactions  between  the  flames  of  adjacent  pseudo- 
propellants, The  case  for  flame  interactions  was  further  demonstrated 
by  the  value  for  the  diameter  exponent  in  the  equation  which  relates, 
geometrically,  the  amount  of  fuel  associated  with  a particular  oxidizer 
particle  to  the  diameter  of  that  particle,  A value  of  3.0  for  the 
diameter  exponent  was  necessary  to  obtain  agreement  of  the  theoretical 
results  with  the  experimental  results.  Thus,  although  on  a geometric 
basis  small  particles  have  proportionally  more  fuel  associated  with  them 


than  do  the  large  particles  in  an  oxidizer  particle  distribution,  much 
of  that  fuel  is  given  up  to  the  larger  particles  in  the  distribution 
during  the  gas  phase  mixing  process. 

The  effect  of  the  oxidizer  particle  size  on  the  steady  state 
burning  rate  has  long  been  '^ecognized.  However,  the  effect  of  the 
oxidizer  particle  size  distribution  has  not.  The  results  of  calcula- 
tions with  the  PEM  clearly  indicate  that  the  particle  size  distribution 
should  be  considered  during  the  tailoring  of  propellants  to  meet 
specific  burning  rate  and  pressure  exponent  requirements.  As  will  be 
discussed  in  a subsequent  section,  the  advantages  of  employing  the 
particle  size  distribution  in  the  tailoring  process  becomes  even  more 
pronounced  when  the  nonsteady  burning  characteristics  are  also  to  be 
considered  in  the  tailoring  process. 

Some  genera]  conclusions  concerning  the  effect  of  particle  size  on 
the  burning  rate  and  pressure  exponent  follow: 

1.  The  burning  rate  of  small  particles  approaches  an  upper  limit 
determined  by  the  premixed  kinetically  controlled  primary  flame  between 
the  oxidizer  and  the  fuel  binder. 

2.  The  burning  rate  of  large  particles  approaches  a lower  limiting 
value  determined  by  the  kinetically  controlled  AP  monopropellant  flame. 

3.  The  burning  rate  of  the  intermediate  particles  lies  somewhere 
between  the  above  extremes  and  is  determined  by  a competitive  process 
between  the  three  flames  above  a pseudopropellant. 

4.  The  pressure  exponent  approaches  a limiting  value  of  approxi- 
mately half  the  kinetically  controlled  primary  flame  reaction  order  for 
small  particles. 
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5.  The  pressure  exponent  approaches  a limiting  value  of  approxi- 
mately half  the  kinetically  controlled  AP  monopropel lant  flame  reaction 
order  for  large  particles. 

6.  The  pressure  exponent  of  intermediate  oxidizer  particles  is 
determined  by  the  degree  of  dominance  of  diffusion  in  the  controlling 
mechanism  for  combustion.  The  higher  the  degree  of  dominance  of 
diffusion  the  lower  will  be  the  pressure  exponent. 

The  Erosive  Ourning  Rate 

In  the  previous  section  it  was  demonstrated  that  the  erosive  PEM 
predicts  the  following  experimentally  observed  trends  in  erosive  burning: 

1.  Threshold  velocities  are  usually  obsen^ed. 

2.  Slower  burning  propellants  are  more  strongly  affected  by  cross 
flows  than  higher  burning  rate  formulations. 

3.  The  threshold  velocity  decreases  as  pressure  is  increased. 

The  agreement  of  the  theoretically  predicted  trends  with  the  ex- 
perimentally observed  trends  indicates  that  the  proposed  theory  for 
erosive  burning  is  quite  plausible. 

Several  conclusions  concerning  the  physical  mechanism  for  the 
threshold  velocity  can  be  stated: 

1.  At  low  cross  flow  velocities  or  high  burning  rates  the  blowing 
parameter  is  very  high.  Consequently,  the  shear  stress  or  level  of 
turbulence  in  the  flame  zone  is  very  low  and  the  contribution  of  tur- 
bulence to  the  transport  properties  in  the  gas  phase  reaction  zone  is 
negligible. 

2.  At  higher  cross  flow  velocities  or  very  low  burning  rates  the 
blowing  parameter  is  low.  Consequently,  the  shear  stress  or  level  of 
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turbulence  in  the  flame  zone  enhances  the  transport  properties  in  that 
flame  zone. 

3.  The  threshold  velocity  results  since  the  turbulence  cannot 
penetrate  into  the  gas  phase  reaction  zone  (due  to  the  blowing)  utitil 
the  cross  flow  velocity  reaches  a certain  level. 

4.  The  physical  effect  of  burning  rate  on  the  threshold  velocity 
is  due  solely  to  the  effect  of  blowing  on  the  turbulent  boundary  layer. 

5.  The  effect  of  pressure  on  the  threshold  velocity  is  due  to  the 
effect  of  pressure  on  the  turbulent  boundary  layer.  Even  though  the 
blowing  parameter  increases  with  increased  pressure,  the  shear  stress 
in  the  gas  phase  reaction  zone  increases  with  pressure  at  a faster  rate 
than  the  increased  blowing  can  decrease  the  shear  stress. 

6.  Due  to  the  complicated  interaction  between  the  blowing  para- 
meter, the  turbulent  boundary  layer  and  the  burning  rate;  the  erosive 
burning  is  not  an  additive  effect. 

The  erosive  PEM  has  also  made  possible  several  unique  observations 
concerning  the  effect  of  the  oxidizer  particle  size  on  the  erosive 
burning  characteristics  of  the  polydisperse  propellant.  They  are  as 
follows: 

1.  As  the  cross  flow  velocity  is  increased,  it  is  the  large 
particles  v/hich  are  first  affected  as  the  turbulence  extends  into  the 
final  flame  zone  of  the  large  particles. 

2.  Further  increases  in  the  cross  flow  velocity  extends  the  tur- 
bulence into  the  flame  zone  of  smaller  oxidizer  particles. 

3.  The  dominant  mechanism  for  the  combustion  of  the  large  particles 
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shifts  from  the  kinetically  controlled  AP  nionopropellant  flame  to  a 
combination  of  that  flame  and  the  diffusion  controlled  final  flame. 

The  Pressure  Coupled  Response 

The  two  techniques  employed  to  calculate  the  pressure  coupled 
response,  the  ZM  technique  and  the  Cohen-Denison  and  Baum  technique, 
both  predicted  the  correct  order  of  magnitude  of  the  response  in  the 
comparisons  to  the  experimental  data.  For  frequencies  below  4000  hertz 
the  Cohen-Denison  and  Baum  technique  prediction  of  the  magnitude  of  the 
response  was  in  good  agreement  with  the  experimental  results.  The  ZN 
technique,  on  the  other  hand,  underpredicted  the  magnitude  of  the 
response  in  every  case.  The  ZN  technique  did,  however,  predict  the 
approximate  frequency  where  the  absolute  maximum  in  the  response 
function  occurred  in  each  case. 

Multiple  peaks  in  the  response  function  were  predicted  by  both 
methods. 

Several  conclusions  concerning  the  effect  of  oxidizer  particle  size 
on  the  pressure  coupled  response  follow: 

1.  At  low  frequencies  the  large  particles  are  resonant;  however, 
their  response  magnitude  is  low. 

2.  As  frequency  increases  the  resonance  moves  to  smaller  particles. 
Moreover,  the  magnitude  of  the  resonant  response  increases. 

3.  Addition  of  fine  oxidizer  shifts  the  resonant  frequency  of  the 
propellant  to  a higher  frequency  while  increasing  the  magnitude  of  the 
peak.  Conversely,  addition  of  coarse  oxidizer  lowers  the  resonant 
frequency  and  the  response  magnitude. 
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4.  The  shape  of  the  pressure  coupled  response  curve  can  be 
adjusted  so  that  relative  minimums  in  the  response  oc« ur  at  harmonic 
frequencies  in  the  rocket  motor  chamber  simply  by  adjusting  the  oxidi/er 
particle  size  distribution. 

The  Velocity  Coupled  Response 

The  velocity  coupled  response  calculations  showed  that  for  a 
linearized  analysis,  the  velocity  coupled  response  is  a very  complicated 
function  of  the  pressure  coupled  response  and  the  erosive  strength  of  the 
individual  pseudopropellants  in  the  overall  oxidizer  particle  size 
distribution. 

The  conclusions  regarding  the  effect  of  the  oxidizer  particle  size 
on  the  velocity  coupled  response  follow: 

1.  At  low  perturbation  frequencies  the  response  is  dominated  by 
the  large  oxidizer  particles  and  is  negative. 

2.  As  frequency  increases  the  velocity  coupled  response  becomes 
positive  due  to  the  influence  of  the  intermediate  sized  oxidizer 
particles. 

3.  At  high  frequencies  (on  the  order  of  5000  hertz,  but  dependent 
on  the  propellant  burning  rate)  the  velocity  coupled  response  vanishes. 

4.  In  general,  the  velocity  coupled  response  will  be  enhanced 
by  the  addition  of  coarse  oxidizer  and  conversely  the  response  will 
be  depressed  by  the  addition  of  fine  oxidizer. 

The  theoretical  results  for  velocity  coupling  are  suspect  for  two 
reasons: 

1.  A linear  analysis  was  employed  so  that  the  velocity  and  pressure 
coupled  response  could  be  assumed  to  be  independent  of  each  other. 
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2.  The  blocking  effect  due  to  injection  at  the  surface  was.  for 
the  sake  of  simplicity,  ignored  in  the  derivation  of  the  relation  be- 
tween the  velocity  coupled  and  pressure  coupled  responses. 
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